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1.  INTRODUCTION 


The  ultimate  goal  of  this  research  is  to  develop  a  new  non-invasive  imaging  modality  for 
breast  imaging  that  is  capable  of  detecting  small  microcalcifications.  Our  present  goal  is 
to  further  develop  an  experimental  system  based  on  Ultrasound  Stimulated  Acoustic 
Emission  (US  AE)  for  in  vitro  imaging  of  breast  tissue  specimens  and  test  its  performance 
in  detecting  microcalcifications.  The  general  theory  of  USAE,  later  named  vibro- 
acoustography,  is  described  in  [1-4],  and  some  of  its  applications  are  explained  in  [5-7]. 
This  research  includes  two  tasks.  The  first  task  of  this  research  is  focused  on  the 
development  of  the  theory,  system  development,  and  optimization  of  the  method  for 
detection  of  beast  microcalcifications.  The  second  task  is  centered  on  experiments  and 
breast  tissue  imaging.  This  final  report  addresses  the  activities  related  to  both  tasks. 


2.  BODY 

This  research  includes  two  main  tasks:  (1)  System  development,  including  theoretical 
developments,  system  design  and  optimizations.  (2)  Detection  of  microcalcifications  in 
breast  tissue  samples  (experiments).  Here  we  report  activities  on  these  two  tasks. 

2.1.  System  and  Theory  Developments: 

In  this  section,  we  report  activities  relating  to  the  first  task  of  this  project.  These  activities 
include  theoretical  developments,  experimental  system  development  and  optimization. 

2.1.1 .  Development  of  the  vibro-acoustography  imaging  theory: 

We  reported  the  principles  and  theory  of  USAE  (or  vibro-acoustography)  in  [1-4]. 

This  theory  includes  modeling  USAE  imaging  in  four  parts: 

1 .  Generation  of  localized  radiation  force  (vibrating  force)  by  means  of  two 
ultrasound  beams 

2.  A  model  for  object  vibration  in  response  to  the  radiation  force 

3.  A  model  for  the  acoustic  emission  field  in  the  surrounding  medium  due  to  object 
vibration 

4.  A  model  for  image  formation  and  system  point  spread  function  (PSF). 

Based  on  this  theory,  a  USAE  (vibro-acoustography)  system  was  designed  for  imaging 
tissue  samples  [7].  We  have  refined  this  theory  in  several  aspects  during  the  past  3  years. 
Because  the  theory  plays  a  central  role  in  understanding  USAE,  here  we  present  a 
summary  of  the  revised  version  of  USAE  theory. 

Radiation  force  is  generated  by  a  change  in  the  energy  density  of  an  incident  acoustic 
field.  Consider  a  collimated  ultrasound  beam  interacting  with  an  object.  The  radiation 
force  arising  from  this  interaction  has  a  component,  F,  in  the  beam  direction.  This 
component  is  proportional  to  the  time  average  energy  density  of  the  incident  wave  <E> 
and  the  projected  area  of  the  object,  S,  as: 
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F  =  drS  <E>  (1) 

where  dr  is  the  drag  coefficient  (or  the  radiation  force  function),  and  is  a  function  of  the 
scattered  and  absorbed  power  by  the  object  [1,2].  For  the  simple  case  of  a  reflecting  plane 
target,  dr  is  proportional  to  the  power  reflection  coefficient.  In  vibro-acoustography  [1, 

2],  this  force  is  used  for  imaging.  This  is  accomplished  by  probing  the  object  point-by- 
point.  This  technique  ideally  requires  the  stress  field  to  be  confined  to  a  point,  while  its 
amplitude  oscillates  at  selected  frequencies.  To  generate  a  localized  oscillatory  stress 
field,  two  intersecting  CW  focused  ultrasound  beams  of  different  frequencies  are  used.  It 
is  only  in  the  intersection  region  that  the  ultrasound  field  energy  density  is  sinusoidally 
modulated,  and  hence,  the  field  can  generate  an  oscillatory  radiation  force  by  interacting 
with  the  object.  For  this  purpose,  we  use  a  two-element  focused  annular  array,  consisting 
of  a  ring  and  a  concentric  disc  element,  both  focused  at  the  same  focal  point.  We  refer  to 
this  arrangement  as  “confocal”  beams. 

We  consider  a  two-element  spherically  focused  annular  array,  consisting  of  a 
central  disc  with  the  radius  of  au  and  an  outer  ring  with  the  inner  radius  of  a2  and  the 
outer  radius  of  ci2.  The  common  focal  length  of  the  elements  is  zo-  We  also  assume  that 
the  elements  are  excited  by  two  CW  signals  at  frequencies  cox  and  co2  =  cox  +  Act) .  Let  us 
assume  that  the  beams  are  propagating  in  the  z  direction  with  the  joint  focal  point  at  z-0. 
The  resultant  field  at  a  point  defined  by  vector  7  on  the  z-0  plane  may  be  written  as: 

p(7,  t)  =  Pl  (7)  cos[t»,t  +  yzx  (r )]  +  P2  (7)  cos  [co2t  +  if/2  (r )]  (2) 

where  fJ(F)  and  P2(F)  are  the  pressure  field  radial  profiles  across  the  focal  plane  with 
\j/x (7)  and  y/2(r)  being  the  associated  phase  functions.  It  is  noted  from  the  diffraction 
theory  that  Px(7)  and  P2(F )  are  determined  by  the  two-dimensional  (2D)  Fourier 
transforms  of  the  corresponding  transducer  aperture  [2,  3].  For  the  confocal  transducer 
described  above,  these  amplitude  functions  can  be  written  as  [2]: 

Px  (7)  =  pcU0l  y1 ~ jinc  (3) 

AjZ0  AjZ0  j 

and 


where  /l(  =  2k j (Oi  for  i=l,  2,  are  ultrasound  wavelengths  at  the  i-th  transducer  element, 
Uoi,  i-1,  2,  is  the  velocity  amplitude  at  the  corresponding  element,  and  r  =  |F| . 

Also,  jinc(X)  =  7,  ( X)/kX  ,  where  Ji(X)  is  the  first  order  Bessel  function  of  the  first  kind. 
For  well-focused  beams,  Pt(7)  and  P2(F)  diminish  quickly  away  from  the  origin.  The 
phase  functions,  y/i (F)  =  -Kr2  /  A*z0 ,  for  i-1,  2,  are  conveniently  set  to  be  zero  at  the 


origin. 

It  can  be  shown  that  the  short  time  average  of  the  acoustic  energy  density  in  the 
intersection  region  has  slow  variations  at  frequency  A  co  about  its  long  time  average 
(mean).  Denoting  this  low  frequency  component  by  eAo)(7,t ) ,  we  can  write: 


e^(r,t)=  ^(r)^(r)cos[A6)t+A^(F)]  (5) 

pc 


where  p  is  the  density,  c  is  the  sound  speed,  and 

Ays  =  iff 2  (f)  -  ty\  (r)  =  -r1  Aco!  2cz0  (6) 

Now,  consider  a  planar  target  on  the  focal  plane.  Referring  to  (1)  and  considering  that  the 
average  energy  density  is  position  and  time  dependent  in  this  case  we  may  replace  <E> 
with  eAo)(7,t) .  The  normal  component  of  the  time  varying  force  on  a  unit  area  element 
dS=l  at  7  with  a  unit  local  drag  coefficient  d,(ro)=l  is: 

A =  /‘(r)4— cos[A6;t+A^(r)]  (7) 

pc 

This  function  also  represents  the  distribution  of  force,  or  the  stress  field,  on  the  focal 
plane.  It  is  convenient  to  define  the  complex  amplitude  of  the  stress  field  as: 

FAco{7)  =P'(r)Pfr)exV[jAy,(7)\  (8) 

pc 

When  an  oscillating  force  is  applied  to  the  object,  the  object  vibrates  at  the  frequency  of 
the  stress  field.  Let  us  assume  an  object  on  the  focal  plane.  (To  simplify  the  solution,  here 
we  consider  only  two-dimensional  flat  objects.  In  reality,  the  stress  field  may  be  exerted 
on  a  volume  of  the  object.)  The  radiation  force  on  an  area  element  dS  of  this  object  at 
point  7  is  given  by  FA(0(7)dr (7)dS  . 

Object  vibration  due  to  this  excitation  is  a  function  of  object  size,  shape,  and  mechanical 
impedance.  The  relationship  between  the  applied  force  and  object  motion  may  be  defined 
in  terms  of  a  function  called  the  acoustic  outflow,  which  describes  the  volume  of  the 
surrounding  fluid  displaced  per  unit  time  because  of  object  vibrations.  An  example  of  the 
acoustic  outflow  function  for  a  circular  disc  object  is  given  in  [2],  Here  we  denote  the 
acoustic  outflow  of  the  object  due  to  a  force  at  frequency  Aco applied  on  an  area  element 
dS  at  7  by  QAa(7) . 

Vibration  of  the  particles  in  the  surrounding  medium  produces  an  acoustic  field  called  the 
“acoustic  emission.”  The  characteristic  of  the  acoustic  emission  field  is  partly  a  function 
of  the  mechanical  properties  of  the  vibrating  object  and  partly  a  function  of  the 
propagation  medium.  We  may  write  the  incremental  acoustic  emission  pressure  field 
observed  at  a  point  located  at  l  in  the  surrounding  medium  as: 

dP^  =  pc2M6.(,T)QiJ7)dr(T)F,.(r)dS  (9) 

where  M  (l )  is  the  transfer  function  of  the  medium,  describing  the  effect  of  medium  on 

acoustic  wave  propagating  from  the  object  to  the  observation  point  l  .  The  subscript  Aco 
indicates  the  frequency  of  the  acoustic  emission  field. 

To  calculate  the  total  acoustic  emission  filed,  we  must  integrate  dPAco  over  S.  Before 
integrating  the  acoustic  emission  field,  we  may  take  a  few  steps  to  simplify  our  treatment 
of  the  problem. 

In  many  practical  cases,  the  wavelength  of  this  acoustic  field,  A/l  =  2k c  /  Aco ,  is  much 
larger  than  the  effective  dimension  of  the  vibrating  part  of  the  object.  Hence,  the  object 
approximately  acts  as  a  dipole  or  an  almost  omni-directional  source,  and 
thus  MAoj(1  )may  be  considered  independent  of  the  object. 


9 


We  also  note  that  the  only  terms  in  (9)  that  depend  on  the  object  are  dr(r)  and  QAo)(r) . 
Hence,  we  may  define  the  object  function  as 

GAc0(7)  =  QAa>(7)dr(7)  (10) 

This  function  defines  object  characteristics  at  frequency  Act). 

Now,  to  generalize  the  problem,  we  may  assume  that  the  focus  of  the  ultrasound  beam  is 
moved  from  the  origin  to  an  arbitrary  location  F  on  the  focal  plane  (the  beam  axis  still 
being  parallel  to  the  z-axis).  To  find  the  complex  amplitude  of  the  total  acoustic  emission 
field,  we  now  integrate  dPAa)  over  the  object: 

Pi.Cr)  =  pc‘Ul.tnjGt.<f)FlM(f-r)<tS  (11) 

S 

where  7  is  the  integration  variable.  The  integral  in  the  above  equation  constitutes  a  two- 
dimensional  convolution  of  the  object  function  and  the  stress  field,  and  may  be 
summarized  as: 

PA0)(r)  =  pc2MAjT)GAo)(7)  *  FAco(j)  (12) 

where  *  represents  the  convolution  operation. 

To  produce  an  image,  we  move  (scan)  the  focal  point  of  the  ultrasound  beam  across  the 
object  on  the  focal  plane  and  record  the  amplitude  of  the  acoustic  emission  field  at  each 
beam  position.  The  image  is  formed  by  modulating  the  brightness  of  each  point  on  the 
image  plane  proportional  to  PAco(r)  (or  its  magnitude)  recorded  from  the  corresponding 
point  on  the  object. 


Commonly,  an  imaging  system  is  characterized  by  its  point-spread-function  (PSF),  which 
is  defined  as  the  image  of  a  point  object.  To  determine  the  PSF  of  the  system,  we 
consider  a  point  object  defined  as: 


GAoj(T)  =  S(7) 


(13) 


The  resulting  image  according  to  (12)  is  PAcop{r)  =  pc2MA(0{l  )FAa(r) ,  where  the 

subscript  p  is  added  to  indicate  the  point  object.  It  is  convenient  to  define  the  normalized 
PSF  of  the  system  as: 


HAa>(7)  =  PAap(7)/PAo)p(0) 


A  cop ' 


=  FAJ7)/FAC0(  0) 


(14) 


Referring  to  (3  and  4)  we  can  write: 


Ha„(7)=- 


1 


Cln 


n 


Jinc 


ra{ 

\z0 


( 


a2  jinc 


ra2 


>2 


a2  jinc 


ra2 


exp[yAyr(r)]  (15) 


The  above  equation  indicates  that  the  PSF  is  a  complex  function.  However,  for  a 
reasonably  focused  transducer,  |//Aa)(r)|  is  significant  only  within  a  small  radius.  Within 

this  radius  and  for  Ary  □  ry, ,  i=l,2,  we  may  assume  Ay/ (7)  ~  0  as  one  may  conclude 
from  (6).  Thus,  the  PSF  may  be  considered  as  a  real  function  for  practical  purposes.  For  a 
given  object  GAa(7 ) ,  the  image  may  be  written  as: 

PAco(7)  ~  KGAoj(7)  *  HAoj(7)  (16) 
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where  K  is  a  constant  multiplier. 

Application  of  the  theory  to  microcalcification  imaging: 

Breast  microcalcifications  are  normally  in  sub  millimeter  ranges,  and  they  are  often 
smaller  than  the  ultrasound  beam  width.  Hence  we  may  consider  them  as  “point  objects”. 
Therefore,  mathematically  we  may  model  a  microcalcification  as: 

GAm(r)  ~  AS(r)  (17) 

Were  A  is  a  constant  multiplier.  Referring  to  Eq.  (16)  the  image  of  a  small 
microcalcification  may  be  written  as: 

PAoj(T)~AKHAo,(r )  (19) 

Therefore,  image  of  small  microcalcifications  would  be  displayed  as  the  PSF  of  the 
system.  Typical  PSF  image  of  a  USAE  system  is  presented  in  [2,3]. 

2.1.2.  Method  of  vibrating  spheres: 

In  an  effort  for  better  understanding  of  the  dynamic  behavior  of  microcalcifications  in 
response  to  the  radiation  force  of  ultrasound,  we  have  modeled  microcalcification  as  an 
elastic  sphere  in  a  fluid  medium.  The  sphere  is  exposed  to  two  plane  ultrasound  waves  at 
slightly  different  frequencies.  The  vibration  velocity  of  the  sphere  and  the  acoustic 
emission  field  has  been  calculated  in  terms  of  sphere  size  and  elastic  properties  [8-11]. 

An  interesting  result  of  this  theory  is  that  under  some  conditions,  the  vibrating  sphere 
exhibits  a  resonance  effect.  The  resonance  frequency  of  the  sphere  is  a  function  of  the 
sphere  parameters,  such  as  its  size  and  elasticity,  as  well  as  the  mechanical  properties  of 
the  surround  medium.  In  particular,  it  can  be  shown  that  knowing  the  sphere  parameters, 
the  sheer  modulus  of  the  surrounding  medium  can  be  evaluated  from  the  frequency 
response  of  the  vibrating  sphere.  This  phenomenon  is  illustrated  in  Figure  1  (see 
appendix  1).  This  figure  shows  that  the  frequency  response  of  the  sphere  exhibits  a  clear 
peak,  and  the  resonance  frequency  of  the  sphere  increases  as  the  sheer  elasticity  of  the 
surrounding  medium  (gel)  increases. 

To  validate  this  finding,  experiments  were  conducted  on  steel  spheres  embedded  in  gel 
medium,  and  the  velocity  of  the  spheres  was  measured  directly  by  a  laser  interferometer 
system  (vibrometer).  We  also  measured  the  acoustic  emission  field  by  a  calibrated 
hydrophone.  This  arrangement  is  used  to  model  vibration  of  a  microcalcification  in  breast 
tissue.  The  experimental  result  was  used  to  calculate  the  sheer  modulus  of  the  gel.  The 
resulting  values  of  the  sheer  modulus  were  then  validated  by  comparing  with 
independently  measured  values.  For  more  detail,  refer  to  reference  [8-11]. 

The  above  method  is  significant  in  quantitative  analysis  of  microcalcifications  using 
USAE  in  the  following  two  ways: 

1 .  Knowing  the  parameters  of  the  soft  tissue,  one  can  use  this  method  to  evaluate  the 
size  of  the  microcalcification. 

2.  Knowing  the  parameters  of  the  microcalcification,  sheer  elasticity  of  the 
surrounding  tissue  can  be  remotely  evaluated.  This  information  is  especially 
important  in  tissue  characterization  and  diagnosing  lesions  in  breast. 
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The  above  results  have  opened  a  new  way  for  breast  tissue  characterization. 


2.1.3.  System  Improvements: 

Noise  Reduction: 

The  acoustic  emission  from  microcalcifications  is  normally  very  small  and  often 
comparable  to  the  noise  level.  To  improve  the  signal  quality,  we  reduced  system  noise  by 
various  means.  These  methods  include: 

1 .  Quite  room:  This  room  is  designed  with  walls  and  ceiling  covered  by  sound 
absorbing  material  to  reduce  the  ambient  noise.  We  used  this  room  to  conduct  our 
tissue  experiments. 

2.  Air  table:  To  reduce  transmission  of  structural  vibrations  to  the  experiment  setup,  we 
started  to  use  an  air  table  to  support  the  water  tank.  This  air  table  provides  vertical 
isolation  efficiency  90-99%  at  10Hz.  It  has  improved  the  signal-to-noise  ratio 
significantly. 

3.  Modified  water  tank  design:  To  reduce  multipath  and  reverberation  effects  we 
designed  a  new  tank.  The  walls  of  the  new  tank  is  covered  with  sound  absorbing 
layers,  and  the  bottom  is  covered  by  a  sound  absorbing  layer  positioned  at  15  degree 
angle  to  reduce  reverberation.  We  also  added  sound  absorbing  wedges  in  the  tank  to 
reduce  the  possibility  of  tank  resonance. 

4.  To  reduce  the  noise  from  the  scanning  motors,  we  used  a  customized  scanning  system 
based  on  analog  servomotors  instead  of  previous  digital  step  motors.  The  new  analog 
motors  produce  much  less  noise. 

High-speed  Scanning: 

In  order  to  speed  up  the  scanning  process,  we  must  scan  the  object  in  a  continuous 
motion  while  the  data  are  collected.  The  problem  in  this  method  is  that  the  audio  noise 
generated  by  the  scanning  motors  interferes  with  the  acoustic  emission  signal  from  the 
object,  which  dramatically  reduces  signal  quality.  Therefore,  one  must  use  a  method  that 
stops  the  noise  from  the  scanning  motors  while  they  are  running.  One  way  to  do  this  is  by 
increasing  the  acoustic  emission  frequency  beyond  the  spectral  bandwidth  of  the  scanner¬ 
generated  noise,  such  that  the  noise  can  be  filtered  out  without  losing  the  acoustic 
emission  signal. 

Based  on  this  strategy,  we  have  performed  high  speed  scanning  at  30 — 60  kHz.  A 
programmable  filter  was  used  to  remove  the  scanner  noise.  The  noise  bandwidth  in  our 
system  is  about  20  kHz.  In  this  method,  the  transmitting  transducer  continuously 
transmits  the  ultrasound  while  the  scanner  scans  the  object  without  stopping,  and  at  the 
same  time,  the  hydrophone  continuously  receives  the  resulting  acoustic  emission  signal. 
Using  this  method  an  object  can  be  scanned  in  a  fraction  of  the  time  that  would  be 
required  for  scanning  by  the  conventional  point-by-point  scanning.  The  resulting  images 
have  high  signal-to-noise  ratio.  The  only  drawback  of  this  method  is  that  it  cannot  be 
used  at  acoustic  emission  frequencies  lower  than  20  kHz. 
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Tone  Burst  Scanning: 

Tone  burst  scanning  is  a  method  we  have  used  to  avoid  the  ultrasound  standing  wave 
effect  and  acoustic  emission  multipath  effect  in  the  water  tank.  Ultrasound  standing  wave 
effect  is  caused  by  multiple  reflection  of  the  sound  between  a  highly  reflective  object  and 
the  transducer.  Acoustic  emission  multipath  effect  is  caused  by  multiple  reflection  of  the 
acoustic  emission  by  the  walls  or  other  reflectors  within  the  tank.  Both  these  phenomena 
reduce  the  quality  of  the  final  image  and  can  introduce  false  patterns  on  the  image.  In  this 
method,  we  excite  each  transducer  element  with  tone  bursts  of  continuous  waves.  The 
length  of  the  bursts  is  chosen  to  be  at  least  several  times  the  acoustic  emission  period,  but 
short  enough  to  avoid  any  overlapping  with  its  multiple  reflections  from  the  nearby 
structures  such  as  the  tank  walls. 

There  are  two  major  drawbacks  for  this  method.  The  first  problem  is  limitation  on 
acoustic  emission  frequency.  To  satisfy  the  above  conditions  we  often  have  to  use  a  high 
acoustic  emission  frequency.  In  our  experimental  setup,  the  acoustic  emission  frequency 
needs  to  be  about  40  kHz  or  more.  The  second  problem  is  production  of  unwanted 
acoustic  noise.  Applying  a  sudden  burst  of  energy  to  a  transducer  usually  produces  a 
transient  vibration  of  the  transducer  case.  This  vibration  introduces  a  wideband  noise  in 
the  entire  experimental  setup.  To  reduce  this  effect,  we  can  use  a  gradual  ascending  and 
descending  tone  burst,  which  is  described  in  the  next  section. 


Gating  Tone  Bursts: 

A  "Soft  Gate"  was  constructed  using  a  raised  cosine  shaping  circuit  to 
Reduce  rise  time  and  fall  time  of  tone  burst  edges.  This  approach  greatly  reduces  the 
shock  and  the  acoustic  noise  from  transducer  element,  thus  improving  the  detection 
threshold  of  the  system.  To  test  this  method,  we  scanned  the  RTV  glass-bead  phantom 
with  and  without  the  soft  gate.  Without  this  gate,  the  noise  from  the  transducer  shock 
overwhelmed  the  weak  acoustic  emission  from  the  bead.  The  image  quality  was  greatly 
improved  and  the  bead  became  readily  detectable  after  using  the  soft  gate.  Experiments 
on  the  RTV  phantom  with  holes  also  showed  the  same  improvement. 


Receiver  Lock-in  Amplifier: 

Our  previous  receiver  included  a  per-amplifier  followed  by  a  narrow-band  programmable 
filter.  This  configuration  operates  well  as  long  as  the  signal  level  at  the  hydrophone  is 
high  enough  for  subsequent  amplification  and  detection.  However,  when  imaging  small 
microcalcifications,  the  signal  level  can  be  very  low  in  comparison  to  the  ambient  noise 
within  the  pass-band  of  the  narrow-band  filter.  There  are  several  sources  of  noise  in  the 
system,  including  airborne  noise  in  the  room,  structural  vibrations,  acoustic  noise  from 
the  scanning  motors,  and  the  noise  from  equipment  fans.  There  is  also  some  electrical 
noise  in  the  system. 

To  improve  the  signal-to-noise  ratio,  we  modified  the  receiver  by  replacing  the 
preamplifier  and  the  programmable  filter  with  a  lock-in  amplifier  (EG&G  Instruments, 
Model  7264  DSP).  The  lock-in  amplifier  uses  a  reference  signal  to  set  the  lock  frequency. 
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The  frequency  of  the  reference  signal  must  be  exactly  equal  to  the  frequency  of  the  input 
signal  to  the  amplifier  (hydrophone  output).  This  system  locks  to  the  frequency  of  the 
reference  signal  (Am)  and  tracks  the  input  signal  by  an  active  phase-lock  loop.  In  effect, 
the  lock-in  amplifier  acts  as  a  very  narrow  band  filter  with  a  pass  band  of  only  a  few  Hz 
centered  at  the  difference  frequency  Am.  This  system  allows  the  acoustic  emission  signal 
to  pass  and  rejects  the  noise  outside  the  pass  band.  The  reference  signal  (at  Ago)  is 
provided  by  down-mixing  the  signals  from  the  two  RF  generators  driving  the  two 
elements  of  the  transducer. 

Improved  System  Control: 

To  improve  the  scanning  process  a  GPIB  control  card  has  been  added  to  the  scanner 
workstation  (Sun  Sparc  Station  5).  This  card  allows  all  sub-systems,  such  as  RF 
generators,  lock-in  amplifier,  digitizer,  and  the  new  scanner  system  with  servomotors,  to 
be  controlled  by  software.  This  system  operates  under  LabView  software. 

Direct  Displacement  Measurement  by  Laser  Vibrometer: 

We  started  using  a  laser  vibrometer  for  direct  measurement  of  velocity  and  displacement 
in  the  order  of  nano-meters  at  frequencies  up  to  250  kHz  range.  This  system  allows  us  to 
directly  measure  the  vibration  that  we  introduce  in  an  object  by  the  radiation  force  of 
ultrasound.  Hence,  we  can  validate  our  acoustic  measurements  by  direct  optical 
measurement. 

Direct  Displacement  Measurement  by  Scanning  Vibrometer 

The  main  function  of  US  AE  is  to  vibrate  the  object  by  the  radiation  force  of  ultrasound.  It 
is  therefore  a  fundamental  issue  to  understand  object  vibration  and  its  role  in  generating 
the  acoustic  emission  field.  In  order  to  analyze  object  vibration  and  verify  our  acoustic 
measurements  we  need  to  measure  such  vibrations  directly.  For  this  purpose  we  have 
acquired  a  scanning  vibrometer  (Polytec  VibraScan  PSV-300F)  that  enables  us  to 
measure  object  vibrations  in  two  dimensions,  in  Angstrom  range,  and  produce  an  image 
(or  a  movie)  that  shows  displacement  (or  particle  velocity)  at  every  point  on  the  objet. 
(This  instrument  is  in  addition  to  the  laser  vibrometer  mentioned  above,  which  measures 
vibration  at  only  one  point  at  a  time.)  This  tool  will  help  us  to  further  develop  USAE  for 
breast  imaging.  The  funds  or  this  instrument  is  provided  by  a  Major  Research 
Instrumentation  grant  from  the  National  Science  Foundation. 


Anechoic  Chamber: 

Because  USAE  experiments  are  performed  at  low  acoustic  intensities,  it  is  essential  to 
keep  the  environmental  noise  and  structural  vibration  as  low  as  possible.  For  this 
purpose,  we  have  designed  a  12’xl0’x7’  anechoic  chamber  with  a  150  Hz  cutoff 
frequency.  This  room  includes  acoustically  insulated  walls  and  floor,  with  the  interior 
covered  with  absorbing  wedges.  This  room  significantly  reduces  the  noise  level  and 
improve  the  quality  of  data  in  our  experiments.  Funding  for  the  anechoic  chamber  is 
provided  by  a  Major  Research  Instrumentation  grant  from  the  National  Science 
Foundation.  Mayo  Clinic  provides  the  space  for  this  room. 


14 


Software  Development: 


1.  Data  Acquisition:  The  software  for  controlling  the  new  scanning  system  and  the  data 
acquisition  system  was  developed. 

2.  Interactive  Imaging  and  Object  Analysis:  New  software  has  been  developed  for 
interactive  scanning  and  data  collection.  An  important  feature  of  the  new  software  is 
its  capability  to  produce  the  image  of  the  object  as  it  is  scanned,  and  allows  the  user 
to  select  a  point  on  the  image  and  reposition  the  transducer  beam  to  the  selected 
points  in  the  object.  Once  at  the  new  position,  the  system  can  collect  the  acoustic 
emission  signal  from  that  particular  point  of  the  object.  Another  important  feature  of 
this  software  is  its  capability  in  measuring  the  frequency  response  of  a  selected  point 
in  the  object.  This  feature  allows  the  user  to  point  the  ultrasound  beam  at  a  selected 
point  in  the  object  and  sweep  the  vibration  frequency  Ao)  while  the  system  records  the 
acoustic  emission  signal.  By  employing  a  built-in  fast  Fourier  transform  algorithm, 
this  system  then  calculates  and  displays  the  frequency  response  (spectrum)  of  the 
object.  This  feature  is  particularly  useful  in  imaging  and  analyzing  the  frequency 
response  of  microcalcifications.  Employing  the  method  of  vibrating  spheres  described 
earlier,  we  may  be  able  to  use  the  frequency  domain  signatures  of  microcalcifications 
to  characterize  the  soft  tissue  or  microcalcifications  in  breast.  This  would  open  the 
way  for  using  USAE  as  a  “non-invasive  biopsy”  tool. 


2.2.  Tissue  Experiments 

The  hypothesis  is  that  ultrasound  stimulated  acoustic  emission  can  be  used  to  image 
tissue  and  detect  microcalcifications.  To  test  this  hypothesis,  we  scanned  75  excised 
human  breast  tissue  and  compared  the  results  with  the  x-ray  mammography  and 
histology. 

2.2.1.  Procedure 

Experiments  were  conducted  on  human  breast  tissue  samples  in  the  following  manner: 

(a)  Initial  tissue  selection 

(b)  Initial  tissue  preparation 

(c)  Screening  mammography 

(d)  Reading  mammography 

(e)  Final  tissue  selection 

(f)  Mounting  on  Scanning  Bracket 

(g)  Base  Mammography 

(h)  USAE  Scanning 
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(i)  Image  evaluation  (USAE  vs.  x-ray) 

0)  Histological  study 

2.2.2.  Results: 

Experiments  were  conducted  on  breast  tissue  samples  collected  from  94  patients.  We 
mounted  219  tissue  samples  in  total.  From  this  collection,  207  tissue  samples  were  x- 
rayed  using  high-resolution  mammography.  12  tissue  samples  were  excluded  from  the 
study  because  they  were  not  usable  (lost  their  water  content  or  became  too  old  for 
imaging).  Based  on  the  x-ray  images,  we  found  microcalcifications  in  75  from  207  tissue 
samples,  and  another  22  samples  were  suspicious  of  having  some  microcalcifications.  We 
conducted  USAE  (vibro-acoustography)  imaging  of  those  tissue  samples  that  were 
microcalcifications  positive  in  their  x-rays.  Out  of  the  75  sample,  one  was  excluded  from 
the  study  because  of  the  tissue  was  damaged  in  the  process  of  mounting.  In  the  USAE 
images,  we  found  44  samples  to  be  microcalcification  positive  (a  ratio  of  0.6),  out  of 
which  42  were  also  histologically  positive.  Microcalcifications  as  small  as  1 10  pm  in 
diameter  could  be  identified  in  USAE  images.  We  also  studied  appearance  of  soft  tissue 
in  USAE  images.  Tissue  structure  appearance  in  the  USAE  images  were  rated  perfect  or 
fair,  as  compared  to  the  x-ray.  We  found  49  perfect  and  25  fair  tissue  structures  in  the 
total  74  samples  imaged.  More  details  on  experiments  can  be  found  in  [5-7]. 


Discussion: 

The  theory  presented  here  clarifies  the  fundamental  relationships  between  the  object 
parameters  and  its  USAE  image.  This  theory  indicates  that  the  acoustic  emission  from  an 
object  is  related  to  its  mechanical  admittance  and  the  size  of  the  object. 

In  the  case  of  small  microcalcifications,  detection  of  the  acoustic  emission  is  often  a 
challenge.  To  improve  the  detectability  by  USAE,  the  acoustic  noise  in  the  medium  must 
be  minimized.  The  acoustic  noise  includes  the  ambient  noise,  noise  from  the  scanner, 
vibrations  produced  by  the  transducer  body,  etc.  The  improvements  that  we  have  made 
in  the  laboratory  setup,  including  using  the  quiet  room,  airtable,  absorbing  material,  soft 
gate,  etc.,  have  enabled  us  to  considerably  increase  the  signal-to-noise  ratio  of  our  data. 
Now,  with  our  low-noise  setup,  we  can  detect  smaller  glass  beads  in  our  control  test 
phantoms.  Using  the  tone  burst  method  has  enabled  us  to  eliminate  the  standing  wave 
problem  at  the  expense  of  working  at  a  higher  acoustic  emission  frequency. 

Driving  a  particle  at  its  resonance  could  provide  useful  information  about  the  particle  and 
its  surrounding  medium.  In  addition,  it  can  enhance  the  detectability  and  selectivity.  That 
is  because  at  resonance  particle  motion  is  at  its  maximum,  hence  produces  a  stronger 
acoustic  emission.  Also,  because  the  resonance  frequency  of  a  particle  depends  on  its  size 
and  mass,  this  method  can  be  used  to  selectively  detect  particles  of  particular  size  or 
mass.  Based  on  this  concept  we  developed  the  theory  of  vibrating  spheres  in  a 
viscoelastic  medium.  An  interesting  result  of  this  method  is  that  by  measuring  the 
resonance  frequency  of  such  sphere  it  is  possible  to  calculate  the  sheer  elasticity  of  the 
medium. 
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Experimental  USAE  images  presented  in  this  report  demonstrate  two  important  facts:  (1) 
USAE  imaging  method  is  capable  of  detecting  small  microcalcifications,  as  well  as 
calcifications  in  vessels,  in  breast  tissue;  (2)  Microcalcifications  can  be  delineated  form 
other  tissue  structures.  The  results  shown  here  also  indicate  that  such  images  have  high 
spatial  resolution,  good  contrast,  and  high  signal  to  noise  ratio.  Microcalcifications  shown 
in  our  results  are  only  a  few  hundred  microns  in  diameter,  and  the  calcifications  in  the 
vessel  wall  were  not  more  than  100  microns  thick. 

A  comparison  between  the  USAE  image  (Figure  4  of  reference  [7])  and  the 
corresponding  mammography  reveals  an  interesting  capability  of  USAE  method.  High 
tissue  density  in  the  region  surrounding  the  microcalcification  absorbs  a  great  deal  of  x- 
ray  energy,  hence  the  x-ray  image  shows  very  low  contrast  within  this  region,  and  as  a 
result,  microcalcifications  are  barely  visible.  However,  the  USAE  image  delineates  these 
calcifications  with  high  contrast.  This  is  because  the  acoustic  emission  is  sensitive  to  the 
elastic  properties  of  the  object. 

The  histology  of  the  tissue  samples  indicates  that  USAE  can  image  breast 
microcalcifications  in  a  wide  range  of  pathological  conditions  (refer  to  2001  Annual 
Report).  This  improves  our  prospects  for  future  use  of  USAE  for  in-vivo  breast  imaging. 
In  most  cases,  tissue  inhomogenieties  seem  not  to  interfere  with  detection  of 
microcalcifications. 

Acoustic  images  reported  here  were  obtained  at  relatively  high  vibration  frequencies  (20- 
-40  kHz).  At  these  frequencies  we  are  able  to  acquire  the  data  at  much  higher  speed, 
reducing  the  acquisition  time  to  a  few  minutes  per  image  (for  more  details,  refer  to  the 
Annual  Report  2000).  This  is  because  at  these  frequencies  the  noise  level  from  the 
scanner  mechanism  is  negligible;  hence,  we  can  keep  the  scanning  motor  running  while 
collecting  the  data.  Scanning  speed  is  an  important  factor  in  our  experiments  and 
eventually  in  application  of  USAE  on  humans. 


3.  Key  Research  Accomplishments 

•  Development  of  the  theoretical  model  describing  the  relation  between  the  acoustic 
emission  field  and  the  mechanical  parameters  of  the  object  and  the  surrounding 
medium. 

•  Development  the  theory  describing  the  point  spread  function  of  the  imaging 
system. 

•  Improvement  of  the  sensitivity  of  the  system  to  be  able  to  detect  small  particles  in 
a  viscoelastic  medium. 

•  Identifying  causes  of  image  artifacts,  such  as  ultrasound  standing  wave, 
multipath,  and  reverberation  of  acoustic  emission. 

•  Developing  methods  for  reducing  or  eliminating  acoustic  emission  multipath  and 
reverberation  problems. 
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•  Developing  methods  for  reducing  acoustic  and  vibrational  noise  in  the  system 

•  Demonstrated  that  small  microcalcification  (as  small  as  1 10  microns  in  diameter) 
and  vessel  calcifications  can  be  detected  by  USAE  imaging  method. 

•  Demonstrated  that  in  USAE  images  microcalcification  can  be  delineated  from 
other  tissue  structures. 

•  Demonstrated  that  microcalcification  can  be  imaged  at  high  vibrational 
frequencies. 

•  Demonstrated  that  USAE  can  detect  microcalcifications  in  various  pathological 
conditions. 

•  Developed  the  theory  of  vibrating  sphere  in  a  viscoelastic  medium 

•  Demonstrated  that  a  small  sphere  in  a  viscoelastic  medium  can  exhibit  a 
resonance.  The  resonance  frequency  depends  on  the  parameters  of  the  sphere  and 
the  mechanical  properties  of  the  medium. 

•  Demonstrated  that  sheer  elasticity  of  the  medium  can  be  evaluated  based  on  the 
resonance  of  the  sphere.  This  method  provides  the  possibility  of  characterizing  the 
soft  tissue  surround  the  microcalcification. 
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5.  Conclusions 

Small  microcalcifications  (as  small  as  1 10  microns  in  diameter)  and  vessel  calcifications 
in  breast  tissue  can  be  detected  by  USAE  imaging  method.  In  the  images 
microcalcifications  appear  with  high  contrast  with  respect  to  the  soft  tissue,  hence  can  be 
delineated  from  other  tissue  structures.  Image  quality  and  system  sensitivity  is  improved 
by  using  various  noise  reduction  methods.  It  is  also  shown  that  USAE  may  be  used  to 
determine  the  mechanical  properties  of  the  medium  surrounding  a  small  sphere.  This 
finding  opens  the  way  for  characterization  of  the  lesion  surrounding  a  breast 
microcalcification. 

5.1.  So  What? 

The  body  of  knowledge  collected  in  this  research  may  be  used  to  development  of  a  new 
class  of  non-invasive  imaging  tool  for  breast  imaging  and  visualization  of  breast 
microcalcifications.  The  results  of  the  present  research  on  breast  tissue  samples  warrant 
further  research  to  explore  the  utility  of  USAE  in  detection  of  microcalcifications  and 
lesions  in  human  breast.  This  research  will  be  continued  at  a  more  advanced  level  within 
the  framework  of  two  funded  research  grants.  The  first  grant  on  in-vivo  breast  vibro- 
acoustography  (USAE)  is  funded  by  Komen  Breast  Cancer  Foundation. 

The  second  grant,  which  is  funded  by  NIH,  is  focused  on  development  of  a  vibro- 
acoustography  system  that  will  be  designed  for  in  vivo  applications. 
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Appendix  1 


Appendix  1  (Figure  1) 


Velocity(mm/s)  -  ,u=ucp* 


FIG.  1.  Simulation  of  a  sphere  vibrated  in  gel.  Parameters  for  the  sphere:  radius=0.59 
mm,  density  =  7667 kg/m3.  The  gel  has  a  shear  viscosity  of  0.1  Pa  s,  and  shear 
elasticity  of  1, 3,  10,  and  30  kPa,  respectively. 
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Vibro- Acoustic  Tissue  Mammography 
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Abstract — A  novel  method  for  detection  and  imaging  of  micro¬ 
calcifications  in  breast  tissue  is  presented.  The  method,  called 
vibro- acoustography,  uses  the  radiation  force  of  ultrasound  to 
vibrate  tissue  at  low  (kHz)  frequency  and  utilizes  the  resulting 
response  to  produce  images  that  are  related  to  the  hardness  of  the 
tissue.  The  method  is  tested  on  human  breast  tissues.  The  resulting 
vibro-acoustographic  images  are  in  agreement  with  corresponding 
X-ray  mammography  images  of  the  specimens.  The  existence  of 
microcalcifications  in  locations  indicated  by  vibro-acoustography 
is  confirmed  by  histology.  Microcalcifications  as  small  as  110  fim 
in  diameter  are  detected  by  this  method.  Resulting  vibro-acous¬ 
tographic  images  show  microcalcifications  with  high  contrast 
with  respect  to  the  background  soft  tissue.  Structures  such  as 
dense  sclerotic  tissue  do  not  seem  to  interfere  with  detection  of 
microcalcifications. 

Index  Terms — Breast,  mammography,  microcalcifications, 
ultrasound,  vibro-acoustography. 


I.  Introduction 

X-RAY  mammography  is  currently  the  only  imaging 
modality  clinically  used  for  detection  of  breast  microcal¬ 
cifications.  The  widespread  use  of  screening  mammography 
has  resulted  in  the  increased  detection  of  microcalcifications 
[1].  A  wide  spectrum  of  breast  lesions  is  associated  with 
microcalcifications,  ranging  from  benign  (fibrocystic  changes, 
vascular  changes,  fat  necrosis)  to  malignant  [2].  The  mam- 
mographic  detection  of  microcalcifications  frequently  results 
in  additional  diagnostic  studies,  including:  spot  magnification 
mammographic  views,  stereotaxic  biopsy  and  wire-localized 
biopsies.  Follow-up  of  microcalcifications  that  are  not  biop- 
sied,  or  that  are  benign  by  stereotaxic  biopsy,  may  require 
additional  mammograms,  e.g.,  at  six-month  intervals.  This 
results  in  significant  additional  costs,  discomfort,  anxiety, 
and  some  increase  in  X-ray  exposure  [1],  Recent  changes  in 
American  Cancer  Society  (ACS)  and  National  Cancer  Institute 
(NCI)  screening  guidelines  for  women  in  their  forties  will 
increase  both  the  number  of  mammograms  and  the  number  of 
detected  microcalcifications,  with  resultant  increased  costs  and 
morbidity.  The  ACS  recommends  annual  mammography  for 
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women  beginning  at  age  40.  In  women  at  “high”  risk  for  breast 
cancer  surveillance  may  begin  at  age  35  [3]. 

There  exist  subsets  of  women  for  whom  mammographic 
screening  and  follow-up  is  limited  by  reduced  mammographic 
sensitivity.  Fibrocystic  change  is  the  most  common  benign 
condition  of  the  breast,  with  the  greatest  incidence  in  women 
30-50  years  old  [4].  The  ACS  estimates  that  about  half  of 
all  women  in  America  have  fibrocystic  breasts  and  many  of 
them  have  very  radiologically  dense  breast  tissue.  For  this 
group,  the  failure  rate  to  detect  cancer  on  their  mammogram 
may  be  as  high  as  15%-25%  [5],  [6].  Mammography  is  more 
effective  in  detecting  occult  malignancy  as  age  increases  and 
the  breast  becomes  more  fatty.  It  is  the  density  of  the  breast 
on  film  screen  mammography  that  most  directly  affects  the 
risk  for  false-negative  and  false-positive  interpretation  [7],  [8]. 
Several  factors  affect  breast  density.  Age  is  the  most  significant 
factor,  with  more  premenopausal  women,  especially  in  the 
40-  to  49-year  age  group,  having  breasts  that  are  denser.  For 
premenopausal  women  mammographic  sensitivity  may  be 
increased  by  scheduling  evaluation  the  first  two  weeks  of 
their  cycle  corresponding  with  the  follicular  phase  [9],  [10]. 
Pregnant  or  lactating  women  may  have  extremely  dense  breast 
tissue,  which  may  affect  the  sensitivity  of  mammography; 
also  these  conditions  are  a  relative  contraindication  for  X-ray 
mammography  and,  thus,  limit  its  usefulness.  Some  studies 
have  shown  hormone  replacement  therapy,  commonly  used  in 
postmenopausal  patients,  reduces  the  sensitivity  of  mammo¬ 
graphic  screening  [11]— [13].  The  sensitivity  or  accuracy  of 
film  screen  mammography  is  also  influenced  by  the  experience 
of  the  radiologist,  with  experienced  radiologists  having  the 
highest  sensitivity  in  diagnosing  breast  cancer  [14].  Finally,  the 
ionizing  nature  of  the  X-ray  mammography  limits  its  frequent 
use. 

To  overcome  some  of  these  problems  associated  with  mam¬ 
mography,  a  number  of  new  technologies  are  currently  being 
explored.  Investigators  have  utilized  alternative  breast  imaging 
methods,  such  as  magnetic  resonance  imaging  (MRI)  and 
conventional  ultrasound  for  imaging  breast  microcalcifications. 
Scintimammography  and  contrast-enhanced  MRI  are  not 
sufficiently  sensitive  to  detect  clinically  significant  microcal¬ 
cifications  [15].  Similarly,  conventional  medical  ultrasound 
imaging  is  not  currently  considered  sufficiently  reliable  to 
visualize  breast  microcalcifications  [16]— [18].  Microcalcifi¬ 
cation  visualization  by  ultrasound  is  limited  by  a  number  of 
factors  specific  to  this  technique,  such  as  speckle  noise  [18]. 
Although  ultrasound  can  be  used  as  an  effective  adjunct  to 
mammography  in  detecting  breast  cancer  [19],  this  technique 
will  miss  occult  cancers  where  microcalcification  may  be  the 
only  sign  of  malignancy  [20].  Thus,  a  sensitivity  of  only  63% 
was  found  in  a  study  by  Kasumi  [21].  However,  Yang,  et  al 
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Fig.  I.  Vibro-acoustography  system  diagram.  The  system  includes  a  focused  confocal  transducer,  (details  shown  on  the  left)  consisting  of  a  center  disk  and  an 
outer  ring.  Two  continuous  wave  generators  drive  these  elements  at  slightly  different  frequencies.  The  transducer  is  focused  on  the  object,  with  the  beams  interacting 
at  the  joint  focal  point  to  produce  an  oscillating  radiation  force  on  the  object  at  the  difference  frequency.  This  force  causes  the  object  to  vibrate  and  as  a  result  an 
acoustic  emission  field  is  produced  in  the  surrounding  medium.  This  field  is  detected  by  the  hydrophone  and  filtered  by  a  bandpass  filter  centered  at  the  difference 
frequency.  The  amplitude  of  the  resulting  signal,  detected  by  the  detector,  is  used  to  modulate  the  intensity  of  the  image  at  a  point  corresponding  to  the  position  of 
the  beam  on  the  object.  The  image  is  formed  by  raster  scanning  of  the  object.  The  experiments  take  place  in  a  water  tank  (not  shown)  containing  the  transducer, 
hydrophone  and  the  object.  The  object,  such  as  a  tissue  specimen,  is  secured  on  a  latex  sheet  attached  to  a  bracket  for  the  scanning  process. 


[22]  reported  that  using  high-resolution  ultrasound  and  X-ray 
mammography  as  the  gold  standard  resulted  in  a  sensitivity 
of  95%  and  specificity  of  87%  and  accuracy  of  91%  in  the 
detection  of  microcalcifications  when  these  calcifications  are 
within  a  mass  lesion  [22],  It  should  be  noted  that  the  presence  of 
a  mass  lesion  and  its  appearance  in  the  ultrasound  image  could 
have  been  instrumental  in  the  localization  and  identification  of 
microcalcifications  in  such  studies. 

Microcalcifications  are  mainly  composed  of  hydroxyapatite 
[18],  [23],  which  is  a  very  hard  material  compared  with  breast 
tissues.  Therefore,  an  imaging  modality  that  is  sensitive  to  the 
elastic  properties  of  tissues  would  likely  be  suitable  for  detec¬ 
tion  of  microcalcifications.  This  group  of  imaging  techniques  is 
called  elasticity  imaging  [24]— [37]. 

The  general  approach  in  elasticity  imaging  is  to  measure  the 
response  of  tissue  to  an  excitation  force.  An  interesting  strategy 
for  producing  the  necessary  excitation  is  to  use  the  radiation 
force  of  ultrasound.  This  approach  allows  one  to  generate  the 
force  directly  inside  an  organ.  Several  researchers  have  previ¬ 
ously  investigated  the  use  of  ultrasound  radiation  pressure  for 
tissue  characterization  purposes.  Sugimoto  etal  [38]  presented 
a  method  to  measure  tissue  hardness  by  using  the  radiation  force 
of  a  single  ultrasound  beam.  In  this  method,  the  impulsive  radi¬ 
ation  force  was  used  to  generate  localized  deformation  of  tissue. 
Resulting  transient  deformation  was  measured  as  a  function  of 
time  by  the  ultrasound  Doppler  method.  Nightingale  et  al  [39] 
have  studied  the  radiation  force  of  ultrasound  in  a  method  named 
“remote  palpation”  to  evaluate  lesions  in  breast.  They  also  used 
the  ultrasound  pulse-echo  technique  to  detect  tissue  displace¬ 
ment.  None  of  the  above  methods  [38],  [39]  are  specifically  ap¬ 
plied  to  the  detection  of  breast  microcalcifications.  Walker  [40] 
evaluated  the  deformation  of  soft  tissues,  including  breast,  by 
the  radiation  force  of  ultrasound.  It  should  be  noted  that  motion 
detection  by  ultrasound  techniques  (at  diagnostic  frequencies) 
becomes  technically  difficult  if  the  motion  is  less  than  1  //,m. 
Walker  predicted  that  the  displacement  produced  by  absorption 
of  a  continuous  wave  ultrasound  beam  of  moderate  intensity 


(1  W /  cm2)  in  breast  tissue  is  too  small  (less  than  0.01  //m)  for 
consistent  ultrasound  detection. 

Vibro-acoustography  is  a  new  imaging  method  based  on  the 
radiation  force  of  ultrasound  [41],  [42].  This  method  is  par¬ 
ticularly  useful  in  detecting  hard  inclusions  in  soft  material. 
For  example,  vibro-acoustography  has  been  used  to  image  cal¬ 
cifications  in  human  arteries  [43]-[45],  A  comparative  study 
of  vibro-acoustography  with  other  radiation  force  methods  for 
tissue  elasticity  imaging  is  presented  in  [46],  The  spatial  res¬ 
olution  of  vibro-acoustography  is  in  the  submillimeter  range, 
which  makes  the  technique  suitable  for  imagine  microcalcifica- 
tions  [47],  [48]. 

In  this  paper,  we  present  an  application  of  vibro-acous¬ 
tography  for  detection  of  breast  microcalcifications.  We  may 
call  the  new  method  vibro-acoustic  tissue  mammography.  We 
compare  images  of  human  breast  tissue  specimens  obtained 
by  vibro-acoustography  with  high-resolution  X-ray  tissue 
mammography  and  validate  the  presence  and  position  of 
calcifications  in  the  tissue  histologically. 

II.  Methods  and  Materials 
A.  Theoty  of  Radiation  Force 

Radiation  force  is  generated  by  a  change  in  the  energy  density 
of  an  incident  acoustic  field.  Consider  a  collimated  ultrasound 
beam  interacting  with  an  object.  The  radiation  force  arising  from 
this  interaction  has  a  component,  F,  in  the  beam  direction.  This 
component  is  proportional  to  the  time  average  energy  density  of 
the  incident  wave  ( E )  and  the  projected  area  of  the  object,  5,  as 

F  =  drS(E)  (1) 

where  dr  is  the  drag  coefficient  and  is  a  function  of  the  scat¬ 
tered  and  absorbed  power  by  the  object.  For  the  simple  case 
of  a  reflecting  plane  target,  dr  is  proportional  to  the  power  re¬ 
flection  coefficient.  In  vibro-acoustography  [17],  this  force  is 
used  for  imaging.  This  is  accomplished  by  probing  the  object 
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point-by-point.  This  technique  ideally  requires  the  stress  field 
to  be  confined  to  a  point,  while  its  amplitude  oscillates  at  se¬ 
lected  frequencies. 

To  generate  a  localized  oscillatory  stress  field,  two  inter¬ 
secting  continuous- wave  focused  ultrasound  beams  of  different 
frequencies  are  used.  It  is  only  in  the  intersection  region  that 
the  ultrasound  field  energy  density  is  sinusoidally  modulated 
and  hence,  the  field  can  generate  an  oscillatory  radiation  force 
by  interacting  with  the  object. 

The  ultrasound  beams  can  be  shaped  in  a  variety  of  ways  for 
this  purpose.  An  interesting  configuration  that  results  in  a  radi¬ 
ally  symmetric  modulated  field  is  obtained  when  two  coaxial, 
confocal  transducer  elements  are  used  (Fig.  1).  In  this  case, 
we  consider  a  two-element  spherically  focused  annular  array, 
consisting  of  a  central  disc  with  the  radius  of  a\  and  an  outer 
ring  with  the  inner  radius  of  a'2  and  the  outer  radius  of  (12.  The 
common  focal  length  of  the  elements  is  z$.  We  also  assume  that 
the  elements  are  excited  by  two  CW  signals  at  frequencies  ui 
and  uj2  =  loi  +  A u.  Let  us  assume  that  the  beams  are  propa¬ 
gating  in  the  z  direction  with  the  joint  focal  point  at  z  —  0.  The 
resultant  field  on  the  z  =  0  plane  may  be  written  as 


p(t)  -  Fi(r)  cos[u;it  +  ^i(r)]  +  P2(r)  cos[<j02t  +  4>2 (r)]  (2) 

where  P\  (r)  and  P2  (r)  are  the  pressure  radial  profiles  in  r  direc¬ 
tion  with  ipi (r)  and  ^ 2{r )  being  the  associated  phase  functions 
across  the  focal  plane.  For  the  transducer  described  above,  these 
amplitude  functions  can  be  written  as  [31] 


Pi(r)  —pcUoi  y— '-jinc 

AlZ0 
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(4) 


where  A*  —  27r/a;*  for  i  =  1,2,  are  ultrasound  wave¬ 

lengths  at  the  ith  transducer  element  Uoi,  i  =  1,2,  is  the 
velocity  amplitude  at  the  corresponding  element.  Also, 
jinc(W)  =  Ji(X)/nX,  where  J\(X)  is  the  first-order  Bessel 
function  of  the  first  kind.  For  well-focused  beams,  P\(r)  and 
P2O*)  diminish  quickly  away  from  the  origin. 

It  can  be  shown  that  the  short  time  average  of  the  acoustic 
energy  density  in  the  intersection  region  has  slow  variations  at 
frequency  A uj  about  its  long  time  average  (mean).  Denoting  this 
low-frequency  component  by  ca^  (t) ,  we  can  write 


(t)  =  cos[Au;t  +  A'0(ro)]  (5) 

pc 

where  p  is  the  density,  c  is  the  sound  speed  and  A ^  =  ip2(r)  - 
ipi  (r).  Now,  consider  a  planar  target  on  the  focal  plane.  Refer¬ 
ring  to  (1)  and  considering  that  the  average  energy  density  is 
position  and  time  dependent  in  this  case,  the  normal  component 


of  the  time- varying  force  on  an  area  element  dS  at  r  =  ro  with 
a  local  drag  coefficient  dr(r0)  is 

fAu(ro,t)  =  •Pl(r°)^>2(ro)  cos[A wt  +  Aip(r0)]dr(r0)dS. 
pc 

(6) 

This  function  also  represents  the  distribution  of  force,  or  the 
stress  field,  on  the  focal  plane.  The  total  force  on  the  object  can 
be  found  by  integrating  the  above  incremental  force  over  the 
projection  area  S  on  the  focal  plane 

fAoj(ro ,  t )  =  |Faw|  cos(Au;*  +  AV>)  (7) 

where  FAuj  is  the  complex  amplitude  of  the  total  force  and  A$ 
is  the  associated  phase.  For  a  well-focused  beam  S  is  very  small 
and  hence,  the  force  can  be  thought  of  as  an  oscillating  point-like 
force  applied  to  the  object  at  the  origin. 


B.  Vibro-Acoustography  System 

To  explain  the  imaging  method,  we  consider  an  oscillating 
point  force,  FAoj,  applied  to  a  point  in  the  object.  This  force 
vibrates  the  object.  Object  vibrations  result  in  the  emission  of 
an  acoustic  field  in  the  surround  energy  that  can  be  detected  by 
a  microphone  (or  hydrophone  in  water).  The  system  diagram  is 
shown  in  Fig.  1.  The  complex  amplitude  of  acoustic  emission 
pressure  field,  Pau,  can  be  written  as  [17] 

Paoj  =  PC2Hau(1)Qau>Fau>  (8) 


where  Qaoj  is  a  complex  function  representing  the  mechanical 
frequency  response  of  the  object  at  this  point  Hau(1)  represents 
the  combined  frequency  response  of  the  propagation  medium 
and  the  microphone  located  at  distance  l  from  the  object.  It  is 
assumed  that  Hau>(1)  is  unchanged  for  any  target  point  in  the 
object.  Now,  by  scanning  the  object  at  a  fixed-frequency  A uj  and 
recording  the  acoustic  emission  signal,  one  can  obtain  the  spatial 
distribution  of  QaujFalo  (within  a  constant  multiplier)  which 
can  be  mapped  into  an  image,  displaying  object  morphology. 

The  point-spread-function  (PSF)  of  this  system  is  defined  as 
the  image  of  a  point  target  located  on  the  focal  plane  [17].  This 
function  is  proportional  to  the  radiation  stress  field  on  the  focal 
plane.  This  function  can  be  found  in  an  analytical  form  by  com¬ 
bining  (3),  (4),  and  (6).  The  amplitude  of  this  function  can  be 
written  as  [17] 


exp(jA^). 


(9) 


Assuming  a  3 -MHz- transducer  with  the  outer  diameter  of  45 
mm,  inner  disk  diameter  of  29.6  mm,  and  a  gap  of  2  mm  between 
the  inner  disk  and  the  outer  ring  and  a  focal  length  of  70  mm,  the 
amplitude  of  the  resulting  PSF  function  for  A /  =  Alj/2'k  — 
13  kHz  is  shown  in  Fig.  2.  The  spatial  resolution  of  the  system, 
defined  by  the  diameter  of  the  central  lobe,  is  about  0.7  mm. 
The  depth  resolution  (or  the  focal  depth)  of  this  transducer  was 
experimentally  determined.  For  this  purpose,  a  small  glass  bead 
was  scanned  at  different  depths.  The  depth  resolution,  defined 
as  the  -6-dB  width  of  the  response  in  the  axial  direction,  was 
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found  to  be  about  9  mm.  The  peak  response  occurred  when  the 
bead  was  located  at  the  focal  distance  from  the  transducer. 

C.  Experimental  Setup 

The  experimental  setup  is  shown  in  Fig.  1.  The  experiments 
were  conducted  in  a  water  tank.  (In  a  system  designed  for  in  vivo 
imaging  the  transducers  can  be  placed  in  contact  with  the  skin 
instead  of  using  water.  The  hydrophone  would  also  be  placed 
in  contact  with  the  skin.  Because  at  low  frequencies  the  sound 
wave  propagates  almost  uniformly  in  all  directions,  hydrophone 
position  is  not  critical,  as  long  as  it  is  relatively  close  to  the 
exposure  site  but  not  in  the  ultrasound  path.)  A  two-element 
confocal  ultrasound  transducer  array  was  positioned  such  that 
the  beams  meet  the  object  at  their  joint  focal  point.  Transducer 
specifications  are  as  those  in  the  previous  section.  The  elements 
were  driven  by  two  stable  radio-frequency  (RF)  synthesizers 
(HP  33  120A)  at  frequencies  of  3  MHz  and  3  MHz  +  A/  Sound 
produced  by  the  object  vibration  was  detected  by  a  submerged 
hydrophone  (ITC  model  680)  placed  within  the  water  tank.  The 
received  signal  was  filtered  and  amplified  by  a  programmable 
filter  (Stanford  Research  Systems,  SR650)  to  reject  the  noise, 
then  digitized  by  a  12-bits/sample  digitizer  (National  Instru¬ 
ments  VXI-1000)  at  a  rate  sufficiently  higher  than  the  Nyquist 
rate.  Data  are  recorded  on  a  computer  disc. 

D.  Experiment  Procedure 

To  evaluate  the  capability  of  the  system  for  imaging  small 
particles  about  the  size  of  a  common  breast  microcalcification, 
we  constructed  a  test  object  comprised  of  four  small  glass  beads, 
ranging  from  260  ^m  to  400  fim  in  diameter  imbedded  in  a 
block  of  tissue-  mimicking  gel.  The  resulting  vibro-acoustog- 
raphy  image  is  shown  in  Fig.  3,  demonstrating  the  capability 
of  the  system  in  detecting  and  imaging  beads  that  are  at  least 
260  jum  in  diameter. 

Tissue  experiments  were  conducted  on  excised  human  breast 
tissue  samples.  These  tissues  were  obtained  postsurgically  and 


Fig.  3.  Image  of  glass  beads  in  gel.  This  image  is  obtained  by  scanning  the 
gel  phantom.  This  image  covers  an  area  of  20  mm  by  20  mm,  scanned  at  0.2 
mm/pixel.  The  beads  are  located  about  10  mm  deep  inside  the  gel.  The  diameter 
of  the  top  bead  is  about  260  // m  and  the  diameter  of  the  largest  bead  is  about 
400  // m.  The  vibration  frequency  was  set  at  53.5  kHz. 

were  fixed  in  formaldehyde  for  at  least  three  days  and  not 
more  than  six  months  at  the  time  of  experiments.  It  is  known 
that  formaldehyde  tends  to  harden  soft  tissue  and,  thus,  may 
introduce  an  additional  variable  in  the  experiments.  However, 
because  our  main  target  (microcalcifications)  is  much  harder 
than  the  fixed  tissue,  we  did  not  consider  tissue-hardening 
phenomenon  as  a  major  problem  in  our  experiments.  Patients’ 
records  were  reviewed  to  identify  tissues  with  high  probability 
of  having  microcalcifications.  Selected  tissues  were  cut  into 
(approximately)  3  x  3  cm  and  3-  to  5-mm-thick  pieces  and 
imaged  using  a  high-resolution  X-ray  mammography  machine. 
Therefore  microcalcifications  were  at  most  5  mm  from  the 
front  surface  of  the  sample.  The  X-ray  mammograms  were  then 
read  to  identify  presence  of  microcalcifications.  Tissue  pieces 
identified  with  microcalcification  were  each  mounted  flat  with 
a  few  droplets  of  glue  on  a  latex  sheet  and  secured  in  a  scanning 
bracket.  This  bracket  is  designed  to  hold  the  tissue  piece  in 
the  water  for  acoustic  scanning.  Glue  drops  were  carefully 
placed  on  areas  away  from  the  region  with  microcalcifications. 
Small  pieces  of  suture  were  sewed  to  the  tissue  and  used  as 
identification  marks.  For  this  purpose,  small  knots  are  placed  at 
different  positions  on  the  tissue  away  from  microcalcifications. 
These  markers  could  be  seen  in  vibro-acoustography,  X-ray 
(often  very  dim)  and  photographic  images. 

Specimen  X-ray  mammographic  images  were  obtained  from 
each  mounted  sample.  Each  X-ray  image  was  used  as  a  ref¬ 
erence  for  comparison  with  the  corresponding  vibro-acoustog¬ 
raphy  image.  The  system  illustrated  in  Fig.  1  was  used  to  ob¬ 
tain  vibro-acoustography  images  of  tissue  samples.  The  scan¬ 
ning  process  was  performed  by  a  raster  motion  of  the  specimen, 
the  ultrasound  beam  being  perpendicular  to  the  large  flat  surface 
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Fig.  4.  Vibro-acoustography  of  the  first  specimen,  (a)  Photograph  of  a  breast  specimen  mounted  on  the  scanning  bracket,  (b)  The  X-ray  mammogram  of  the 
selected  tissue  pieces  mounted  on  the  scanning  bracket.  Microcalcifications  can  be  seen  on  the  top  left  of  this  X-ray,  (c)  Vibro-acoustography  image  of  the  tissue 
at  25  kHz.  Microcalcifications  can  be  seen  in  this  image  as  bright  spots  on  the  top  left  and  (d)  Histology  of  the  tissue  specimen  around  the  region  with  some 
microcalcifications. 


Fig.  5.  Vibro-acoustography  of  the  second  specimen,  (a)  Photograph  of  the  tissue  specimen  mounted  on  the  scanning  bracket,  (b)  X-ray  mammogram  of  (a) 
showing  some  microcalcifications  at  the  center  of  the  image,  (c)  Vibro-acoustography  of  (a).  Microcalcifications  can  be  seen  as  bright  spots  at  same  location  as  in 
the  X-ray. 


of  the  tissue  layer  and  passing  through  the  3-  to  5-mm-thick- 
ness  of  the  tissue.  The  specimen  was  placed  on  the  focal  plane, 
such  that  its  entire  thickness  was  within  the  focal  depth  of  the 
transducer.  Experiments  were  conducted  in  a  water  tank  at  room 
temperature  (about  21  °C).  Sections  of  the  tissue  identified  with 
microcalcification  were  then  cut  for  histologic  study  to  identify 
microcalcifications  and  validate  the  imaging  results. 

III.  Results 

Fig.  4(a)  shows  the  photograph  of  a  breast  specimen  mounted 
on  the  scanning  bracket.  Sutures  are  knotted  on  the  tissue  for 
identification  purpose.  Fig.  4(b)  shows  the  X-ray  mammogram 
of  the  selected  tissue  pieces  mounted  on  the  scanning  bracket. 
The  sutures  are  barely  visible.  Microcalcifications  can  be  seen 
on  the  top  left  of  this  X-ray  in  a  high-density  sclerotic  region. 
This  sample  was  then  scanned  in  the  water  tank.  The  vibro- 
acoustography  image  of  the  tissue,  at  A /  =  25  kHz,  is  shown  in 
Fig.  4(c).  Microcalcifications  can  be  seen  in  this  image  as  bright 
spots.  Note  that  the  number  of  microcalcifications  matches  with 
the  corresponding  spots  in  the  X-ray  mammography.  Also  no¬ 
ticeable  is  that  dense  sclerotic  tissue  appears  dim  and  does  not 
interfere  with  identification  of  calcification  in  the  vibro- acoustic 
image.  Microcalcifications  identified  by  vibro-acoustography 
were  verified  histologically.  Fig.  4(d)  shows  the  histology  of  the 
region  around  the  calcifications.  The  breast  parenchyma  shows 
atypical  lobular  hyperplasia  and  microcalcifications  in  the  lob¬ 
ules.  We  measured  the  size  of  the  microcalcification  at  the  lower 


left,  by  measuring  the  spot  size  in  the  X-ray  mammogram  and 
found  to  be  approximately  110  (xm  in  diameter. 

Fig.  5(a)  shows  the  photograph  of  the  second  breast  tissue 
specimen  mounted  on  the  scanning  bracket.  Fig.  5(b)  is  the 
X-ray  mammogram  of  this  tissue  and  shows  the  microcal¬ 
cifications  at  the  center  of  the  image.  Fig.  5(c)  shows  the 
vibro-acoustic  image  of  this  tissue  at  22  kHz.  Microcalcifica¬ 
tions  can  be  seen  as  bright  spots  at  the  same  location  as  in  the 
X-ray.  The  breast  parenchyma  histologically  showed  comedo 
type  ductal  carcinoma  in  situ ,  with  calcifications  in  regions  of 
comedo  necrosis. 

Fig.  6(a)  shows  the  X-ray  mammogram  of  another  breast 
tissue  specimen  with  a  large  microcalcification.  Fig.  6(b)  shows 
the  corresponding  vibro-acoustic  image,  with  the  microcalcifi¬ 
cation  shown  as  a  bright  spot  within  the  relatively  dark  back¬ 
ground  of  soft  tissue. 

IV.  Discussion 

Experimental  vibro-acoustography  images  presented  in  this 
report  demonstrate  two  important  facts:  1)  The  vibro-acoustog¬ 
raphy  imaging  method  is  capable  of  detecting  small  microcalci¬ 
fications  in  breast  tissue;  2)  Microcalcifications  can  be  delin¬ 
eated  from  within  dense  and  sclerotic  tissues.  Results  shown 
here  also  indicate  that  such  images  have  high  spatial  resolution 
(700  fx m),  the  capability  of  detecting  small  (about  1 10  /im  in  di¬ 
ameter  and  above)  microcalcifications,  no  speckles,  good  con¬ 
trast  and  high  signal  to  noise  ratio. 
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Fig.  6.  Vibro-acoustography  of  the  third  specimen,  (a)  X-ray  mammogram  of  the  tissue  sample  showing  a  large  microcalcification  and  (b)  Vibro-acoustography 
image  of  (a). 


A  comparison  between  the  vibro-acoustography  image 
[Fig.  4(b)]  and  the  corresponding  X-ray  mammography 
[Fig.  4(c)],  reveals  an  interesting  capability  of  vibro-acous¬ 
tography  method.  High-density  sclerotic  tissue  in  the  region 
surrounding  the  microcalcification  absorbs  a  great  deal  of 
X-ray  energy;  hence,  the  X-ray  image  shows  very  low  contrast 
within  this  region  and  as  a  result  microcalcifications  are  barely 
visible.  However,  the  vibro-acoustography  image  delineates 
these  calcifications  with  high  contrast.  This  demonstrates  that 
some  tissues  that  are  radiologically  opaque  may  be  transparent 
in  vibro-acoustography  images  giving  the  potential  to  detect 
calcification  in  radiologically  dense  breasts. 

The  halo  seen  around  the  images  of  each  microcalcification 
is  an  artifact  of  the  vibro-acoustography  system.  This  effect  is 
seen  more  clearly  in  the  PSF  image  shown  in  Fig.  2.  Multiple 
rings  around  the  center  spot  are  due  to  the  sidelobes  of  the  Bessel 
functions  that  appear  in  the  PSF  function  (9)  as  discussed  in  the 
previous  section. 

Another  artifact  is  due  to  the  glue  used  to  bond  the  tissue 
specimens  to  the  latex  sheet.  We  used  small  droplets  of  instant 
glue  (Krazy  Glue  brand,  Elmers  Products,  Inc.  Columbus,  OH) 
for  this  purpose.  This  type  of  glue  produces  a  relatively  strong 
acoustic  emission  making  it  visible  in  vibro-acoustography 
images.  Fig.  7  shows  a  vibro-acoustography  image  of  three 
droplets  of  this  glue  on  a  latex  sheet,  resembling  images  of 
breast  microcalcification.  To  avoid  mistaking  glue  droplets 
for  microcalcifications,  we  marked  the  glued  region  of  tissue 
and  took  extra  precautions  to  locate  these  regions  on  the 
vibro-acoustography  images. 

The  choice  of  vibration  frequency  A  f  used  in  each  exper¬ 
iment  was,  to  some  extent,  determined  by  the  resonance  fre¬ 
quencies  of  the  water  tank.  The  acoustic  emission  field  resulting 
from  vibrations  of  microcalcifications  is  very  small.  To  facilitate 
the  detection  process,  we  often  chose  the  vibration  frequency 
to  correspond  to  one  of  the  resonance  frequencies  of  the  water 
tank.  This  way,  the  water  tank  acts  as  an  amplifier  to  improve 
the  signal-to-noise  ratio  (SNR).  Because  the  human  body  does 
not  have  a  regular  geometric  shape  with  parallel  walls,  such  res¬ 
onance  phenomenon  is  not  likely  to  take  place.  Therefore,  the 
choice  of  A /  may  no  longer  be  based  on  the  resonance  effect. 

Ultrasound  standing  waves  between  the  transducer  and  the 
object  may  cause  artifacts  in  the  image.  This  may  happen  if 
the  object  reflects  a  significant  portion  of  the  incident  beam 
back  to  the  transducer.  The  transducer  surface  can  then  act  as 
a  mirror  to  re-reflect  the  reflected  field  toward  the  object.  This 
process  may  occur  more  than  once.  The  resulting  field  at  the 
object  is  the  sum  of  the  incident  and  all  the  reflected  fields.  De- 


Fig.  7.  Vibro-acoustography  of  the  instant  glue  drops.  The  glue  drops  (about 
1  mm  in  diameter)  are  placed  on  a  latex  sheet. 

pending  on  the  distance  between  the  transducer  and  the  object 
and  the  ultrasound-  wavelength,  those  components  may  interact 
constructively  or  destructively.  Hence,  the  field  intensity  and 
the  resulting  image  can,  therefore,  become  excessively  sensi¬ 
tive  to  the  transducer-to-object  distance.  The  effect  of  standing 
wave  phenomena  on  the  image  diminishes  as  the  ultrasound  bath 
loss  increases.  Therefore,  it  is  expected  that  standing  wave  ar¬ 
tifacts  would  be  less  of  a  problem  in  clinical  applications  of 
vibro-acoustography  for  in  vivo  breast  imaging. 

In  the  present  study,  we  imaged  tissue  samples  only  in  two 
dimensions  (i.e.,  a  single  slice).  However,  the  ultrasound  beam 
used  here  produces  a  stress  field  that  is  confined  in  three  di¬ 
mensions.  Therefore,  it  is  possible,  in  principle,  to  selectively 
scan  several  slices  of  the  object  at  different  depths  to  produce  a 
three-dimensional  (volume)  image  of  the  object. 

Application  of  the  vibro-acoustography  for  breast  imaging  is 
both  promising  and  challenging.  An  expected  difficulty  in  using 
vibro-acoustography  on  patients  is  achieving  a  high  SNR.  Sev¬ 
eral  factors  that  can  reduce  the  SNR  below  what  we  have  in  the 
experimental  system.  These  factors  include  limitation  on  the  ul¬ 
trasound  intensity,  body  attenuation  of  the  ultrasound  beam,  de- 
focusing  of  the  ultrasound  beam  due  to  phase  aberrations,  and 
the  biological  noise  generated  by  the  body.  Imaging  speed  is 
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another  concern  in  in  vivo  vibro-acoustography.  The  scanning 
mechanism  used  in  the  experiment  presented  in  this  paper  scan 
the  object  one  point  at  a  time,  making  the  data  collection  a  rel¬ 
atively  lengthy  process  (up  to  a  few  minutes  per  image).  A  long 
imaging  time  is  not  desirable  for  breast  imaging  because  body 
motions  within  this  period  can  introduce  “motion  artifact  “  in 
the  images. 

There  are  several  benefits  in  using  vibro-acoustography 
for  in  vivo  breast  imaging  as  an  alternative,  or  in  parallel,  to 
X-ray  mammography.  One  advantage  of  vibro-acoustography 
stems  from  the  fact  that  it  uses  ultrasound  energy,  which  can 
easily  penetrate  in  dense  tissues.  Such  tissues  are  normally 
radio-opaque  and  difficult  to  image  with  X-ray.  Additionally, 
the  present  method  is  relatively  safe  and  hence,  may  be  useful 
in  applications  where  the  use  of  X-ray  is  limited  to  its  ionizing 
effects.  Vibro-acoustic  detection  and  follow-up  of  microcalcifi¬ 
cations  has  several  clinical  potentials,  including  the  following. 

•  Detection  of  microcalcifications  in  the  breast  of  women 
who  are  younger  and/or  have  radiologically  very  dense 
breast.  We  showed  that  ultrasound  can  penetrate  radiolog¬ 
ically  dense  tissue  and  detect  calcification. 

•  Detection  of  microcalcification  in  pregnant  or  lactating 
women. 

•  Detection  of  microcalcifications  in  women  with 
radio-opaque  breast  implants. 


V.  Summary 

In  this  paper,  we  present  a  noninvasive  imaging  method  for 
detecting  microcalcifications  in  breast  tissue.  This  method  uses 
ultrasound  in  a  fundamentally  new  way  to  image  the  tissue  at 
low  (kHz  range)  frequencies,  producing  high  signal-to-noise 
and  speckle-free  images. 

The  spatial  resolution  of  vibro-acoustic  images  is  determined 
by  the  focal  beamwidth  of  the  two-element  transducer,  which 
was  in  the  order  of  0.7  mm  for  the  transducer  used  here.  Ex¬ 
periments  were  conducted  on  human  breast  tissue  specimens. 
X-ray  tissue  mammography  images  were  used  for  initial  iden¬ 
tification  of  microcalcifications  and  also  as  reference  images 
for  comparison  with  the  corresponding  vibro-acoustography 
images.  Tissue  histology  was  also  used  to  validate  presence  of 
microcalcification. 

Results  indicate  that  microcalcifications  associated  with  a  va¬ 
riety  of  breast  pathology,  including  those  located  in  high-density 
or  sclerotic  regions  of  tissue,  can  be  detected  by  vibro-acoustog¬ 
raphy.  Microcalcifications  as  small  as  1 10  fim  in  diameter  were 
detected  by  vibro-acoustic  tissue  mammography. 

Further  development  of  vibro-acoustic  tissue  mammography 
may  lead  to  a  novel  imaging  tool  for  applications.  This  method 
can  be  considered  as  a  nonionizing  alternative  to  conventional 
X-ray  mammography  and  may  be  useful  for  such  cases  as 
pregnant  patients,  where  the  use  of  ionizing  radiation  is  not 
allowed.  Also,  vibro-acoustography  is  potentially  useful  for 
imaging  dense  breast  where  the  use  of  conventional  X-ray 
mammography  may  be  limited. 
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ABSTRACT  We  describe  theoretical  principles  of  an  imag¬ 
ing  modality  that  uses  the  acoustic  response  of  an  object  to  a 
highly  localized  dynamic  radiation  force  of  an  ultrasound  field. 
In  this  method,  named  ultrasound-stimulated  vibro-acoustogra¬ 
phy  (USVA),  ultrasound  is  used  to  exert  a  low-frequency  (in  kHz 
range)  force  on  the  object.  In  response,  a  portion  of  the  object 
vibrates  sinusoidally  in  a  pattern  determined  by  its  viscoelastic 
properties.  The  acoustic  emission  field  resulting  from  object 
vibration  is  detected  and  used  to  form  an  image  that  represents 
both  the  ultrasonic  and  low-frequency  (kHz  range)  mechanical 
characteristics  of  the  object.  We  report  the  relation  between  the 
emitted  acoustic  field  and  the  incident  ultrasonic  pressure  field 
in  terms  of  object  parameters.  Also,  we  present  the  point- spread 
function  of  the  imaging  system.  The  experimental  images  in  this 
report  have  a  resolution  of  about  700  /mm,  high  contrast,  and  high 
signal-to-noise  ratio.  USVA  is  sensitive  enough  to  detect  object 
motions  on  the  order  of  nanometers.  Possible  applications  in¬ 
clude  medical  imaging  and  material  evaluation. 


The  study  of  objects  in  terms  of  their  mechanical  response  to 
external  forces  is  of  considerable  interest  in  material  science 
and  medical  diagnosis.  Elastic  constants  are  closely  connected 
to  the  thermodynamic  properties  of  materials  and  can  be 
related  to  a  wide  range  of  physical  parameters.  Elastic  con¬ 
stants  can  be  determined  by  measuring  deformation  in  re¬ 
sponse  to  an  applied  force.  Although  a  static  force  can  be  used 
for  this  purpose,  using  a  dynamic  force  is  preferred  if  one  is 
interested  in  measuring  the  dynamic  characteristics  of  the 
material  (1). 

Changes  in  elasticity  of  soft  tissues  are  often  related  to 
pathology.  Palpation  is  a  traditional  example  of  estimating 
mechanical  parameters  for  tissue  characterization,  where  a 
static  force  is  applied  and  a  crude  estimation  of  the  tissue 
elasticity  is  obtained  through  the  sense  of  touch.  In  palpation, 
force  is  exerted  on  the  body  surface,  and  the  result  is  an 
accumulative  response  of  all  the  tissues  below.  Physicians  can 
sense  abnormalities  if  the  response  to  palpation  of  the  suspi¬ 
cious  tissue  is  sufficiently  different  from  that  of  normal  tissue. 
However,  if  the  abnormality  lies  deep  in  the  body,  or  if  it  is  too 
small  to  be  resolved  by  touch,  then  the  palpation  method  fails. 

Elasticity  imaging,  a  subject  extensively  investigated  in  re¬ 
cent  years,  is  a  quantitative  method  that  measures  the  me¬ 
chanical  properties  of  tissue.  The  general  approach  is  to 
measure  tissue  motion  caused  by  an  external  (or,  in  some 
methods,  internal)  force/ displacement  and  use  it  to  recon¬ 
struct  the  elastic  parameters  of  the  tissue.  Some  investigators 
have  used  static  force  to  compress  the  tissue  and  measured  the 
resulting  strain  by  ultrasound  (2, 3).  Others  have  used  external 
mechanical  vibrators  to  vibrate  the  tissue  and  detected  the 
resulting  displacement  in  tissue  by  Doppler  ultrasound  (4-7). 
For  a  review  of  elasticity  imaging  methods,  refer  to  ref.  8.  A 
recently  developed  method  uses  an  actuator  to  vibrate  the  body 
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payment.  This  article  must  therefore  be  hereby  marked  “ advertisement ”  in 
accordance  with  18  U.S.C.  §1734  solely  to  indicate  this  fact. 

PNAS  is  available  online  at  www.pnas.org. 


surface  and  then  measures  the  strain  waves  with  phase- 
sensitive  MRI  (9). 

Most  of  the  elasticity  imaging  methods  are  based  on  an 
external  source  of  force  resulting  in  a  spatially  wide  stress-field 
distribution.  This  requires  the  stress  field  to  pass  through  the 
superficial  portion  of  an  object  before  reaching  the  interior 
part.  Analysis  of  the  object  response  can  be  complicated 
because  the  stress-field  pattern  changes,  often  unpredict  ably, 
at  different  depths  before  it  reaches  the  region  of  interest 
within  the  object.  An  alternative  strategy  is  to  apply  a  localized 
stress  directly  in  the  region  of  interest.  One  way  to  accomplish 
this  is  to  use  the  radiation  pressure  of  an  ultrasound  source(s). 
Based  on  this  strategy,  Sugimoto  et  al.  (10)  presented  a  method 
to  measure  tissue  hardness  by  using  the  radiation  force  of  a 
single  focused  ultrasound  beam.  In  this  method,  impulsive 
radiation  force  was  used  to  generate  localized  deformation  of 
the  tissue.  Resulting  transient  deformation  was  measured  as  a 
function  of  time  by  an  ultrasound  Doppler  method.  Radiation 
force  has  also  been  used  to  generate  shear  elastic  waves  in 
tissues  (11). 

In  this  paper,  we  describe  the  principles  of  an  imaging 
technique  that  produces  a  map  of  the  mechanical  response  of 
an  object  to  a  force  applied  at  each  point.  The  method  uses 
ultrasound  radiation  force  to  remotely  exert  a  localized  oscil¬ 
lating  stress  field  at  a  desired  frequency  within  (or  on  the 
surface  of)  an  object.  In  response  to  this  force,  a  part  of  the 
object  vibrates.  The  size  of  this  part  and  the  motion  pattern 
depend  on  object  viscoelastic  characteristics.  The  acoustic 
field  resulting  from  object  vibration,  which  we  refer  to  as 
“acoustic  emissions, ”t  is  detected  by  a  sensitive  hydrophone 
and  used  to  form  the  image  of  the  object.  This  method  benefits 
from  the  high  spatial  definition  of  ultrasound  radiation  force 
and  high  motion-detection  sensitivity  offered  by  the  hydro¬ 
phone.  We  call  this  technique  ultrasound-stimulated  vibro- 
acoustography  (USVA).  Some  general  aspects  of  this  method, 
including  some  experimental  results,  have  been  outlined  by  the 
authors  in  ref.  12.  Here  we  present  the  theoretical  foundations 
of  USVA. 

METHODS 

Our  aim  is  to  image  an  object  based  on  its  mechanical 
characteristics.  This  is  achieved  by  vibrating  the  object  by 
applying  a  highly  localized  oscillating  force  to  each  point  of  the 
object.  The  localized  force  is  produced  by  modulating  the 
intensity,  and  thereby  the  radiation  force,  of  the  ultrasound  at 
low  frequencies  (normally  in  kHz  range).  The  resulting  sound 
emitted  by  the  object  is  a  function  of  object  mechanical 


Abbreviations:  USVA,  ultrasound-stimulated  vibro-acoustography; 
PSF,  point  spread  function;  CW,  continuous  wave. 
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Research,  Mayo  Clinic,  e-mail:  fatemi.mostafa@mayo.edu. 
tThe  term  “acoustic  emission”  is  used  to  describe  the  acoustic  field  in 
response  to  a  cyclic  vibration  of  the  object.  Similar  terminology  is  also 
used  in  the  field  of  nondestructive  testing  of  materials  and  in 
optoacoustic  imaging  to  describe  a  different  phenomenon,  usually  the 
acoustic  field  resulting  from  structural  deformation. 
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characteristics  and  the  location  of  the  excitation  point.  The 
image  is  produced  by  mapping  the  amplitude  or  phase  of  this 
sound,  which  is  detected  by  a  sensitive  hydrophone,  vs.  posi¬ 
tion.  Fig.  1  illustrates  this  method.  In  the  following  section,  we 
describe  the  relationship  between  the  USVA  image  and  the 
properties  of  the  object. 

THEORY 

Generation  of  a  Dynamic  Radiation  Force  on  a  Target.  The 

acoustic  radiation  force  is  the  time-average  force  exerted  by  an 
acoustic  field  on  an  object.  This  force  is  an  example  of  a 
universal  phenomenon  in  any  wave  motion  that  introduces 
some  type  of  unidirectional  force  on  absorbing  or  reflecting 
targets  in  the  wave  path.  Radiation  force  is  produced  by  a 
change  in  the  energy  density  of  an  incident  acoustic  field.  For 
a  review  of  this  phenomenon,  refer  to  ref.  1 3.  Consider  a  plane 
ultrasound  beam  interacting  with  a  planar  object  of  zero 
thickness  and  arbitrary  shape  and  boundary  impedance  that 
scatters  and  absorbs.  The  radiation  force  vector,  F,  arising 
from  this  interaction  has  a  component  in  the  beam  direction 
and  another  transverse  to  it.  The  magnitude  of  this  force  is 
proportional  to  the  average  energy  density  of  the  incident  wave 
(E)  at  the  object,  where  { )  represents  the  time  average  and  S , 
the  projected  area  of  the  object  (14) 

F  =  c !,$<£>,  [1] 

where  dr  is  the  vector  drag  coefficient  with  a  component  in  the 
incident  beam  direction  and  another  transverse  to  it.  The 
coefficient  dr  is  defined  per  unit  incident  energy  density  and 
unit  projected  area.  For  a  planar  object,  d,  is  numerically  equal 
to  the  force  on  the  object.  Physically,  dr  represents  the  scat¬ 
tering  and  absorbing  properties  of  the  object  and  is  given  by 
(14) 

dr  =  1  (n„  +  I1A  -  /  y  cos ctsdS)  +  q.S~!  J  Y  sin  ctsdS,  [2] 

where  p  and  q  are  the  unit  vectors  in  the  beam  direction  and 
normal  to  it,  respectively.  The  quantities  IIfl  and  Ilv  arc  the 
total  absorbed  and  scattered  powers,  respectively,  and  y  is  the 
scattered  intensity,  all  expressed  per  unit  incident  intensity. 
Also,  as  is  the  angle  between  the  incident  and  the  scattered 
intensity,  and  dS  is  the  area  clement.  The  drag  coefficient  can 
also  be  interpreted  as  the  ratio  of  the  radiation  force  magni¬ 
tude  on  a  given  object  to  the  corresponding  value  if  the  object 
were  replaced  by  a  totally  absorbing  object  of  similar  size.  This 
is  because  |dr|  =  1  for  a  totally  absorbing  object.  This  coeffi¬ 
cient  can  be  determined  for  objects  of  different  shapes  and 
sizes.  For  simplicity,  we  assume  a  planar  object  normal  to  the 
beam  axis.  In  this  case,  the  transverse  component  vanishes, 
thus  the  drag  coefficient  (force)  will  have  only  a  component 
normal  to  the  target  surface,  which  we  denote  by  scalar  dr  (F). 
Values  of  dr  for  spheres,  in  terms  of  the  diameter  and  the 
wavelength,  are  given  in  ref.  14. 

To  produce  a  dynamic  radiation  force,  one  can  use  an 
amplitude-modulated  beam  (15).  Consider  an  amplitude- 
modulated  incident  (ultrasonic)  pressure  field,  /?(/),  as 

pit)  -  PMtt  cos(Awf/2)  coscuf,/,  [3] 


Focused 

Ultrasound  Image 


Fig.  1.  Principle  of  ultrasound-stimulated  vibro-acoustography. 


where  PW{),  Aw/2,  and  wo  are  the  pressure  amplitude,  modu¬ 
lating  frequency,  and  center  frequency,  respectively.  In  our 
analysis  and  experiments,  we  assume  that  the  condition  Aw  <$C 
o)0  holds.  In  such  a  case,  the  energy  density  of  the  incident  field 
has  slow  variations  in  time.  To  discriminate  the  slow  time 
variations  of  a  function,  let  us  define  the  short-term  time  average 
of  an  arbitrary  function  £(0  over  the  interval  of  T  seconds  at 
time  instance  t,  as  <£(/)>-/•  =  1/7  I'Xrri  which  is  a 

function  of  t.  The  long-term  time  average  (or  simply  the  time 
average)  is  obtained  by  setting  T  -»  sc.  To  compute  the 
short-term  time  average  of  the  acoustic-energy  density  rele¬ 
vant  to  field  variations  at  Aw/2,  we  choose  T  longer  than  the 
ultrasound  wave  period  but  much  shorter  than  the  modulation 
period,  that  is  2tt/w{)  <$C  T  <$:  47t/Aw.  Under  this  condition, 
the  short-term  time  average  of  p2(t)  is  (p2(t))r  =  (P^/4)  (1  + 
cos  A  wf).  The  energy  density  is  given  by p2(t) /pc2,  where  p  and 
c  arc  the  density  and  propagation  speed  in  the  medium  (16). 
We  arc  interested  in  the  time-varying  component  of  the 
short-term  time  average  of  the  energy  density.  Denoting  this 
component  by  e±0,(t),  we  can  write:  e±(0(t)  =  (P^jApc2) 
cosAo )t.  This  component  of  the  energy  density  produces  a 
time-varying  radiation  force  on  the  target  (Eq.  1)  at  frequency 
Aw.  The  amplitude  of  this  force,  F± w,  is 

FlM  =  P\Sdr/4pc\  [4] 

This  equation  states  that  the  time-varying  force  amplitude  is 
proportional  to  the  square  of  incident  ultrasound  pressure,  or 
equivalently,  to  the  incident  power.  If  the  object  moves  in 
response  to  this  force,  then  the  high-frequency  ultrasound 
energy  would  convert  to  low-frequency  mechanical  energy. 

Acoustic  Emission  from  a  Target  Caused  by  a  Dynamic 
Force.  The  radiation  force  F±to  vibrates  the  target  object  at 
frequency  Aa>.  Object  vibration  results  in  an  acoustic  field  in 
the  medium  (acoustic  emission).  This  field  is  related  to  object 
shape,  size,  and  viscoelastic  properties.  To  present  a  conclusive 
analysis  of  this  relationship,  we  have  to  assume  an  object  with 
specific  characteristics.  Here  we  assume  that  the  vibrating 
object  has  a  circular  cross-section  of  radius  b  and  uniformly 
vibrates  back  and  forth  like  a  piston.  This  choice  allows  us  to 
illustrate  the  concept  in  a  simple  form.  We  also  consider  an 
area  S  ^  i rb2  of  the  piston  surface  to  be  projected  normally  by 
the  beam.  Similar  solutions  can  be  carried  out  for  other 
objects.^ 

The  steady-state  normal  velocity  amplitude  of  a  piston,  U 
caused  by  a  harmonic  force  F±0)  at  frequency  Aw,  can  be 
described  in  terms  of  the  mechanical  impedance  Z^w, 

=  [5] 

where  ZAw  =  Z*m  +  Zr  is  comprised  of  the  mechanical 
impedance  of  the  object  in  vacuum  Z’m ,  and  the  radiation 
impedance  of  the  object  Zr,  all  defined  at  Aw.  Modeling  the 
object  as  a  mass-spring  system,  Z,'„  can  be  written  in  terms  of 
Aw  as  (16,  17) 

Z;„  -  R’m  -  j(m So)  -  Kr /Ao>),  [6] 

where  m ,  R'm,  and  K'  are  the  mass,  mechanical  resistance,  and 
spring  constants  of  the  object,  respectively.  The  radiation 
impedance  of  the  piston  can  be  written  (17)  as  Zr  =  irb2(Rr  - 
jXr)  [7],  where  Rr  =  pc[  1  -  (c/Awfc)/i(c/2Aw&)],  [8],andYr 
=  (4pc/7r)  Jo/2  sin  [(26Aw/c)  cos  a]  sin2  ad  a  [9],  where  7i(*) 
is  the  first-order  Bessel  function  of  the  first  kind.  In  many 
applications  of  our  interest,  the  wavelength  is  much  greater 
than  the  object  size,  hence,  (Aw/c)6  0.  In  such  cases  Zr 

assumes  a  simpler  form  as  Zr  =  7rfc3pAw(£>Aw/2c  -  j$/3ir) 


$Thc  theory  can  be  extended  to  include  arbitrary  vibrating-part  shapes 
and  nonuniform  displacement  of  the  object.  Nonuniform  displace¬ 
ment  would  be  an  important  issue  when  the  vibration  wavelength  in 
the  object  material  is  smaller  than  2b. 
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[10].  The  mechanical  impedance  of  the  piston  object  can  now 
be  written  as 

ZAw  =  ( R'm  +  'rrb^'RO  -  j(m  Aw  -  K' /Aw  +  tt b2Xr) 

~  (R'm  +  rrpb4  Aw2 /2c)  -  j(m Aw  -  K' /Aw 

+  Spb3Aw/3),  fc(Aa>/c)  0.  [11] 


Once  we  calculate  UAw  we  can  calculate  the  pressure  field 
it  produces  in  the  medium.  We  assume  that  the  acoustic 
emission  signal  propagates  in  a  free  and  homogenous  medium. 
The  farfield  acoustic  pressure  caused  by  a  piston  source  of 
radius  b  set  in  a  planar  boundary  of  infinite  extent  is  given  by 
(17), 


Ptuo  =  -jAcop 


exp(/Aw//c) 
47 tI 


2J1[(bAw/c)  sin  tf]  ^  cos 


(bAw/c)  sin  ft  cos  +  fiB] 


X  (2t7^2EaJ, 


[12] 


where  /  is  the  distance  from  the  observation  point  to  the  center 
of  the  piston,  ft  is  the  angle  between  this  line  and  the  piston 
axis,  and  is  the  specific  acoustic  admittance  of  the  boundary 
surfaced  The  factor  of  two  comes  from  the  presence  of  the 
boundary  wall,  which  would  be  replaced  by  unity  if  the 
boundary  wall  were  not  present  (16).  The  acoustic  emission 
field  resulting  from  object  vibration  can  be  written  in  terms  of 
the  incident  ultrasound  pressure  by  combining  Eqs.  4,  5,  and 
12,  as 

{  Aw  exp{jAwl/c)  [U^bSw/c)  sin  ft]  cos  ft  11 
J  c2  X  477 i  (bAw/c)  sin  ft  X  cos  ft  +  fiB_  J 

X  {1  m;n  +  77 b2Rr)  —  j{m Aw  -  K! /Aw  +  Tib2XM2*b2)PlQSdr. 

[13] 

For  wavelengths  long  compared  to  the  object  size,  i.e.,  when 
b  A(o/c  — »  0,  the  term  in  the  first  brace  approaches  a  constant, 
hence  we  may  consider  the  contents  of  the  first  brace  to  be  an 
object-independent  function  (the  specific  acoustic  admittance 
ft b  relates  to  the  surrounding  boundary  surface).  Under  these 
conditions,  the  first  brace  in  the  above  equation  represents  the 
effect  of  the  medium  on  the  acoustic  emission  field,  which  we 
may  call  the  medium  transfer  function,  and  denote  it  by 

'IJ/fbAw/c)  sin  tf]  cos  ft 
(bAw/c)  sin  ft  cos  ft  +  fiB_ ' 

[14] 


Ao)  exp(/Aa J/c) 

HlJl}  ~J  c2  4itI 


The  second  brace  in  Eq.  13  is  1/Z Aa),  or  the  mechanical 
admittance  of  the  object  at  the  frequency  of  the  acoustic 
emission  (Aw),  and  we  denote  it  by  YA(0.  It  is  convenient  to 
combine  this  term  with  the  next  term  (277 b2)  in  Eq.  13,  as 
Gao>  =  27tZ?2Ya<0  =  27rZ>2/ZAtu,  which  is  the  total  acoustic 
outflow  by  the  object  per  unit  force  (acoustic  outflow  is  the 
volume  of  the  medium  (e.g.,  the  fluid)  in  front  of  the  object 
surface  that  is  displaced  per  unit  time  because  of  object 
vibration.).  Function  QA(0  represents  the  object  characteristics 
at  the  acoustic  frequency.  We  may  thus  rewrite  Eq.  13  in  a 
more  compact  form  as 

[15] 


Eq.  15  indicates  that  the  acoustic  emission  pressure  is  propor¬ 
tional  to:  (i)  the  square  of  ultrasound  pressure  P^;  (ii)  the 
ultrasound  characteristics  of  the  object,  dr,  in  the  projected 


area  S ;  (iii)  the  acoustic  outflow  by  this  object,  QA6>,  repre¬ 
senting  the  object  size  b  and  its  mechanical  admittance  at  the 
acoustic  frequency,  YAo);  and  (zv)  the  transfer  function  of  the 
medium  at  the  acoustic  frequency,  HAw(l)  .  The  above  equation 
illustrates  the  basic  nonlinear  relationship  between  the  ultra¬ 
sound  and  acoustic  emission  pressure  amplitudes.  Note  that 
neither  the  medium  nor  the  object  needs  to  be  nonlinear  for 
this  relationship  to  hold.  It  is  interesting  to  note  that  the 
projection  area  S  and  the  vibrating  area  irb2  play  different 
roles.  The  projection  area  determines  the  extent  of  the  force 
applied  to  the  object  (Eq.  4).  The  vibrating  area,  however, 
influences  the  total  acoustic  outflow  in  the  medium  caused  by 
object  vibration.  The  mechanism  of  object  vibration  is  some¬ 
what  analogous  to  that  of  a  loudspeaker,  where  the  electro¬ 
motive  force  is  exerted  at  a  small  area  of  the  membrane 
(usually  at  the  center),  causing  the  entire  membrane  surface  to 
vibrate.  In  our  method,  the  size  of  the  vibrating  area  depends 
on  the  object  structure.  For  a  free  suspended  point  object, 
smaller  than  the  beam  cross-section,  the  vibrating  area  would 
be  the  same  as  the  projected  area.  For  a  large  stiff  plate, 
however,  the  vibrating  area  could  be  much  larger  than  the 
projected  area  (similar  to  a  loudspeaker).  In  some  cases,  it  is 
more  convenient  to  write  the  acoustic  emission  field  in  terms 
of  the  applied  force  FAw.  Referring  to  Eq.  4,  we  can  rewrite  Eq. 
15  as 


Fao.  =  4pc2//Aw(/)gA(1/'Aa).  [16] 

Again  in  analogy  to  a  loudspeaker,  FA(t>,  QAb>,  and  HA(0(l) 
represent  the  electromotive  force,  dynamic  characteristics  of 
the  membrane,  and  propagation  medium  transfer  function. 

Beam  Forming.  To  probe  an  object  with  the  dynamic 
radiation  force  at  high  spatial  resolution,  it  is  ideal  to  confine 
the  dynamic  stress  field  to  a  very  small  region  in  three- 
dimensional  space.  We  may  define  the  resolution  cell  of  the 
system  as  the  volume  within  which  the  amplitude  of  the 
modulated  field  is  high  enough  to  produce  a  stress  field  on  a 
target.  The  purpose  of  beam  forming  is  to  produce  a  resolution 
cell  as  small  as  possible.  An  amplitude  modulated  single- 
focused  beam  can  provide  a  resolution  cell  that  is  small  in 
diameter  but  long  in  depth  direction.  A  superior  strategy  that 
can  achieve  a  small  resolution  cell  in  all  dimensions  is  to  use 
two  unmodulated  focused  beams  at  slightly  different  frequen¬ 
cies  and  allow  them  to  cross  each  other  at  their  focal  regions. 
This  is  accomplished  by  projecting  two  coaxial  confocal  con¬ 
tinuous-wave  (CW)  ultrasound  beams  on  the  object.  An  am¬ 
plitude-modulated  field  is  produced  only  at  the  interference 
region  of  the  two  unmodulated  beams  around  their  focal  areas, 
resulting  in  a  small  resolution  cell.  For  this  purpose,  elements 
of  a  two-element  spherically  focused  annular  array  (consisting 
of  a  central  disc  with  radius  a  i  and  an  outer  ring  with  the  inner 
radius  of  a  '2  and  outer  radius  of  a2)  are  excited  by  separate  CW 
signals  at  frequencies  =  w0  -  Aco/2  and  co2  =  cuo  +  Aco/2. 
We  assume  that  the  beams  are  propagating  in  a  lossless 
medium,  in  the  +z  direction  of  a  Cartesian  coordinate  system 
( x ,  y ,  z),  with  the  joint  focal  point  at  z  =  0.  The  resultant 
pressure  field  on  the  z  =  0  plane  may  be  written  as 

Pit)  =  Pi(r)  cos(w1t  +  fa (r))  +  Pz(r)  cos (w2t  +  ifr2W),  [17] 

where  r  =  Vx2-hy2  is  the  radial  distance.  The  amplitude 
functions  are  (16,  18) 

Pi(r)  =  pcUoi  (m^/AjZo)  jincb^/A^o),  [18] 


and 


P2(r)  =  pcEo2(7r/A2Zo)  [fl|  jinc(ra2/A2z0)  -  a'22  jinc(ra^/A2z0)],  [19] 


§The  specific  acoustic  admittance  is  j8  b  =  pc /Zb,  where  Zb,  the 
acoustic  impedance  of  the  boundary,  represents  the  ratio  between  the 
pressure  and  normal  fluid  velocity  at  a  point  on  the  surface. 


where  A/  =  lir/wt,  i  =  1,  2,  is  the  ultrasound  wavelength,  Uot 
is  the  particle  velocity  amplitude  at  the  i-th  transducer  element 
surface,  and  jinc  (X)  =  Ji(2irX)/ ttX.  The  phase  functions, 
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ipi(r)  ~  -  7T7'2/ A/Zo,  for  /  =  1,2,  arc  conveniently  set  to  be  zero 
at  the  origin. 

Now,  we  define  a  unit  point  target  at  position  (xo,  yo)  on  the 
focal  plane  with  a  drag  coefficient  distribution  as 

dr(x,  y)  =  &(x  -  x0.  y  -  y0),  [20] 

such  that  dr( x,  y)dxdy  is  unity  at  (x0,  y0)  and  zero  elsewhere. 
This  equation  is  merely  used  as  a  mathematical  model  because 
dr  is  physically  finite.  In  this  case,  the  projected  area  can  be 
considered  to  be  S  =  dxdy.  We  replace  drS  in  Eq.  1  with  dr(x , 
y)dxdy  and  follow  the  steps  similar  to  those  outlined  in  Eqs.  3 
and  4  for  the  pressure  field  expressed  by  Eq.  17,  then  the 
complex  amplitude  of  the  normal  component  of  the  force  on 
the  unit  point  target  can  be  found  as 

FAwCv(i,yo)  =  pUmUoz  ( T7a\/4\xz 0)  jinc(r0tf  i /A jinc(r„fl2/A2z()) 
-  (ira22/\2z{})  jincfr0fl2/A:z(,)]exp( -//^Aw/2 rz0),  [21] 


acoustic  outflow  by  the  object  when  unit  force  is  applied  at 
point  (x,y). 

Commonly,  an  imaging  system  is  studied  through  its  point- 
spread  function  (PSF),  which  is  defined  as  the  image  of  a  point 
object.  To  determine  the  PSF  of  our  system,  we  consider  a  unit 
point  target  at  the  origin  with  unit  mechanical  response, 
£?AW(*>y)  =  1*  Hence,  referring  to  Eqs.  20  and  22,  we  can  write 
g(x,y)  =  S(x,y).  To  obtain  the  PSF,  we  move  this  point  object 
to  every  possible  position  (x0,  y0)  on  the  z  =  0  plane  and  form 
the  image  using  the  resulting  acoustic  emission  field,  PAJx^ 
yn)-  Because  x0  andy0  are  now  being  treated  as  variables,  we 
may  replace  them  by  variables  x  andy,  respectively.  We  define 
the  normalized  PSF  of  the  coherent  imaging  system  as  the 
complex  function 

/*(A%y)  =  PAjA%y)/PAJ0,0).  [23] 

Division  by  Pau>(0,  0)  cancels  the  constant  multipliers.  Refer¬ 
ring  to  Eqs.  16  and  21,  we  can  write 


where  the  argumentsxo  andy0  are  added  to  denote  the  position 
of  the  point  target  and  r0  =  Vxl  +  y(2}.  Eq.  21  describes  the 
spatial  distribution  of  the  force  (the  stress  field).  This  equation 
shows  that  the  stress  field  is  confined  to  the  regions  near  the 
beam  axis  (r0  =  0)  and  decays  as  the  radial  distance  r{) 
increases.  The  lateral  extent  of  the  stress  field,  and  hence  the 
resolution  cell  diameter,  would  be  smaller  at  higher  ultrasound 
frequencies  (smaller  A]  and  A2).  One  can  calculate  the  total 
force  on  an  arbitrary  object  by  integrating  the  force  over  the 
projected  area.  The  axial  extent  of  the  resolution  cell  (depth 
resolution  or  the  depth  of  field)  can  be  determined  by  calcu¬ 
lating  the  force  FAu}  as  a  function  of  the  depth  variable,  in  a 
fashion  as  outlined  in  Eqs.  17  to  21.  For  conciseness,  we  will 
present  only  the  measured  values  for  the  depth  resolution  in 
Results. 

Loss  in  the  propagation  path  would  attenuate  both  ultra¬ 
sound  beams,  thus  less  radiation  force  would  be  generated  by 
the  remaining  ultrasound  energy.  In  the  case  of  soft  tissues,  the 
force  attenuation  factor  is  A(z{))  =  exp[oz0(wi  +  w2)],  where 
a  is  the  attenuation  coefficient  of  the  tissue.  Energy  loss  in  the 
medium  would  also  result  in  generation  of  a  separate  radiation 
stress  on  the  medium  along  the  ultrasound  paths.  However, 
because  the  two  beams  propagate  along  separate  paths  in  the 
CW  form,  they  exert  mainly  steady  radiation  stresses  to  the 
medium,  which  does  not  cause  object  or  medium  vibrations. 
Dynamic  radiation  force  is  produced  only  in  the  interference 
region  around  the  focal  area,  which  is  another  advantage  of 
using  two  unmodulated  beams  over  a  modulated  single  beam. 

Image  Formation.  To  produce  an  image,  we  scan  the  object 
in  a  plane  and  record  the  complex  amplitude  of  the  acoustic 
emission,  PAo„  at  different  positions.  In  this  process,  we  keep 
Aw  fixed.  For  transverse  view  images,  the  scan  plane  is  the  focal 
plane  (x-y).  Alternatively,  for  the  parallel  view  the  scan  plane 
is  the  x-z  plane.  In  the  conventional  ultrasound  imaging 
context,  these  two  views  are  called  the  C-scan  and  B-scan, 
respectively.  Our  main  focus  here  is  the  transverse  view 
imaging.  In  this  case,  the  acoustic  emission  data  obtained  by 
vibrating  the  object  at  point  (x,  y)  are  assigned  to  the  corre¬ 
sponding  point  (x,  y)  in  the  image. 

Before  defining  the  image,  we  need  to  define  the  function 
that  represents  the  object.  Referring  to  Eq.  15,  the  terms  that 
are  object  dependent  are  the  drag  coefficient  d,  and  the 
function  QAb)  (assuming  that  HAlo(l)  is  object  independent). 
The  object  function  g(x,  y)  is  thus  defined  as  the  spatial 
distribution  of  these  terms, 


g(x,  y)  =  Q±Jx,  y)dr(x ,  y).  [22] 

Variablesx  andy  are  added  to  denote  the  dependency  of  dr  and 
Qau  on  position.  In  particular,  0a(1)(a’,  y)  implies  the  total 


/i(x,y)  =  (a\  -  a‘22yl  jinc(ra,/AiZ0)[<7?  jincfrtf2/A2Zo) 

-  a 22  jinc(rtf  2/A2z0)]  exp(-y>2A<o/2cz()).  [24] 

This  equation  illustrates  that  the  system  PSF  is  a  circularly 
symmetric  function  with  the  peak  at  the  origin  and  decaying 
amplitude  with  increasing  the  radial  distance  r.  Amplitude 
decays  faster  for  higher  ultrasound  frequency.  This  function 
will  be  discussed  further  in  the  next  section. 

EXPERIMENTS 

The  experimental  setup  is  shown  in  Fig.  2.  The  confocal 
transducer  is  constructed  by  using  a  spherical  piezoelectric  cap. 
The  two  elements  are  constructed  by  dividing  the  back  elec¬ 
trode  of  the  piezoelectric  wafer  into  a  central  disc  and  the  outer 
ring,  such  that  the  elements  have  identical  beam  axes  and  focal 
lengths.  Radii  of  the  elements  are  a\  =  14.8  mm,  a2  =  22.5 
mm,  and  a2  =  16.8  mm,  and  the  focal  distance  is  70  mm. 
Transducer  elements  were  driven  by  two  stable  radio  fre¬ 
quency  synthesizers  [Hewlett-Packard  33120A  and  Analogic 
2045  (Peabody,  MA)]  at  frequencies  of  /0  -  A//2  and  /0  + 
A//2,  where /o  =  3  MHz,  and  the  value  of  A /  =  Aw/27t  is 
stated  separately  for  each  experiment.  The  object  was  placed 
at  the  focal  plane  of  the  ultrasound  beams  in  a  water  tank. 
Sound  produced  by  the  object  vibration  was  detected  by  an 
audio  hydrophone  (International  Transducer,  Santa  Barbara, 
CA,  model  680,  sensitivity  -154  dB  re  lV/pPa)  placed  within 
the  water  tank.  The  received  signal  was  filtered  and  amplified 
by  a  programmable  filter  (Stanford  Research  Sunnyvale,  CA, 
SR650)  to  reject  noise,  then  digitized  by  a  12  bits/sample 
digitizer  (Hewlett-Packard  E1429A)  at  a  rate  sufficiently 
higher  than  the  Nyquist  rate  for  the  particular  A/ used.  Data 
were  recorded  on  a  computer  disc.  For  coherent  imaging, 
which  requires  the  phase  information,  the  reference  signal 
(i.e.,  cos  Aw/)  was  obtained  by  electronic  downmixing  of  the 
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Fig.  2.  Ultrasound-stimulated  vibro-acoustography  system.  The 
confocal  ultrasound  annular  array  transducer  with  two  elements  is 
shown  on  the  left. 
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two  driving  signals  and  was  recorded  along  with  the  hydro¬ 
phone  signal.  The  relative  phase  of  the  acoustic  emission  data 
was  then  calculated  at  each  point  by  discrete  Hilbert  transform. 
We  conducted  two  experiments.  The  first  experiment  is  de¬ 
signed  to  verify  the  relationship  between  the  acoustic  emission 
pressure  and  the  ultrasound  pressure  (Eq.  15),  and  the  second 
is  designed  to  experimentally  measure  the  PSF  stated  in  Eq.  24. 
These  experiments  are  presented  to  prove  the  principles  of  the 
method.  Further  experiments,  illustrating  applications  of  the 
method  for  evaluation  of  object  mechanical  properties  and 
tissue  imaging,  are  presented  in  ref.  12. 

RESULTS 

Acoustic  Emission  vs.  Ultrasound  Pressure.  Eq.  15  states 
that  the  acoustic  emission  field  amplitude  is  linearly  propor¬ 
tional  to  the  square  of  the  incident  ultrasound  pressure,  or 
equivalently,  to  the  incident  power.  To  test  this  hypothesis,  a 
calibrated  1-mm-diameter  ultrasound  needle  hydrophone, 
with  its  tip  facing  the  ultrasound  beam,  was  used  to  measure 
the  ultrasound  field  at  the  focal  point.  The  tip  of  the  hydro¬ 
phone  also  served  as  an  object  to  generate  the  acoustic 
emission.  Here,  A/ was  set  at  40  kHz,  and  the  radial  distance 
from  the  tip  of  the  needle  hydrophone  to  the  audio  hydro¬ 
phone  was  about  65  mm.  The  result  is  shown  in  Fig.  3.  The 
slope  of  the  acoustic  emission  intensity  vs.  ultrasonic  intensity 
indicates  that  the  intensity  of  the  acoustic  emission  field  is 
proportional  to  the  square  of  the  ultrasound  intensity,  or 
equivalently,  the  acoustic  emission  pressure  amplitude  is  lin¬ 
early  proportional  to  the  ultrasound  power,  as  predicted  by  Eq. 
15.  In  another  experiment,  a  450-jum-diameter  glass  bead  was 
used  as  a  point  object.  In  this  case,  A/ was  set  at  either  7  or  40 
kHz.  The  radial  distance  from  the  glass  bead  to  the  audio 
hydrophone  was  about  50  mm.  Again,  the  data  indicate  a 
quadratic  relationship  between  the  acoustic  emission  and 
ultrasonic  intensities.  These  glass-bead  data  also  show  that 
increasing  the  frequency  increases  the  acoustic  intensity.  This 
can  be  better  understood  by  investigating  the  theoretical  model 
presented  in  Eq.  13.  If  we  assume  that  the  mechanical  admit¬ 
tance  of  the  object  is  approximately  constant  at  these  frequen¬ 
cies,  then  the  object  behaves  almost  as  a  point  source,  and  the 
acoustic  pressure  amplitude  is  proportional  to  Aw.  Hence  the 
acoustic  emission  intensities  IAb)1  and  IAb)2  at  frequencies  Ao>i 
and  Aco2,  respectively,  are  related  by  IAa>2/IA(0l  =  (Aa^/Awi)2. 
Now,  letting  Aaq  and  A w2  correspond  to  7  and  40  kHz, 
respectively,  we  can  write  the  intensity  ratio  as:  [I AftiJ I Awi]dB  = 
20  log  (Aw2/Awi)  =  15  dB.  The  mean  value  of  the  intensity 
ratio  calculated  from  the  glass-bead  data  at  7  and  40  kHz  is  16 
dB,  which  is  in  close  agreement  with  the  theoretical  result. 

PSF  Measurement.  To  demonstrate  the  image-formation 
process  and  support  the  theoretical  derivation  of  the  PSF  (Eq. 
24),  we  evaluated  this  function  experimentally.  For  this  pur¬ 
pose,  we  used  a  380-pm  diameter  glass  bead  as  a  model  for  a 
point  and  placed  it  on  a  thin  latex  sheet.  The  latex  sheet 
produces  only  a  small  change  in  the  incident  energy,  and  hence 
does  not  produce  significant  radiation  force  or  acoustic  emis¬ 
sion.  The  entire  object  was  placed  in  a  water  tank  and  the  latex 
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Fig.  3.  Acoustic-emission  field  intensity  vs.  the  combined  ultra¬ 
sound  intensity. 
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Fig.  4.  USVA  images  of  a  380-ju,m  glass  bead:  (4)  in-phase,  (. B ) 
quadrature,  (C)  phase,  and  ( D )  magnitude.  The  phase  in  C  ranges  from 
—  7r  radians  (black  regions)  to  +7T  radians  (white  regions),  and  was 
normalized  to  be  zero  at  the  center  of  the  glass  bead.  (Modified  with 
permission  from  ref.  12,  copyright  1998,  American  Association  for  the 
Advancement  of  Science.) 

sheet  surface  was  scanned  in  a  raster  format  at  0.2-mm 
increments  in  either  direction  at  A f  -  7.3  kHz.  The  amplitude 
and  phase  of  the  acoustic  emission  signal  were  calculated  at 
each  point  relative  to  the  reference  signal  data.  The  phase  was 
normalized  to  the  phase  value  at  the  center  of  the  bead.  The 
resulting  in-phase,  quadrature,  phase,  and  magnitude  images 
are  shown  in  Fig.  4A-D.  Transverse  image  resolution,  defined 
as  the  —  6-dB  width  of  the  bead  image,  is  approximately  700 
pm  for  the  in-phase  image  in  either  dimension  (refer  to  Fig.  5). 
To  compare  the  experimental  results  with  those  of  the  theory, 
we  calculated  the  profile  of  the  PSF  for  the  transducer 
parameters  used  in  this  experiment  according  to  Eq.  24.  Fig. 
5  shows  the  theoretical  PSF  profile  and  the  glass-bead  profile 
obtained  from  the  experiment  (in-phase  image  profile,  shifted 
to  center  at  zero).  This  figure  shows  excellent  agreement 
between  theory  and  experiment  for  amplitudes  above  20%  of 
the  peak.  The  experimental  data  show  some  background  offset 
about  12%  of  the  peak.  We  believe  that  this  background  offset 
is  caused  at  least  by  the  following  sources:  (i)  acoustic  emission 
by  the  latex  sheet;  ( ii )  the  background  acoustic  noise  in  the 
experimental  setup  caused  by  equipment  fans  and  some  struc¬ 
tural  and  building  vibrations;  (Hi)  nonlinearity  of  water  that 
can  produce  a  nonlinear  mix  of  the  two  beams  even  in  the 
absence  of  the  object;  and  (iv)  streaming  as  a  result  of  energy 
absorption  by  water  (14),  which  in  turn  can  vibrate  the  latex 
sheet  blocking  the  stream.  To  evaluate  the  depth  resolution  (or 
the  slice  thickness),  we  placed  the  glass  bead  on  the  beam  axis 
and  scanned  it  in  the  z-direction  about  the  focal  point.  The 
depth  resolution,  defined  as  the  distance  between  the  points 
where  the  amplitude  of  the  acoustic  emission  field  drops  to  -6 
dB  of  its  peak,  was  9  mm. 

DISCUSSION 

System  Properties.  Spectral  characteristics.  In  general, 
USVA  images  represent  object  characteristics  at  two  ends  of 
the  spectrum:  the  drag  coefficient  at  the  ultrasound  frequency 
and  the  mechanical  admittance  at  the  low  acoustic  frequency 
A/.  The  ultrasound  frequency  is  usually  set  at  a  value  suitable 
to  form  the  beam,  whereas  A/  can  vary  in  a  wide  range 
depending  on  the  application.  If  the  two  beams  are  produced 
by  similar  ultrasound  transducer  elements,  then  the  practical 


Fig.  5.  The  theoretical  PSF  profile  of  the  USVA  system  according 
to  Eq.  24  and  the  glass-bead  in-phase  image  profile  (Fig.  4.4)  obtained 
from  the  experiment. 
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upper  limit  for  Af  is  about  equal  to  the  transducer  bandwidth. 
The  lower  limit  of  Af  is  zero. 

Sensitivity.  System  sensitivity  in  detecting  very  small  displace¬ 
ments  is  an  important  practical  issue,  especially  when  the  allow¬ 
able  ultrasound  power  is  limited  (for  example,  in  medical  imag¬ 
ing).  Motion  measurement  with  ultrasound  pulse  echo  has  been 
used  previously  to  study  “stiffness”  of  tissues  (10).  However,  the 
sensitivity  of  ultrasound  pulse  echo  to  motion,  at  common 
medical  ultrasound  frequencies,  is  limited  to  several  micrometers. 
An  advantage  of  USVA  is  its  high  displacement  sensitivity.  Cyclic 
displacement  of  100  nm  at  10  kHz  produces  an  acoustic  intensity 
of  about  3.0  X  10“3  watts/cm2.  Hydrophones  similar  to  the  one 
used  in  these  experiments  are  sensitive  to  as  little  as  10“ 15 
watts/cm2  and  therefore  are  capable  of  detecting  very  small  cyclic 
displacements.  For  instance,  the  hydrophone  detected  an  acoustic 
pressure  of  about  15  X  10~3  Pa  at  a  distance  of  5  cm  from  the  glass 
bead  shown  in  Fig.  4.  Under  assumptions  of  isotropic  vibration, 
this  pressure  would  be  produced  by  a  similar-sized  sphere  vibrat¬ 
ing  with  a  displacement  amplitude  of  about  6  nm.  Sensitivity 
increases  at  higher  frequencies  because  the  acoustic  emission 
pressure  is  proportional  to  frequency  for  constant  mechanical 
admittance  (Eq.  13). 

Comparisons  with  pulse-echo  systems.  Some  contrasting 
features  of  USVA  with  respect  to  the  conventional  ultrasound 
pulse-echo  imaging  (B-mode  and  C-mode)  arc:  (/)  A  pulsc- 
echo  image  represents  object  microstructure  by  displaying  its 
ultrasonic  reflectivity  distribution.  The  acoustic  emission  sig¬ 
nal  in  a  USVA  system  is  proportional  to  the  drag  coefficient 
dr(x,  y),  which  is  a  local  ultrasound  parameter,  and  to  the 
function  Q\lo(x,y),  which  represents  the  bulk  response  of  the 
object  at  acoustic  frequency  Aio.  Hence,  a  USVA  image,  in 
general,  represents  both  the  microstructure  and  macrostruc¬ 
ture  of  the  object.  (//)  The  echo  signal  in  a  pulse-echo  system 
is  a  linear  function  of  the  incident  ultrasound  pressure  ampli¬ 
tude  and  the  amplitude  reflection  coefficient  of  the  object.  In 
USVA,  however,  the  acoustic  emission  signal  is  proportional 
to  the  ultrasound  power  and  the  power  reflection  coefficient  of 
the  object.  (Hi)  Pulse-echo  systems  are  not  directly  sensitive  to 
medium  absorption.  Absorption  is  indicated  as  relative 
changes  in  the  amplitude  of  the  echoes  resulting  from  the 
scatterers  within  the  medium.  In  a  USVA  system,  the  acoustic 
emission  can  be  produced  directly  as  a  result  of  energy 
absorption  by  the  medium,  even  if  the  medium  is  homogeneous 
(refer  to  Eq.  2).  (iv)  Pulse-echo  systems  are  generally  broad¬ 
band.  The  USVA  method  presented  here  is  basically  a  nar¬ 
rowband  technique,  (v)  Pulse-echo  systems  achieve  high 
depth  resolution  by  transmitting  short  wideband  pulses.  A 
USVA  system  gains  its  depth  resolution  by  tailoring  beam 
geometry  to  limit  the  depth  of  the  region  where  the  two  beams 
interfere.  A  USVA  system  does  not  require  a  wide  bandwidth 
signal  to  achieve  a  high  depth  resolution,  (vz)  The  data  of 
USVA  images  are  acquired  one  point  at  a  time,  which  resem¬ 
bles  the  data  acquisition  in  C-mode  pulse-echo  systems. 
B-mode  pulse-echo  systems,  however,  require  much  less  ac¬ 
quisition  time  because  the  data  are  collected  one  line  at  a  time. 

Applications.  USVA  promises  applications  in  two  general 
areas:  medical  imaging  and  material  evaluation. 

Medical  applications.  USVA  can  be  used  to  image  tissues  and 
evaluate  their  mechanical  characteristics.  To  use  USVA  for  in 
vivo  applications,  one  must  take  into  account  limitations  such  as 
safe  (ultrasound)  power  limit,  tissue  attenuation,  body  noise,  and 
phase  aberration.  The  ultrasound  power  required  to  generate  a 
detectable  acoustic  emission  depends  on  the  object,  acoustic 
noise,  and  receiver  sensitivity.  Experimental  results  shown  in  Fig. 
3  demonstrate  that  ultrasound  intensities  as  low  as  30  mW/cm2 
are  sufficient  to  detect  the  acoustic  emission  from  a  1-mm- 
diameter  object  with  our  hydrophone  in  the  water  tank.  This 
power  value  is  much  smaller  than  the  FDA  limit  for  safe 
diagnostic  ultrasound  applications.  Tissue  attenuation  reduces 
the  ultrasound  intensity  at  the  target,  and  hence  the  acoustic 


emission  [by  the  factor  A(z{))].  Attenuation  limits  the  usable 
ultrasound  frequency,  and  hence,  lowers  the  resolution.  It  also 
limits  the  signal-to-noisc-ratio  (SNR)  because  of  the  loss  in  the 
acoustic  emission  energy  as  a  result  of  ultrasound  attenuation  by 
tissue.  Direct  effect  of  tissue  attenuation  on  the  acoustic  emission 
signal  is  probably  negligible  because  attenuation  of  the  compres- 
sional  waves  at  frequencies  in  the  order  of  a  few  kHz  in  soft  tissues 
is  low.  Sources  of  biological  noise  of  the  human  body  include 
cardiovascular  and  respiratory  systems  and  muscle  movements. 
Body  noise  is  usually  concentrated  below  1  KHz  and  can  be 
filtered  out  if  Af  is  above  this  value.  The  SNR  can  be  improved 
by  increasing  the  time  duration  of  the  signal  recorded  at  each 
point  (to  increase  the  signal  energy)  and  by  using  very  narrow- 
band  filters  (to  reject  the  noise),  while  keeping  the  ultrasound 
power  within  the  safe  level.  Phase  aberration  in  tissue  can  reduce 
the  sensitivity  of  the  system  by  decreasing  the  effective  ultrasound 
energy  density  at  the  beam  interaction  region.  One  may  use 
known  phase  aberration  correcting  methods  to  reduce  such  an 
effect.  The  practical  value  of  these  methods  for  USVA  remains 
to  be  studied. 

Material  evaluation.  Another  field  in  which  USVA  can  be 
potentially  useful  is  material  characterization,  including  me¬ 
chanical  parameter  evaluation,  imaging,  and  nondestructive 
testing  of  materials.  USVA  can  be  used  for  detection  and 
imaging  flaws  in  materials.  Also,  one  may  use  USVA  to 
evaluate  the  mechanical  frequency  response  of  an  object  at  low 
frequencies.  In  such  case,  we  are  interested  in  determining 
Q±o>(x,  y)  vs.  frequency.  We  assume  that  the  object  is  uniform 
within  the  projected  area  S.  Then,  the  total  radiation  force  on 
this  object,  F&a„  can  be  calculated  by  integrating  F^fay)  over 

5.  Referring  to  Eq.  21,  one  can  show  that  for  Aa>  <sC  wo,  F±w 
is  virtually  independent  of  Aco.  If  H±tt)(l)  is  known  and  nonzero, 
then  the  function  Qa<.>(*>jO  can  be  estimated  using  Eq.  16  as 
Q±<o(x,  y)  =  Pao>(x,  y)/FA<tfA«>(l).  In  practice,  PA  Jx9  y)  is 
obtained  by  sweeping  A  to  in  the  range  of  interest  and  recording 
the  resulting  acoustic  emission  (12). 
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Ultrasound-Stimulated  Vibro-Acoustic 
Spectrography 

Mostafa  Fatemi  and  James  F.  Greenieaf 

An  ultrasound  method  based  on  radiation  force  is  presented  for  imaging  the  acoustic 
response  of  a  materia!  to  mechanical  excitation.  Acoustic  energy  was  emitted  from  solids 
and  tissues  in  response  to  an  oscillatory  radiation  force  produced  by  interfering  focused 
beams  of  ultrasound.  Frequency  spectra  of  ultrasound-stimulated  acoustic  emission 
exhibited  object  resonances.  Raster-scanning  the  radiation  force  over  the  object  and 
recording  the  amplitude  and  phase  of  the  emitted  sound  resulted  in  data  from  which 
images  related  to  the  elastic  compositions  of  the  acoustically  emitting  objects  could  be 
computed.  Acoustic  emission  signals  distinguished  tuning-fork  resonances,  submilli¬ 
meter  glass  spheres,  and  calcification  in  excised  arteries  and  detected  object  motions 
on  the  order  of  nanometers. 


Th  e  mechanical  response  of  objects  to  ex¬ 
ternal  forces  is  of  considerable  interest  in 
medical  diagnosis,  nondestructive  inspec¬ 
tion  of  materials,  and  materials  science.  An 
applied  force  is  often  used  to  produce  dis¬ 
placement  from  which  elastic  constants, 
like  spring  constants,  can  be  determined.  In 
resonant  ultrasound  spectroscopy,  an  ultra¬ 
sound  source  and  a  detector  are  used  to 
measure  the  resonance  frequencies  of  a  sam¬ 
ple  with  known  size  and  mass.  The  reso¬ 
nances  are  related  to  mechanical  parame¬ 
ters,  including  the  elastic  constants  of  the 


high-frequency  and  low-frequency  varia¬ 
tions  in  energy  density  result  at  the  inter¬ 
section  of  the  two  beams  produced  by  the 
two  elements.  Ultrasound-stimulated  acous¬ 
tic  emission  results  from  the  energy  term 
that  produces  a  low-frequency  vibration. 
The  low-frequency  force  on  a  target  at  the 
focal  point  can  be  computed  by 


wave  (CW)  ultrasound  beams  of  different 
frequencies  at  a  selected  point  on  the  ob¬ 
ject.  Interference  in  the  intersection  region 
of  the  two  beams  produces  sinusoidal  mod¬ 
ulation  of  the  ultrasound  energy  density. 
Modulation  of  the  energy  density  creates  an 
oscillatory  force,  effectively  vibrating  the 
object  at  the  selected  region.  The  resulting 
vibration  of  the  object  produces  an  acoustic 
field  [acoustic  emission  (7)]  that  can  be 
measured  some  distance  away. 

Ultrasound  beams  can  he  constructed  in 
a  variety  of  ways  for  this  purpose  (8).  We 


F,(t)  =  dr  (EUKJt,x,y))Tdxdy 


=  Cdr  cos(Ao)t)  (1) 

where  C  is  a  constant,  Aoj  =  | co t  —  o)2| ,  S 
is  the  area  over  which  E(;)C;i|(t,x',\),  the  total 
energy  density  in  the  focal  plane,  has  sig¬ 
nificant  value,  and  ()T  represents  a  short¬ 
term  time  average  (9).  For  focused  beams, 
the  intersection  region  can  be  small  enough 
that  Fj(t)  can  be  considered  to  be  an  oscil¬ 
lating  point  force  applied  to  the  object  at 
the  focal  intersection  of  the  beams. 

To  produce  an  ultrasound-stimulated  vi- 
bro-acoustic  spectrogram,  we  vibrate  a  small 
region  of  the  object  with  an  oscillating 
radiation  force  of  varying  frequency.  The 
complex  amplitude  of  the  resulting  acoustic 
emission  field  is 

$(Aoj)  =  CdrH(Ato)Q(Aw)  (2) 


material  (/)*  Recently,  a  magnetic  reso¬ 
nance  elastography  method  for  quantita¬ 
tively  measuring  the  displacement  of  tissues 
in  response  to  externally  applied  cyclic  forc¬ 
es  was  reported  by  Muthupillai  et  at.  (2). 
The  method  resulted  in  high-resolution  im¬ 
ages  of  the  shear  modulus  of  normal  and 
pathologic  tissues.  Others  have  used  ultra¬ 
sound  to  measure  tissue  displacement  asso¬ 
ciated  with  externally  applied  compressive 
and  cyclic  forces  (3). 

We  describe  an  imaging  technique  that 


used  two  coaxial,  confocal  transducer  ele¬ 
ments  of  a  spherically  focused  annular  array 
(consisting  of  a  central  disc  and  an  outer 
annulus)  driven  by  two  CW  signals  at 
slightly  different  frequencies  and  <o2 
(Fig.  1).  The  energy  density  at  a  point  in 
this  ultrasound  field,  say  at  the  focus,  is 
proportional  to  the  square  of  the  sum  of  the 
ultrasound  fields  from  the  two  elements. 
Squaring  the  sum  of  two  sines  gives  rise  to 
sum  and  difference  frequency  terms.  Thus, 


where  Q(Ao:>)  is  a  complex  function  repre¬ 
senting  the  mechanical  frequency  response, 
or  admittance,  of  the  object  at  the  selected 
point,  and  H(Agj)  represents  the  combined 
frequency  response,  or  transfer  function,  of 
the  propagation  medium  and  receiver  and  is 
assumed  to  be  fixed  and  known  (iO).  Re¬ 
cording  <KAto)  allows  us  to  obtain  Q(Aoo) 
for  each  point  within  a  constant  multiplier 
(/I).  We  raster-scan  the  radiation  force 
over  the  object  to  produce  data,  which  can 


uses  acoustic  emission  to  map  the  mechanical 
response  of  an  object  to  local  cyclic  radiation 
forces  produced  by  interfering  ultrasound 
beams.  Radiation  force  is  generated  by 
changes  in  the  energy  density  of  an  acoustic 
field  (4)-  For  instance,  a  collimated  ultra¬ 
sound  beam  impinging  normally  on  the  sur¬ 
face  of  an  object  of  arbitrary  shape  and 
boundary  impedance  will  produce  a  radiation 
force.  The  radiation  force  arising  from  this 
interaction  has  a  component  F  =  drs(E)  (5) 
in  the  beam  direction.  This  component  is 
proportional  to  the  time-average  energy  den¬ 
sity  of  the  incident  wave  (E),  the  projected 
area  of  the  object  s,  and  dT  (6),  the  scattering 
and  absorbing  properties  of  the  object. 

We  probe  the  object  by  arranging  the 
intersection  of  two  focused  continuous- 
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Fig.  1.  Experimental 
system  for  ultrasound- 
stimulated  vibro-acous- 
tic  spectrography:  a  two- 
element  confocal  ultra¬ 
sound  annular  array 
transducer,  consisting  of 
a  center  disc  and  an  out¬ 
er  ring.  The  elements  are 
driven  by  two  CW  sourc¬ 
es,  at  frequencies  equal 
to  <j)1  and  o>2  =  oj1  +  Aw, 
where  these  frequencies 


are  very  close  to  the  central  frequency  of  the  elements,  and  Aw  is  much  smaller  than  (<1%)  the  center 
frequency  of  the  ultrasound  transducer.  The  beams  interact  only  in  a  small  region  around  the  joint  focal 
point,  where  the  amplitude  of  the  field  oscillates  at  the  difference  frequency  Aw.  The  region  of  interest  is 
placed  at  the  joint  focal  point  and  is  probed  point-by-point  by  raster  scanning.  The  sound  field  resulting 
from  object  vibrations  at  each  position  is  received  by  a  hydrophone  and  recorded.  The  recorded  signal 
at  one  or  more  difference  frequencies  is  used  to  form  an  image  of  the  object.  The  experiments  were 
conducted  in  a  water  tank.  The  transducer  center  frequency  was  3  MHz;  its  outer  diameter  was  45  mm; 
and  it  was  focused  at  70  mm.  The  difference  frequencies  used  in  each  experiment  are  mentioned  in  the 
corresponding  legends. 
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be  mapped  into  a  pictorial  format.  The 
spatial  resolution  of  the  resulting  image  is 
determined  by  the  region  in  which  signifi¬ 
cant  interference  between  the  ultrasound 
beams  occurs  and  is  of  the  order  of  a  few 
wavelengths  at  the  ultrasound  frequency. 

Experiments  were  conducted  in  a  water 
tank,  which  provided  good  ultrasonic  and 
acoustic  coupling  to  the  object  and  freedom 
of  movement  for  the  prototype  scanner 
mechanism  (Fig.  1).  The  two-element  con- 
focal  ultrasound  transducer  array  was  posi¬ 
tioned  such  that  the  beams  interfered  at  the 
selected  region  of  the  object.  Sound  pro¬ 
duced  by  vibrations  of  the  object  is  approx¬ 
imately  omnidirectional  because  of  the  small 
size  of  the  vibrating  portion  of  the  object 
compared  with  the  wavelength.  This  sound 
was  detected  by  a  submerged  hydrophone 
placed  near  the  object  within  the  water  tank. 

To  test  the  hypothesis  that  ultrasound- 
stimulated  acoustic  emission  is  sensitive  to 
object  mechanical  properties  and  to  show 
how  such  properties  can  be  quantitatively 
evaluated  by  this  method,  we  produced  an 
ultrasound-stimulated  vibro-acoustic  spec¬ 
trogram  of  a  tuning  fork  immersed  in  isopro¬ 
pyl  alcohol  at  two  different  temperatures. 
We  aimed  the  focal  point  of  the  confocal 
transducer  at  a  fixed  position  on  one  of  the 
tines.  The  shear  viscosity  of  alcohol  changes 
with  temperature,  causing  a  slight,  but  de¬ 
tectable,  shift  in  the  spectrogram  (Fig.  2). 
The  shear  viscosity  T]  of  a  liquid  is  deter¬ 
mined  by  measuring  the  resonant  frequency 
/R  and  the  bandwidth  8/R  of  a  tuning  fork 
immersed  in  this  liquid  (12) 

k/r/S/r  b/R0\2 

(3) 

where  /R0  and  8/R0  are  the  resonant  fre- 


Fig.  2.  Vibro-acoustic  spectrograms  of  a  tuning 
fork  immersed  in  isopropyl  alcohol  at  two  different 
temperatures.  A  point  on  a  tine  of  the  fork  was 
vibrated  with  the  use  of  the  system  shown  in  Fig. 
1 .  The  difference  frequency  was  swept  from  1 600 
to  2000  Hz.  A  change  of  temperature  from  1 5°  to 
30°C  decreases  the  shear  viscosity  of  the  alcohol, 
which,  in  turn,  changes  the  resonance  frequency 
and  bandwidth  of  the  tuning  fork. 


quency  and  bandwidth  measured  in  vacu¬ 
um,  respectively,  and  p  is  the  liquid  density. 
The  constant  k  is  determined  experimen¬ 
tally.  The  measured  values  for  /R  and  8/R 
were  1769  and  5.76  Hz  at  15°C  and  1765.7 
and  4.53  Hz  at  30°C.  The  viscosity  of  iso¬ 
propyl  alcohol  is  reported  to  be  2.89  cP  (1 
centipoise  =  1  mPa*s)  at  15°C  (13),  From 
this  value,  the  constant  k  was  calculated. 
The  shear  viscosity  at  30°C  was  found  using 
Eq.  3  to  be  T)  =  1.77  cP,  which  is  the  same 
as  the  published  data  (13). 

We  tested  the  ability  of  ultrasonically 
stimulated  vibro-acoustic  spectrography  to 
image  the  frequency  response  of  different 
objects  with  identical  dr  by  scanning  three 
tuning  forks  with  different  resonant  frequen¬ 
cies.  A  color  acoustic  spectrogram  was  ob¬ 
tained  by  sweeping  the  frequency  of  the  ra¬ 
diation  force,  Ago,  in  a  range  covering  the 
resonant  frequencies  of  all  forks  at  each 
beam  position.  The  acoustic  emission  signal 
was  filtered  by  three  bandpass  filters  cen¬ 
tered  at  different  frequencies.  The  outputs  of 
these  filters  were  used  to  form  a  three-color 
composite  image  (Fig.  3).  The  forks  appear 


with  three  distinct  colors  because  each  fork 
responds  primarily  at  its  own  resonant  fre¬ 
quency.  Because  the  forks  were  made  from 
identical  materials,  other  ultrasound  imaging 
methods  would  not  be  capable  of  distin¬ 
guishing  these  objects. 

We  tested  the  hypothesis  that  the  ampli¬ 
tude  of  acoustic  emission  at  a  single  frequen¬ 
cy  can  be  used  to  detect  small,  highly  reflec¬ 
tive  isolated  objects.  We  scanned  a  380-jxm- 
diameter  glass  bead  placed  on  a  thin  latex 
sheet  and  recorded  the  amplitude  of  acoustic 
emission  (Fig.  4).  The  latex  sheet  produces 
only  a  small  change  in  the  incident  energy 
because  it  is  almost  transparent  to  the  ultra¬ 
sound  beam.  This  experiment  demonstrated 
the  ability  of  the  method  to  detect  isolated 
regions  of  hardness  with  respect  to  a  soft 
background. 

To  test  the  feasibility  of  using  the  tech¬ 
nique  to  image  mechanical  properties  of 
tissues,  we  measured  the  phase  and  ampli¬ 
tude  of  acoustic  emission  from  calcified  and 
noncalcified  excised  human  iliac  arteries. 
The  arteries  were  scanned  in  a  plane  per¬ 
pendicular  to  the  ultrasound  beam  axis. 


Frequency  scale  (Hz) 


Fig.  3.  Vibro-acoustic 
spectrogram  of  tuning 
forks.  (A)  Three  tuning 
forks  made  from  identical 
material  with  identical  tine 
cross  sections  (lengths 

are  different).  Resonance  500  1  ooo  1 500 

frequencies  in  water  are 
407  Hz  (right),  809  Hz 
(middle),  and  1709  Hz 
(left).  The  forks  were 
scanned  in  a  water  tank 
with  the  use  of  the  sys¬ 
tem  shown  in  Fig.1.  The 
scanning  plane  covers 

the  front  tines  at  the  bottom  part  of  the  forks.  At  each  position,  the  difference  frequency  was  swept  from 
250  to  2250  Hz.  The  ultrasound-stimulated  acoustic  emission  was  detected  with  the  hydrophone  and 
filtered  by  three  overlapping  bandpass  filters  with  frequencies  centered  at  500,  1000,  and  1500  kHz, 
respectively.  (B)  Color  acoustic  spectrogram  of  the  forks.  The  outputs  of  the  bandpass  filters  were  used  to 
produce  the  red,  green,  and  blue  image  components.  This  image  displays  two  characteristics  of  the 
object:  shape  and  frequency  response.  The  color  associated  with  each  fork  indicates  its  resonance 
frequency,  which  can  be  deduced  from  the  frequency  scale. 


Fig.  4.  Image  of  ultrasonically  stimulated  acoustic 
emission  amplitude  from  a  380- |xm -diameter 
glass  bead  placed  on  a  thin  latex  sheet.  The  latex 
surface  was  scanned  at  0.2-mm  increments  in 
each  direction.  The  difference  frequency  was  fixed 
at  7.3  kHz.  The  amplitude  of  the  acoustic  emission 
of  the  bead  (in  relative  units)  is  shown  in  gray  < 

scale.  The  latex  sheet  is  almost  transparent  to  the  < 

imaging  system.  The  glass  bead,  however,  pre¬ 
sents  a  large  acoustic  impedance  discontinuity, 
resulting  in  significant  oscillatory  radiation  force. 

The  magnitude  of  the  radiation  force  gives  rise  to 
large-amplitude  acoustic  emission,  thus  yielding  a 
high-contrast  image  of  the  bead.  The  image  re¬ 
veals  a  system  resolution  of  about  700  jxm.  The 
confocal  beam  geometry  leads  to  a  negative-amplitude  ring  around  the  bright,  positive  central  spot.  This 
effect  produces  edge  enhancement,  as  seen  on  artery  walls  in  Fig.  5. 
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Calcifications  within  the  arteries  produced 
distinctive  amplitude  and  phase  values 
when  compared  to  the  normal  arterial  walls 
(Fig.  5).  The  phase  of  the  oscillation  of 
driven  mechanical  systems  relative  to  the 
driving  force  depends  on  the  ratio  of  mass 
to  stiffness  (14).  Calcified  regions  of  the 
diseased  artery,  identified  by  an  x-ray  of  the 
sample,  produced  phase  shifts  in  acoustic 
emission  completely  different  from  that  of 
the  non  calcified,  and  thus  softer  or  less 
dense,  regions.  The  amplitude  images  are 
highly  detailed  and  exhibit  variations  in 
acoustic  emission  from  both  calcified  and 
uncalcified  regions  of  the  diseased  artery. 
These  differences  are  caused  by  variations 
in  the  product  of  the  reflection  properties  dr 
and  the  effective  mechanical  vibration  ad¬ 
mittance  properties  Q(Aoj)  of  the  tissue. 
Thus,  vibro-acoustic  spectrography  is  simi¬ 
lar  to  conventional  pulse -echo  ultrasound 
imaging,  which  is  sensitive  to  the  ultrasonic 
parameters  of  the  object,  but  has  the  advan¬ 
tage  of  also  being  sensitive  to  the  mechan¬ 
ical  admittance  Q(Ato)  at  low  frequencies. 

Motion  induced  by  ultrasound  and  mea¬ 
sured  with  ultrasound  pulse  echo  has  been 
used  previously  to  study  “hardness”  (15). 
However,  the  sensitivity  of  ultrasound  pulse 


Phase  (rad) 


Fig.  5.  Vibro-acoustic  spectrography  of  excised 
human  iliac  arteries.  (A)  X-ray  image  of  normal  (left) 
and  calcified  (right)  excised  human  iliac  arteries 
obtained  from  a  35-year-old  woman  and  a  67- 
year-old  man,  respectively.  Bright  areas  indicate 
calcifications.  (B)  Vibro-acoustic  spectrogram 
amplitude  image  at  a  fixed  difference  frequency  of 
6  kHz.  Calcification  details  appear  bright,  whereas 
the  arterial  walls  are  dim.  (C)  Phase  image.  Calci¬ 
fied  regions  produce  acoustic  emission  of  differ¬ 
ent  phase  with  respect  to  regions  of  the  tissue 
having  little  calcification,  as  indicated  by  the  x-ray. 


echo  to  motion  at  common  ultrasound  fre¬ 
quencies  is  limited  to  several  micrometers. 
The  advantage  of  ultrasound-stimulated  vi¬ 
bro-acoustic  emission  is  its  high  displace¬ 
ment  sensitivity.  Cyclic  displacement  of 
100  nm  at  10  kHz  produces  an  acoustic 
intensity  of  about  3.0  X  10“^  W/cm2.  Hy¬ 
drophones  similar  to  the  one  used  in  these 
experiments  arc  sensitive  to  as  little  as 
10“  H  W/cm2  and,  therefore,  can  detect 
very  small  cyclic  displacements.  For  exam¬ 
ple,  the  hydrophone  detected  an  acoustic 
pressure  of  about  1 5  X  10“  *  Pa  at  a  distance 
of  5  cm  from  the  glass  bead  shown  in  Fig.  4. 
Under  the  assumption  of  isotropic  vibra¬ 
tion,  this  pressure  would  be  produced  by  a 
similarly  sized  sphere  vibrating  with  a  dis¬ 
placement  amplitude  of  about  6  nm.  The 
method  will  be  more  sensitive  for  higher 
frequency  sound  because  acoustic  power  is 
proportional  to  the  square  of  frequency  for 
constant  displacement  amplitude.  The 
practical  upper  limits  for  the  difference  fre¬ 
quency  produced  with  modern  ultrasound 
transducers  is  about  equal  to  their  hand- 
width.  For  modern  transducers,  this  value  is 
80%  or  more  of  the  central  frequency  of  the 
transducer.  For  experiments  like  those  we 
conducted,  emission  frequency  well  in  ex¬ 
cess  of  1  MHz  could  be  produced.  The  lower 
limit  on  the  frequency  of  radiation  pressure 
is  zero,  that  is,  static  pressure. 

Ultrasound-stimulated  vibro-acoustic 
spectrography  has  potential  applications  in 
at  least  two  general  areas.  The  first  is  non¬ 
destructive  evaluation  of  materials,  where 
material  characteristics  and  structural  flaws 
can  be  identified  by  measuring  changes  in 
the  mechanical  response  to  vibration  at  a 
point.  The  object  under  test  could  be  re¬ 
motely  vibrated,  for  instance,  by  beams 
propagating  and  interfering  in  either  water 
or  air,  or  beams  propagating  within  the 
object  could  be  used  to  produce  acoustic 
emission  from  flaws.  For  medical  imaging 
and  detection,  the  technique  appears  par¬ 
ticularly  suitable  for  noninvasive  detection 
of  hard  tissue  inclusions,  such  as  the  imag¬ 
ing  of  arteries  with  calcification,  detection 
of  breast  microcalcifications,  visualization 
of  hard  tumors,  and  detection  of  foreign 
objects.  The  stiffness  of  soft  tissues  is  re¬ 
lated  to  their  composition  (for  example, 
relative  values  of  fibrotic  content),  and  its 
change  is  often  related  to  pathology  or 
therapy.  In  conventional  palpation,  phy¬ 
sicians  estimate  tissue  stiffness  by  feeling 
with  the  fingers.  Because  changes  of  stiff¬ 
ness  alter  the  vibration  frequency  response 
or  damping  of  tissue,  the  present  method 
can  potentially  provide  a  noninvasive,  re¬ 
mote,  high-resolution  “palpation”  tech¬ 
nique  that  can  reach  small  abnormalities 
that  arc  otherwise  untouchable  by  con¬ 
ventional  methods. 
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where  Ha  and  U  are  the  total  absorbed  and  scattered 
powers,  respectively:  y  is  the  scattered  intensity  (”a, 
T7l>t  and  y  are  expressed  per  unit  incident  intensity):  p 
is  the  angle  between  the  incident  and  scattered  inten¬ 
sities:  s  is  the  projected  area  of  the  object:  and  dS  is 
the  area  element.  The  real  and  imaginary  parts  repre¬ 
sent  the  components  of  the  force  parallel  and  perpen¬ 
dicular  to,  respectively,  the  incident  field  momentum. 
In  our  treatment,  we  assume  a  planar  object  normal  to 
the  beam  axis:  hence,  df  is  real,  and  the  force  has  only 
a  normal  component  to  the  object  surface. 
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resulting  from  structural  deformation,  cracking,  or 
thermal  expansion  of  the  object.  We  note  a  funda¬ 
mental  difference  betv/een  the  method  we  present 
and  that  of  opto-acoustic  imaging.  In  ultrasound- 
stimulated  vibro-acoustic  spectrography,  the  ultra¬ 
sound  energy  is  converted  directly  to  low-frequency 
acoustic  energy  by  the  object,  whereas  the  opto- 
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signal  results  in  a  field  that  is  not  spatially  confined, 
producing  a  radiation  force  on  any  object  (including 
the  transducer  itself)  that  happens  to  be  in  the  beam 
path.  The  use  of  two  single -frequency  beams  is  ad¬ 
vantageous  because  field  modulation  occurs  only  in 
a  confined  region,  the  size  of  which  is  controlled  by 
the  intersection  of  the  two  beams. 

9.  We  define  the  short-term  time  average  of  an  arbitrary 
function  of  time  f{t)  around  time  instance  t  as 


T.  2 


The  long-term  time  average  is  obtained  when  T  — >  x. 
To  compute  the  short-term  time  average  of  the 
acoustic  energy  density  relevant  to  acoustic  emis¬ 
sion  at  Aw,  we  choose  T  longer  than  the  ultrasound 
wave  period  but  much  shorter  than  the  acoustic 
wave  period,  that  is,  2 tt/m2  <$r7  <sc2tt AW 

1 0.  In  Eq.  2,  H(Aw)  can  be  position -invariant  if  the  geom¬ 
etry  of  the  propagation  medium  remains  unchanged 
during  the  scan.  In  our  experiments,  we  minimized 
position  dependency  by  fixing  the  position  of  the 
transducer  relative  to  the  hydrophone  and  moving, 
instead,  the  object  in  raster-scanning  motion. 

1 1 .  Changing  Aw  by  shifting  the  frequency  of  the  ultrasound 
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Introduction 

A  common  approach  for  evaluating  some  of  the  characteristics 
of  an  object  is  listening  to  the  sound  that  object  makes.  This 
approach  is  used  for  many  purposes  in  our  daily  life.  Evaluating 
the  quality  of  a  crystal  glass  by  tapping  a  finger  on  it  is  a  simple 
example.  Plucking  a  violin  string  to  adjust  the  tuning  is  another. 

It  is  well  known  that  changes  in  elasticity  of  soft  tissues  are 
often  related  to  pathology.  Physicians  traditionally  use  sound  as 
a  diagnostic  tool,  for  example,  by  tapping  on  the  patient's  body. 
Vibro-acoustography  implements  the  same  concept  but  in 
micro-scale,  and  in  a  way  that  can  be  applied  to  living  tissues. 
This  method  uses  a  small  point-like  tapping  force  that  vibrates  a 
region  within  the  body,  while  a  sensitive  microphone  is  used  to 
listen  to  the  “sound”  of  the  vibrating  region. 

In  analogy  to  the  sound  of  the  crystal  glass,  sound  from  the 
body  carries  some  information  about  the  structure  and  the 
material  properties  of  the  region  in  the  body.  For  example,  the 
loudness  of  the  sound  may  indicate  the  hardness  of  the  body 
region  being  vibrated.  Also,  regions  with  different  elasticity  may 
respond  differently  at  different  frequencies.  Sound  from  body 
vibrations  can  be  used  to  produce  an  image  of  the  body.  To  do 
this,  the  point-like  tapping  force  is  moved  across  the  region  of 
interest  in  a  raster  motion,  and  the  sound  emitted  from  each 
point  is  used  to  determine  the  brightness  of  the  corresponding 
point  on  the  image  display  system.  The  result  is  a  display  of 
elastic  properties  of  the  region  of  interest. 

In  this  paper,  we  first  describe  how  vibro-acoustography  works, 
then  present  results  of  some  experiments  on  biological  tissues. 

Vibro-acoustography 

Vibro-acoustography  produces  a  map  of  the  mechanical 
response  of  an  object  to  a  dynamic  force  applied  at  each  point. 
The  method  utilizes  ultrasound  radiation  force  to  remotely  exert 
a  localized  force  (a  force  that  is  confined  to  a  small  region)  at 
a  desired  frequency  within  (or  on  the  surface  of)  an  object, 
and  records  the  resultant  acoustic  response12.  This  acoustic 
response,  which  is  normally  at  low  kHz  range,  is  a  function  of 
the  viscoelastic  properties  of  the  object  and  can  be  used  to 
produce  an  image  of  the  object3. 

Before  introducing  the  system,  we  shall  briefly  discuss  the 
phenomenon  of  the  radiation  force.  The  phenomenon  of 
acoustic  radiation  force  and  radiation  pressure  has  been  studied 
for  a  century4’5.  The  acoustic  radiation  force  is  an  example  of  a 
universal  phenomenon  in  any  wave  motion  that  introduces  some 
type  of  unidirectional  force  on  absorbing  or  reflecting  targets  in 
the  wave  path.  For  example,  the  radiation  force  of  sunlight 
pushes  the  tail  of  comets  away  from  the  sun.  Similarly,  a  sound 
wave  can  produce  the  radiation  force  on  objects  it  interacts 
with.  Vibro-acoustography  utilizes  an  oscillating  radiation  force 
of  ultrasound  to  vibrate  the  tissue,  and  consequently,  produce 
an  acoustic  emission  from  the  object.  To  produce  an  oscillating 
radiation  force  the  intensity  of  the  incident  ultrasound  must  be 
amplitude  modulated  at  the  desired  low  frequencies.  Using  a 
single  amplitude  modulated  beam  seems  to  be  the  simplest 
means  to  attain  this  purpose.  However,  such  a  beam  will  exert 
a  radiation  force  on  any  object  that  is  present  along  the  beam 


path,  producing  undesirable  acoustic  emission.  To  confine 
oscillations  of  the  radiation  stress  to  the  desired  region,  vibro- 
acoustography  uses  two  unmodulated  continuous  wave  beams 
at  slightly  different  frequencies,  propagating  along  separate 
paths.  The  beams  are  arranged  to  cross  each  other  at  their 
respective  foci,  and  thus  produce  a  modulated  field  at  a 
confined,  small,  cross  over  region.  This  beam  produces  a  force, 
oscillating  at  the  difference  frequency,  on  any  object  located  at 
the  focal  region. 

Figure  (1)  illustrates  a  vibro-acoustography  system.  The 
elements  of  a  two-element  annular  array  transducer  are  driven 
by  two  continuous  wave  signals  at  frequencies  f±  and  f2.  The 
object  to  be  imaged  is  placed  at  the  joint  focal  plane  of  the 
transducer  elements  (the  scanning  plane).  The  ultrasound  field 
produces  a  localized  radiation  stress  at  the  focal  point  on  the 
object  at  frequency  f2  -  fi.  Depending  on  the  elastic  properties  of 
the  object,  the  radiation  force  may  cause  a  portion  of  the  object, 
or  the  entire  object,  to  vibrate  at  this  frequency.  The  acoustic 
emission  resulting  from  object  vibration  is  received  by  a 
hydrophone  (or  microphone)  that  is  located  nearby.  Normally, 
the  wavelength  of  the  vibration  is  large  compared  to  the  object 
size,  hence  the  acoustic  emission  field  is  almost 
omnidirectional.  Therefore,  hydrophone  position  is  not  a  critical 
parameter  in  this  measurement.  The  filter  is  used  to  eliminate 
noise  and  other  interfering  signals.  To  form  an  image,  the  focal 
point  of  the  transducer  is  moved  across  the  scanning  plane 
within  the  object  in  a  raster  pattern.  The  acoustic  emission  is 
received  at  each  position,  and  an  image  is  formed  by  displaying 
the  amplitude  (or  phase)  of  such  signals  at  corresponding 
positions  on  the  image  plane.  The  spatial  resolution  of  this 
imaging  method  is  determined  by  the  ultrasound  beamwidth  at 
the  focal  plane,  which  is  normally  comparable  to  the  incident 
ultrasound  wavelength. 

Microphone 


Fig.  1.  Vibro-acoustography  system.  The  object  is  vibrated  by  the  radiation 
force  of  the  confocal,  two-beam  transducer.  Resuiting  vibrations  are  detected 
by  the  microphone  (or  hydrophone)  and  mapped  into  the  image. 


Experiments 

In  medical  applications,  one  can  use  vibro-acoustography  to 
obtain  images  of  the  human  body  for  diagnostic  purposes.  To 
demonstrate  the  capability  of  vibro-acoustography,  we  used  this 
method  to  image  a  specimen  of  human  iliac  artery  that  included 
calcifications.  A  3  MHz  transducer,  with  a  diameter  of  45mm 
and  focal  length  of  70mm,  was  used  for  this  experiment. 

The  difference  frequency  was  set  at  44  kHz  (this  is  also  the 
frequency  of  object  vibrations).  Figure  (2a)  shows  a  photo  of  the 
artery  that  is  cut  open.  The  areas  with  calcification  are  marked 
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by  letter  “C”.  The  vibro-acoustic  image  of  the  same  artery  is 
shown  in  Fig.  (2b).  Calcium  deposits  are  formed  in  large  plates, 
ranging  in  size  from  a  few  mm  to  a  few  cm.  Such  calcium 
plates,  which  are  much  stiffer  than  the  arterial  wall,  constitute 
efficient  acoustic  radiators  and  can  produce  a  strong  acoustic 
field  when  vibrated  by  the  radiation  force  of  the  incident 
ultrasound,  allowing  the  calcified  regions  to  stand  out  in  the 
vibro-acoustography  image. 


j?-  imaie  of  an  ex-vivo  human  iliac  artery  near  the  bifurcation 

on  the  left  (a)  is  a  photo  of  the  cut-open  artery  specimen.  Calcified  areas  are 

withjetter  "C".  Vibro-acoustic  image  of  the  specimen  is  shown  on  the  rig! 
[bJ-  Calcified  areas  stand  out  as  bright  spots .  This  image  was  obtained  with  a  3 
MHzconfbcal  transducer  as  described  in  the  text.  The  difference  frequency  was  si 
at44  The  gray- level  scale  is  in  arbitrary  units.  (Reproduced  with  permission 
from  Reference  6J 


The  next  example  illustrates  detection  of  micro-calcifications  in 
breast  tissue.  Figure  (3a)  shows  a  photograph  of  a  piece  of 
human  breast  tissue  collected  from  a  57-year-old  subject.  Two 
small  sutures  are  knotted  on  the  tissue  for  identification 
purpose.  Figure  (3b)  shows  the  x-ray  mammogram  of  this  tissue 
piece.  The  sutures  are  barely  visible  in  this  mammogram.  This 
sample  was  scanned  by  a  vibro-acoustography  system  similar  to 
the  one  used  in  the  previous  experiment.  The  resulting  vibro- 
acoustic  image  of  the  tissue  at  25  kHz  is  shown  in  Fig.  (3c). 
Micro-calcifications  can  be  seen  in  this  image  as  bright  spots.  It 
is  interesting  to  note  that  the  number  and  position  of  micro¬ 
calcifications  matches  with  the  corresponding  spots  in  the  x-ray 
mammogram.  Also  noticeable  is  that  tissue  inhomogeneities 
appear  dim  and  do  not  interfere  with  calcification  in  the  acoustic 
image.  Sutures  are  also  visible  in  the  vibro-acoustic  image, 
apparently  because  they  contain  entrapped  gas. 

The  above  experimental  results  show  that  vibro- 
acoustography  can  be  used  for  high-quality  imaging.  The  spatial 
resolution  is  comparable  to  the  ultrasound  beamwidth  at  the 
focal  region.  For  example,  for  the  3  MHz  ultrasound  transducer 
used  in  the  experiments,  the  spatial  resolution  is  about  700  pm. 
These  results  also  show  that  vibro-acoustic  images  have  high 
signal-to-noise  ratio  and  no  speckle.  These  features  allow  vibro- 
acoustography  to  delineate  calcium  deposits  with  high  definition 
and  contrast. 


In  summary,  elastic  properties  of  tissue  appear  to  be  closely 
related  to  tissue  pathology.  These  properties  also  relate  to  our 
perception  of  how  the  tissue  feels.  Therefore,  they  are  readily 
understandable  by  the  physician.  The  value  of  the  information 
on  elastic  properties  of  tissue  cannot  be  overemphasized. 
Methods  aimed  at  imaging  or  evaluating  elastic  features  of 
tissue  promise  a  new  era  in  diagnostic  imaging. 


y.  wcooi  uoouc  zauifjit:.  svu^ru-caicmcavons  are  located 
at  the  top  left  corner  of  the  sample  (indicated  by  the  arrow).  Two  sutures  are 
knotted  on  the  tissue  as  indicators,  (b)  Tissue  x-ray  mammogram  of  the  tissue 
sample.  Thts  sample  includes  four  small  micro-calcifications  at  the  top  left  comer. : 
The  brightness  of  the  image  has  been  reduced  to  make  the  micro-calcifications 
stand  out  from  the  high  density  tissue  at  the  top  left  comer  of  the  sample.  As  a 
result,  most  of  the  tissue  is  not  visible  in  the  x-ray  image,  except  for  the  calcified 
and  sorn?  brigfit  irregular  regions  (containing  no  calcification)  at  the  center 
imagt?ft  me  samp!e  in  (a}  and  {b)-  fora#  clearly  shows  four 
micro-calcifications  as  bright  spots  at  the  top  left  comer  of  the  tissue  sample.  The 
two  sutures  are  also  indicated  in  the  image.  The  location  of  micro-calcifications 
shown  here  matches  well  with  those  shown  in  the  x-ray  mammogram. 
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Abstract.  Elasticity  and  viscosity  of  soft  tissues  are  often  related  to  pathology. 
These  parameters,  along  with  other  mechanical  parameters,  determine  the  dynamic 
response  of  tissue  to  a  force.  Tissue  mechanical  response,  therefore,  may  be  used 
for  diagnosis.  Measuring  and  imaging  of  the  mechanical  properties  of  tissues  is  the 
aim  of  a  class  of  techniques  generally  called  elasticity  imaging  or  elastography.  The 
general  approach  is  to  measure  tissue  motion  caused  by  a  force  or  displacement 
and  use  it  to  reconstruct  the  elastic  parameters  of  the  tissue.  The  excitation  stress 
can  be  either  static  or  dynamic  (vibration).  Dynamic  excitation  is  of  particular  in¬ 
terest  because  it  provides  more  comprehensive  information  about  tissue  properties 
in  a  spectrum  of  frequencies.  In  one  approach  an  external  stress  field  must  pass 
through  the  superficial  portion  of  the  object  before  reaching  the  region  of  interest 
within  the  interior.  An  alternative  strategy  is  to  apply  a  localized  stress  directly 
in  the  region  of  interest.  One  way  to  accomplish  this  task  is  to  use  the  radiation 
force  of  ultrasound.  This  approach  offers  several  benefits,  including:  (a)  safety  - 
acoustic  energy  is  a  noninvasive  means  of  exerting  force;  (b)  adaptability  -  exist¬ 
ing  ultrasound  technology  and  devices  can  be  readily  modified  for  this  purpose; 
(c)  remoteness  —  radiation  force  can  be  generated  remotely  inside  tissue  without 
disturbing  its  superficial  layers;  (d)  localization  —  the  radiation  stress  field  can  be 
highly  localized,  thus  allowing  for  precise  positioning  of  the  excitation  point;  and  (e) 
a  wide  frequency  spectrum.  Several  methods  have  been  developed  for  tissue  probing 
using  the  dynamic  radiation  force  of  ultrasound,  including:  (a)  transient  methods 
which  are  based  on  impulsive  radiation  force;  (b)  shear-wave  methods  which  are 
based  on  generation  of  shear-waves;  and  (c)  vibro-acoustography,  recently  devel¬ 
oped  by  the  authors,  where  a  localized  oscillating  radiation  force  is  applied  to  the 
tissue  and  the  acoustic  response  of  the  tissue  is  detected  by  a  hydrophone.  Here, 
we  focus  on  vibro-acoustography  and  present  a  detailed  description  of  the  theory 
and  the  experimental  results.  We  conclude  with  the  capabilities  and  limitations  of 
these  radiation- force  methods. 


1  Introduction 

It  is  well  known  that  changes  in  elasticity  of  soft  tissues  are  often  related 
to  pathology.  Traditionally,  physicians  use  palpation  as  a  simple  method  for 
estimating  mechanical  properties  of  tissue.  In  palpation,  a  static  force  is  ap¬ 
plied  and  a  crude  estimation  of  tissue  elasticity  is  obtained  through  the  sense 
of  touch.  The  force  is  applied  on  the  body  surface  and  the  result  is  a  collec¬ 
tive  response  of  all  the  tissues  below.  Clinicians  can  sense  abnormalities  if 
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the  response  to  palpation  of  the  suspicious  tissue  is  sufficiently  different  from 
that  of  normal  tissue.  However,  if  the  abnormality  lies  deep  in  the  body,  or  if 
it  is  too  small  to  be  resolved  by  touch,  then  the  palpation  method  fails.  The 
dynamic  response  of  soft  tissue  to  a  force  is  also  valuable  in  medical  diagno¬ 
sis.  For  instance,  rebound  of  tissue  upon  sudden  release  of  pressure  exerted 
by  the  physician’s  finger  on  the  skin  provides  useful  diagnostic  information 
about  the  tissue. 

Quantitative  measurement  of  the  mechanical  properties  of  tissues  and 
their  display  in  a  raster  format  is  the  aim  of  a  class  of  techniques  gener¬ 
ally  called  elasticity  imaging  or  clastography  [1].  The  general  approach  is  to 
measure  tissue  motion  caused  by  an  external  (or,  in  some  methods,  internal) 
force  or  displacement  and  use  it  to  reconstruct  the  elastic  parameters  of  the 
tissue.  The  excitation  stress  can  be  either  static  or  dynamic  (vibration).  The 
dynamic  excitation  is  of  particular  interest  because  it  provides  more  compre¬ 
hensive  information  about  tissue  properties  in  a  spectrum  of  frequencies,  or 
alternatively,  the  transient  behavior  of  the  tissue  could  be  deduced  from  the 
measurements.  In  most  elasticity  imaging  methods  ultrasound  is  used  to  de¬ 
tect  the  motion  or  displacement  resulting  from  the  applied  stress.  Magnetic 
resonance  elastography  is  a  recently  developed  method  [2]  that  employs  a  me¬ 
chanical  actuator  to  vibrate  the  body  surface  and  then  measures  the  strain 
waves  with  phase-sensitive  magnetic  resonance  imaging  (MRI).  The  majority 
of  elasticity  imaging  methods  are  based  on  an  external  source  of  force,  result¬ 
ing  in  a  spatially  wide  stress  field  distribution.  This  requires  the  stress  field  to 
pass  through  the  superficial  portion  of  the  object  before  reaching  the  region 
of  interest  within  the  interior.  This  requirement  can  complicate  the  estima¬ 
tion  of  stiffness  because  the  stress-field  patterns  change  at  different  depths. 
Also,  because  the  stress  field  is  widely  distributed,  the  response  of  the  object 
is  an  integral  sum  of  the  stress  field.  An  alternative  strategy  is  to  apply  a 
localized  stress  directly  in  the  region  of  interest.  One  way  to  accomplish  this 
is  to  use  the  radiation  pressure  of  ultrasound. 

Acoustic  radiation  force  is  the  time-average  force  exerted  by  an  acoustic 
field  on  an  object.  This  force  is  produced  by  a  change  in  the  energy  density  of 
an  incident  acoustic  field  [3],  for  example,  due  to  absorption  or  reflection.  Sev- 
eial  benefits  may  result  from  using  ultrasound  radiation  force  for  evaluating 
tissue  pioperties,  including:  (a)  acoustic  (ultrasound)  energy  is  a  noninvasive 
means  of  exerting  force;  (b)  existing  ultrasound  technology  and  devices  can 
be  readily  modified  for  this  purpose,  thus  eliminating  the  need  for  develop¬ 
ing  a  new  technology;  (c)  radiation  force  can  be  generated  remotely  inside 
tissue  without  disturbing  its  superficial  layers;  (d)  the  radiation  stress  field 
can  be  highly  localized,  thus  allowing  for  precise  positioning  of  the  excitation 
point;  and  (e)  radiation  force  can  be  produced  in  a  wide  range  of  frequencies 
or  temporal  shapes.  These  features  make  radiation-force  methods  more  at¬ 
tractive  than  other,  mostly  mechanical,  excitation  methods  used  in  elasticity 
imaging. 
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Tissue  probing  with  the  radiation  force  produced  by  ultrasound  can  be 
accomplished  in  a  variety  of  techniques,  depending  on  the  excitation  and  de¬ 
tection  methods  used.  We  may  categorize  these  methods  as  follows:  (a)  tran¬ 
sient  methods,  where  an  impulsive  radiation  force  is  used  and  the  transient 
response  of  the  tissue  is  detected  by  Doppler  ultrasound  [4] ;  (b)  shear-wave 
methods,  where  an  oscillating  radiation  is  applied  to  the  tissue  and  the  re¬ 
sulting  shear  wave  is  detected  by  ultrasound  or  other  methods  [5,6,7];  and 
(c)  vibro-acoustography,  recently  developed  by  the  authors,  where  a  localized 
oscillating  radiation  force  is  applied  to  the  tissue  and  the  acoustic  response 
of  the  tissue  is  detected  by  a  hydrophone  [8,9].  Figure  1  illustrates  the  rela¬ 
tionship  between  various  methods  that  have  been  developed  for  evaluating 
or  imaging  the  elastic  properties  of  tissue. 

Dynamic  radiation-force  methods  seem  to  be  evolving  rapidly  as  a  new 
field  in  tissue  characterization  and  imaging.  The  purpose  of  this  chapter  is 
to  systematically  discuss  the  features  and  capabilities  that  are  offered  in  this 
group  of  elasticity  measurement  techniques.  Here,  we  describe  the  general 
approach  used  in  the  transient,  shear- wave,  and  vibro-acoustography  meth¬ 
ods.  We  pay  particular  attention  to  the  last  method,  present  its  theory,  and 
discuss  its  features.  Finally  we  discuss  the  capabilities  and  limitations  of  all 
three  methods. 


Fig.  1.  Tissue  elasticity  evaluation  and  imaging  methods.  These  methods  can  be 
divided  into  two  general  groups.  The  static  group  employs  a  steady  force  to  deform 
the  object.  The  dynamic  group  use  either  a  momentary  or  an  oscillating  force 
to  vibrate  the  object  and  can  be  divided  into  two  subgroups.  Mechanical-drive 
methods  use  a  mechanical  actuator,  whereas  the  radiation-force-drive  methods  use 
the  radiation  force  of  ultrasound  to  excite  the  object.  The  latter  subgroup  is  divided 
into  three  subsubgroups  based  on  the  detection  method  used.  Vibro-acoustography 
uses  a  hydrophone  to  detect  object  response 
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2  Theory  of  the  Radiation  Force 


The  study  of  radiation  force  and  radiation  pressure  dates  hack  nearly  one 
century,  to  the  time  of  Rayleigh.  [10].  Since  then,  this  subject  has  been  under 
continuous  investigation.  A  historical  review  of  radiation  force  and  radiation 
pressure  is  presented  in  [11],  and  a  critical  review  of  the  subject  can  be  found 
in  [3].  Some  recent  analyses  of  radiation  force/pressure  in  an  attenuating 
medium,  which  may  be  applicable  to  biological  tissues,  are  presented  in  ]12l 
and  [13]. 

The  acoustic  radiation  force  is  an  example  of  a  universal  phenomenon 
in  any  wave  motion  that  introduces  some  type  of  unidirectional  force  on 
absorbing  or  reflecting  targets  in  the  wave  path.  Radiation  force  in  fluids  is 
often  studied  in  the  context  of  radiation  pressure.  Depending  on  the  boundary 
conditions,  radiation  pressure  can  be  defined  differently.  Simple  explanations 
of  these  definitions  can  be  presented  by  considering  a  sound  traveling  inside, 
and  along  the  axis  of,  a  cylindrical  container  toward  the  opposite  wall  [11], 
Rayleigh  radiation  pressure  is  the  excess  pressure  produced  on  the  opposite 
wall  when  the  container’s  side  wall  is  confining  the  fluid  inside.  Langevin 
radiation  pressure  is  the  excess  pressure  on  the  opposite  wall  when  we  remove 
the  confining  side  wall,  so  that  the  fluid  is  free  to  move  (Fig.  2).  Here  we  will 
focus  on  Langevin  radiation  pressure  because  the  conditions  for  which  this 
pressure  is  defined  apply  to  our  experimental  situation.  It  can  be  shown  that 
the  Langevin  radiation  pressure  of  a  plane  wave  impinging  normally  on  a 
perfectly  absorbing  wall  is  equal  to  the  total  energy  density  {E),  where  (. . .} 
represents  the  time  average.  If  the  wall  is  partially  reflecting,  this  pressure 
would  be  equal  to  (1  +  R)(E),  where  R  is  the  power  reflection  coefficient  [11]. 
Thus,  in  general,  we  can  write  the  radiation  force  of  a  normally  impinging 
sound  beam  on  a  wall  as 


F  =  drS(E ),  (1) 

where  F  is  the  force  in  the  beam  direction,  and  d,r  is  the  “radiation- force 
function  or  the  drag  coefficient.  This  dimensionless  coefficient  is  defined  per 
unit  incident  energy  density  and  unit  projected  area.  For  a  planar  object,  d.r 
is  numerically  equal  to  the  force  on  the  object.  Physically,  the  drag  coefficient 
represents  the  scattering  and  absorbing  properties  of  the  object  [14],  For  a 
perfectly  absorbing  object  dr  =  1,  and  for  a  perfectly  reflecting  object  dr  =  2. 
In  the  case  of  oblique  incidence,  the  radiation  force  will  have  a  normal  as  well 
as  a  transverse  component.  A  more  detailed  description  of  dr  is  presented 
in  [14], 

To  produce  a  time-varying  radiation  force,  the  intensity  of  the  incident 
beam  can  be  modulated  in  various  ways.  For  example,  a  short  ultrasound 
pulse  can  produce  a  transient  pulsed  radiation  force,  and  a  sinusoidally  mod¬ 
ulated  beam  can  result  in  a  sinusoidally  varying  force. 
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Fig.  2.  Radiation  force  produced  by  projection  of  ultrasound  on  an  object.  Inter¬ 
action  of  an  ultrasound  beam  with  an  object  that  scatters  and/or  absorbs  results 
in  a  force  on  the  object  in  the  beam  direction.  The  magnitude  of  this  force  is  pro¬ 
portional  to  the  time- average  energy  density  of  the  incident  beam  and  a  factor  that 
represents  the  scattering  and  absorbing  properties  of  the  object 


3  Radiation-Force  Methods 

In  transient  methods  the  radiation  force  of  ultrasound  is  used  to  make  a 
minute  deformation  in  the  tissue.  The  transient  recoil  of  the  tissue  resulting 
from  this  deformation  is  measured  and  used  for  evaluation  of  tissue  elas¬ 
tic  properties.  A  method  for  measuring  tissue  hardness,  presented  by  Sugi- 
motoe tal.  [4],  uses  the  radiation  force  of  a  single  focused  ultrasound  beam. 
Ideally,  hardness  may  be  represented  by  the  spring  constant  of  the  object, 
which  is  the  ratio  of  the  applied  force  to  the  displacement.  Principles  of  the 
transient  methods  are  presented  in  Fig.  3. 


Fig.  3.  Principle  of  the  transient  method.  The  excitation  transducer  produces  a  mo¬ 
mentary  radiation  force  which  induces  a  minute  deformation  in  the  object.  Doppler 
methods  are  used  to  measure  this  deformation  using  the  probe  transducer 
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3.1  Transient  Method 

In  this  method  an  ultrasound  pulse  is  used  to  generate  a  short-duration  ra¬ 
diation  force  which  produces  localized  deformation  of  the  tissue.  Immedi¬ 
ately  after  the  force  pulse,  the  resulting  transient  deformation  of  the  tissue 
is  measured  as  a  function  of  time  with  Doppler  ultrasound  using  a  separate 
transducer.  The  deformation  includes  an  initial  rapid  squeeze  of  the  tissue, 
followed  by  a  relaxation  and  possibly  a  rebound.  The  deformation  is  a  func¬ 
tion  of  tissue  viscoelastic  parameters,  as  well  as  the  applied  force.  The  authors 
argue  that,  because  quantitatively  measuring  the  internal  radiation  force  is 
difficult,  it  is  advantageous  to  derive  a  relative  quantity  that  is  representative 
of  tissue  hardness.  To  derive  a  single  relative  quantity  from  the  deformation 
data,  the  relaxation  part  of  the  function  is  approximated  by  a  sum  of  several 
exponential  curves,  and  the  sum  of  the  first-order  derivatives  of  such  expo¬ 
nentials  is  calculated.  Sugimoto  etal.  [4]  show  that  this  quantity  is  correlated 
to  the  spring  constant  of  the  tissue;  thus  it  may  be  used  as  a  measure  of 
tissue  hardness. 

3.2  Shear- Wave  Methods 

The  shear  modulus  is  related  to  the  hardness  or  elasticity  of  the  material. 
It  is  known  that  the  shear  moduli  of  various  soft  tissues  range  over  several 
orders  of  magnitude,  while  the  bulk  modulus,  a  parameter  that  is  associated 
with  the  conventional  pulse  echo  ultrasound  compressional  wave  speed,  varies 
significantly  less  than  an  order  of  magnitude  [15,16].  These  features  indicate 
that  the  shear  modulus  may  be  a  better  parameter  for  tissue  characterization 
than  bulk  modulus.  Tissue  attenuation  of  the  shear  wave  is  very  large,  even 
at  low  kHz  frequencies.  One  way  to  induce  localized  shear  waves  inside  tissue 
is  to  use  the  radiation  force  of  focused  ultrasound  [5],  In  a  method  called 
shear- wave  elasticity  imaging  (SWEI)  [6],  an  amplitude-modulated,  single- 
focused  ultrasound  beam  is  used  to  induce  a  localized  radiation  stress  inside 
the  soft  tissue.  Localization  of  the  stress  field  is  critical  to  the  success  of  the 
method.  To  achieve  a  high  degree  of  localization,  the  method  uses  a  focused 
ultrasound  beam.  It  is  shown  that  the  radiation  stress  exerted  within  a  dissi¬ 
pative  medium  peaks  about  the  focal  region  of  the  highly  focused  transducer. 
Also,  it  has  been  suggested  that  localization  can  be  improved  by  designing  the 
transducer  and  selecting  the  beam  parameters  such  that  a  nonlinear  shock 
wave  is  produced  in  the  focal  region,  increasing  the  magnitude  of  the  stress 
field  in  the  vicinity  of  the  focal  region,  thus  augmenting  localization  of  the 
stress  field. 

Modulation  of  the  ultrasound  beam  can  be  in  the  form  of  an  oscillating 
wave  or  short  pulse.  The  resulting  radiation  force  elicits  a  shear  wave  prop- 
agating  in  the  radial  direction  with  respect  to  the  beam  axis,  with  particle 
motion  parallel  to  the  beam  axis.  Shear  waves  in  soft  tissue  travel  at  very 
low  speed,  typically  around  a  few  meters  per  second;  thus  the  corresponding 
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wavelength  is  much  shorter  than  that  of  the  compressional  waves  for  the  same 
frequency.  Shear  waves  are  also  highly  attenuated  in  soft  tissue,  with  an  at¬ 
tenuation  coefficient  two  or  three  orders  of  magnitude  higher  than  that  of  the 
compressional  waves.  Because  of  high  attenuation  of  shear  waves,  it  is  possi¬ 
ble  to  induce  them  in  a  very  limited  region  in  the  vicinity  of  the  focal  point  of 
the  ultrasound  beam,  hence  avoiding  the  influence  of  tissue  boundaries.  Shear 
parameters  of  the  tissue,  such  as  shear  modulus  and  shear  viscosity,  can  be 
calculated  by  measuring  the  amplitude  and  the  temporal  characteristics  of 
this  wave.  For  example,  the  time  required  for  the  wavefront  to  propagate 
from  one  point  to  another  can  be  used  to  calculate  the  shear- wave  speed,  and 
consequently  the  shear- wave  modulus  \x  as 

M  =  Pci  (2) 

where  p  and  ct  are  the  density  and  the  shear-wave  velocity,  respectively. 

Shear  waves  may  be  detected  optically.  In  this  method  a  laser  source 
and  a  photo  detector  are  used  to  detect  the  displacement  of  particles  due  to 
the  shear  wave  in  a  transparent  phantom  [6].  Because  this  method  requires 
a  transparent  medium,  its  application  in  vivo  is  difficult.  Phase-sensitive 
MRI  [2,6]  is  an  alternative  method  that  can  be  used  to  measure  the  3-dimen¬ 
sional  distribution  of  particle  displacement  in  a  given  direction  versus  time  in 
a  material.  Figure  4  illustrates  a  simplified  system  for  shear-wave  elasticity 
measurement  using  MRI.  In  an  experiment  presented  in  [6],  an  ultrasound 
pulse  of  3.6-ms  duration  produced  by  a  70-mm  diameter  transducer  focused 
at  100  mm  was  transmitted  within  a  cylindrical  rubber  phantom.  The  dis¬ 
placement  was  measured  by  2.0-T  MRI  system  at  two  different  times  after 
the  acoustic  pulse  was  applied.  The  position  of  the  peak  displacement  at  these 
time  points  was  used  to  estimate  the  shear  velocity,  which  was  shown  to  be 
consistent  with  the  independently  measured  value.  Shear  waves  can  also  be 
detected  by  Doppler  ultrasound  [5].  It  has  been  shown  [7]  that  to  achieve 
the  appreciable  displacement  needed  for  Doppler  detection,  most  soft  tissues 
require  high  ultrasound  intensities  which  might  be  beyond  the  safe  limit. 

3.3  Vibro-Acoustography 

This  technique  produces  a  map  of  the  mechanical  response  of  an  object  to 
a  dynamic  force  applied  at  each  point.  The  method  utilizes  an  ultrasound 
radiation  force  to  remotely  exert  a  localized  oscillating  stress  field  at  a  desired 
frequency  within  (or  on  the  surface  of)  an  object,  and  records  the  resultant 
acoustic  response  [8,9].  Figure  5  illustrates  the  principle  of  this  method.  This 
acoustic  response,  which  is  normally  in  the  low  kHz  range,  is  a  function  of 
the  viscoelastic  properties  of  the  object  and  can  be  used  to  produce  an  image 
of  the  object. 
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Fig.  4.  Principle  of  shear-wave  elastography.  The  focused  ultrasound  beam  is  ab¬ 
sorbed  by  the  object,  resulting  a  localized  radiation  stress,  which  squeezes  the  elastic 
object.  This  deformation  propagates  in  the  form  of  shear  waves  through  the  object. 
The  phase-sensitive  MRI  machine  detects  the  spatial  distribution  of  the  resulting 
displacement  in  the  object 


Focused 

ultrasound 


To  Monitor 


Vibration 


Fig.  5.  Principle  of  vibro-acoustography.  The  ultrasound  beam  is  used  to  produce 
a  localized  oscillatory  radiation  stress  to  vibrate  the  object.  Vibration  of  the  object 
produces  an  acoustic  emission  field  in  the  medium.  This  field,  which  is  a  function 
of  the  viscoelastic  properties  of  the  object,  is  detected  by  a  sensitive  hydrophone 

3.3.1  Method 

This  method  operates  on  an  oscillating  radiation  force  to  vibrate  the  tissue 
and  consequently  produce  an  acoustic  emission  from  the  object.  To  produce 
an  oscillating  radiation  force  the  intensity  of  the  incident  ultrasound  must  be 
amplitude  modulated  at  the  desired  low  frequency.  Using  a  single  amplitude- 
modulated  beam  seems  to  be  the  simplest  means  to  attain  this  purpose. 
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However,  such  a  beam  could  exert  a  radiation  force  on  any  object  that  is 
present  along  the  beam  path,  producing  undesirable  acoustic  emission.  To 
confine  the  radiation  stress  to  the  desired  region,  we  use  two  unmodulated 
cw  beams  at  slightly  different  frequency,  propagating  along  separate  paths. 
The  beams  are  positioned  to  cross  each  other  at  their  respective  foci,  and 
thus  produce  a  modulated  field  at  a  confined,  small,  cross-sectional  region. 

3.3.2  Theory 

In  the  beamforming  method  to  be  described,  the  amplitude-modulated  field  is 
obtained  by  the  interference  of  two  unmodulated  ultrasound  beams.  The  main 
advantage  of  this  approach  is  that  the  interference  volume  size  can  be  limited 
to  a  small  region.  Radially  symmetric  interfering  beams  are  obtained  when 
two  coaxial,  confocal  transducers  are  used  [9].  For  this  purpose,  elements  of 
a  two-element  spherically  focused  annular  array  (consisting  of  a  central  disc 
with  radius  a\  and  an  outer  ring  with  the  inner  radius  of  a2  and  outer  radius 
of  a,2 )  are  excited  by  separate  cw  signals  at  frequencies  =  ljo  -  Acj/2  and 
(jj2  —  ujq  T  A  £j/2,  respectively.  We  assume  that  the  beams  are  propagating  in 
the  +z-direction  with  the  joint  focal  point  at  z  —  0,  as  shown  in  Fig.  6.  The 
resultant  field  on  the  z  =  0  plane  may  be  written  as  follows: 

pit )  =  Pi(r)  cos  [(Jit  +  (r)]  +  P2(r)  cos  [oj2t  +  t\)2 (r)\  ,  (3) 

where  r  =  \J x 2  +  y 2  is  the  radial  distance.  The  amplitude  functions 
are  [17] 


x 


Fig.  6.  The  coordinate  system  describing  the  position  of  the  ultrasound  beam  and 
the  observation  point  relative  to  the  object.  The  beam  travels  in  the  ^-direction 
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P2{r)  =  pdJo2- 


a^jinc 


■  a2  jinc 


where  U0i  is  the  velocity  amplitude  at  the  *th  transducer  element  surface,  and 
\i  =  2tt /uj  for  i  —  1.2  are  the  ultrasound  wavelengths.  The  phase  functions, 


for  i  =  1.2,  are  conveniently  set  to  be  zero  at  the  origin.  Also,  jinc(X)  = 
J\  {2ixX)/'nX,  where  J\  (. . .)  is  the  first-order  Bessel  function  of  the  first  kind. 

The  instantaneous  energy  is  E  =  p2(t)/pc2.  Replacing  p(t.)  from  (3),  this 
energy  will  have  a  time-independent  component,  a  component  at  the  differ¬ 
ence  frequency  Au>  =  u>2  -  u>i  which  results  from  the  cross-product  of  the 
two  pressure  fields,  and  high-frequency  components  at  lo\  and  lj2  and  their 
harmonics.  The  energy  component  at  the  difference  frequency  is 

eAu(t.)  =  cos [Auf-A (/>(?•„)]  ,  (7) 


•cos [A ojt.  -  A V'(ro)]  , 


where  A^(r0)  =  (r0)  -  Vh  (»'<>)• 

Now,  we  define  a  unit  point  target  with  an  area  of  dxdy  at  position 
0ro?Vo)  on  the  focal  plane,  and  with  a  drag  coefficient  dr(x{),y0)  such  that 
dT(xO'.  yo)dxdy  =  1  on  the  target  and  zero  elsewhere.  This  equation  is  merely 
used  as  a  mathematical  model.  In  this  case,  the  projected  area  can  be  con¬ 
sidered  to  be  S  =  dxdy.  Therefore',  if  the  projected  area  is  unity,  then 
dr(x.y)  =  1,  which  corresponds  to  a  totally  absorptive  object. 

Referring  to  (1),  and  replacing  dvS  with  unity  and  (E)  with  eAu,(t)  of  (7), 
we  can  write  the  low-frequency  component  of  the  radiation  force  on  the  unit 
point  target  as 

Ja  wfa’th  .Vo;  t)  =  ^2^1  (ro)P2(nj)  COS  [A  Lot  +  A  V>(ro)] ,  (8) 

where  arguments  x.q  and  jjq  are  added  to  denote  the  position  of  the  noint 

,  i  i  /  o  !  o  ^  ^  1 


target,  and  r0  =  \/.rg  +  j/g.  Referring  to  (4)  and  (5),  the  complex  amplitude 
of  the  stress  field  can  be  found  as 


Fau>  (xo  ,  yo )  —  pU o  i  Uq2  — — j  inc 
X  iZo 


7T  a2 

— — jmc 

A2Z{) 


where  A  A  =  2txc/Auj  is  the  wavelength  associated  with  Aw.  The  above 
equation  indicates  that  the  radiation  stress  is  concentrated  at  the  focal  point 
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and  decays  quickly  as  jinc(. .  .)2.  It  should  be  noted  that  because  /aw(xo,  J/oi  t) 
and  Fau(xo,Vo)  are  defined  per  unit  point  target,  they  have  the  dimension 
of  force  per  unit  area,  or  equivalently  the  unit  of  stress.  As  an  example,  we 
consider  a  confocal  transducer  with  dimensions  a\  —  14.8  mm,  a'2  =  16.8  mm, 
a2  =  22.5  mm,  and  the  focal  length  70  mm.  Also,  we  assume  that  the  center 
frequency  is  3  MHz,  and  the  difference  frequency  is  7.3  kHz.  The  radiation 
stress  at  the  focal  plane  of  this  transducer  is  plotted  in  Fig.  7. 

In  medical  applications  the  maximum  ultrasonic  intensity  is  regulated 
for  safety  reasons.  It  is  therefore  useful  to  write  the  stress  field  in  terms  of 
the  peak  ultrasonic  intensity  at  the  focal  point.  The  long-term  average  of 
ultrasonic  intensity  at  the  focal  point  can  be  written  as 


m 


pm±im 

2  pc 


(10) 


where  Pi(0)  and  P2( 0)  can  be  found  from  (4)  and  (5)  as 


Assuming  Uoi  =  U02  =  do,  we  can  write  the  focal  plane  stress  field  in  terms 
of  1(0): 


FAu,(xo,yo) 


<D 


•o 


a 

E 
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Fig.  7.  Normalized  stress  field  of  the  vibro-acoustography  system  with  a  3-MHz 
confocal  transducer 
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The  first  fraction  in  the  above  equation  represents  the  static  radiation  force 
produced  by  the  two  beams  on  a  total  absorber,  the  second  fraction  is  a  con¬ 
stant  factor,  and  the  rest  represents  the  spatial  distribution  of  the  stress  field 
on  the  focal  plane.  If  Aw  o>2,  then  we  may  replace  Aj  and  A2  with  im¬ 
age  spatialAo,  and  simplify  the  expression  for  the  stress  field.  Under  these 
conditions,  the  stress  at  the  focal  point  is 


*A  u,(0, 0) 


2/(0)  a?  (o|  -  «22) 

C  «1  +  («2  -  «22)2 


(13) 


The  fraction  on  the  right  represents  the  effect  of  transducer  dimension  on 
the  stress  field.  For  the  transducer  used  in  the  previous  example,  the  value 
of  this  fraction  is  0.4999.  [This  fraction  can  be  also  calculated  for  a  single¬ 
element  transducer  of  the  same  diameter  excited  by  an  amplitude-modulated 
signal.  For  this  purpose,  we  may  let  a  1  =  n 2  and  a. 2  =  0.  For  these  values, 
the  last  fraction  in  (13)  is  O.5.]  The  resulting  radiation  stress  (assuming  c  = 
1500 m/s  for  water)  is  FAu)( 0,0)  =  6.67  x  10-4/(0)N/m2.  Now,  letting  /( 0)  = 
7200  W/  m2,  which  is  the  intensity  suggested  by  the  FDA  for  safe  in  vivo 
applications,  the  resulting  stress  at  the  focal  point  is  FAu)(0,0)  =  4.80 N/m2. 
Referring  to  Fig.  7,  we  note  that  this  stress  field  is  applied  only  in  a  small 
region  around  the  focal  point  to  the  object. 


Acoustic  Emission  To  explain  the  acoustic  emission  we  consider  an  “ob¬ 
ject  within  an  homogeneous  infinite  medium.  This  model  allows  us  to  sep¬ 
arate  the  roles  played  by  the  parameters  of  the  object  and  the  surrounding 
medium.  Also,  this  model  can  be  fitted  to  various  applications.  When  an  oscil¬ 
lating  stress  field  is  applied  to  the  object,  the  object  vibrates  at  the  frequency 
of  the  stress  field.  Vibrational  energy  of  the  object  is  partly  transferred  to  the 
suriounding  medium,  resulting  in  an  acoustic  emission  field.  Here,  we  calcu¬ 
late  the  acoustic  emission  of  an  object  subjected  to  cyclic  radiation  stress. 

To  explain  the  physics  in  an  analytical  form,  we  consider  a  flat  plate 
facing  the  beam.  Here,  we  assume  that  the  vibrating  object,  has  a  circular 
cross-section  of  radius  b  and  uniformly  vibrates  back  and  forth  like  a  piston. 
We  also  consider  an  area  S  <  ~b2  of  the  piston  surface  to  be  projected 
normally  to  the  beam.  We  can  always  return  to  our  elementary  point  object 
by  reducing  the  area  of  this  disk  to  d.7;dy.  Similar  solutions  can  be  carried 
out  for  objects  of  other  forms.  The  theory  can  be  also  extended  to  include 
arbitrary  vibrating-part  shapes  and  nonuniform  displacement  of  the  object. 
The  total  radiation  force  on  this  object,  FAu,  can  be  found  by  integrating 
the  radiation  stress  over  the  area  of  the  object,.  This  force  vibrates  the  target 
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object  at  frequency  Aw.  The  steady-state  normal  velocity  amplitude  of  a 
piston  at  frequency  A  a>,  Uaw,  due  to  a  harmonic  force  Faw,  can  be  written  as 


Fa  u) 


Fac 

uj 


(14) 


where  Zaw  is  the  mechanical  impedance  of  the  object  at  Aw.  The  mechanical 
impedance  of  the  object  has  two  components,  one  resulting  from  the  inertia, 
friction,  and  the  elasticity  of  the  object  itself,  and  the  other  resulting  from 
the  loading  effect  of  the  surrounding  medium  on  the  vibrating  object.  The 
mechanical  impedance  can  be  interpreted  as  a  measure  of  object  rigidity  and 
how  much  it  yields  to  the  applied  force.  For  example,  for  a  rigid  object, 
is  high  and  hence  resisting  the  force. 

Knowing  (7aw,  we  can  calculate  the  pressure  field  it  produces  in  the 
medium.  We  assume  that  the  acoustic  emission  signal  propagates  in  a  free  and 
homogenous  medium.  The  far-field  acoustic  pressure  due  to  a  piston  source 
of  radius  b  set  in  a  planar  boundary  of  infinite  extent  is  given  by  [17]: 


Fa  u>  = 

„exp(jA ul/c)  ( 2Ji  (^sinfl)  _cosd_\ 

i  47 d  y  ^r^sintf  cos$  +  /3b  J 


(2t Tb2UAu)  ,  (15) 


where  l  is  the  distance  from  the  observation  point  to  the  center  of  the  piston, 
d  is  the  angle  between  this  line  and  the  piston  axis,  and  (3b  is  the  specific 
acoustic  admittance  of  the  boundary  surface  (the  specific  acoustic  admittance 
is  (3b  —  where  Zg,  the  acoustic  impedance  of  the  boundary,  represents 
the  ratio  between  the  pressure  and  normal  fluid  velocity  at  a  point  on  the 
object).  The  factor  two  comes  from  the  presence  of  the  boundary  wall.  It 
would  be  replaced  by  unity  if  the  boundary  wall  were  not  present  [17]. 

The  acoustic  emission  field  resulting  from  object  vibration  can  be  written 
in  terms  of  object  mechanical  impedance  by  combining  (14)  and  (15)  as 


P  -  nr2  V  exPGA^/c)  ( 2Jl  (^Sinfl)  COST?  \ 

c2  47rZ  ^^sini?  cosi9  +  /?By 

2irb2— — 
x-^—  Fau,. 


(16) 


For  wavelengths  long  compared  to  the  object  size,  i.e.,  when  bAuj/c  — >  0, 
the  term  in  the  paranthesis  approaches  a  constant,  hence  we  may  consider  the 
contents  of  the  brackets  to  be  an  object-independent  function  (the  specific 
acoustic  admittance  (3b  relates  to  the  surrounding  boundary  surface).  Under 
these  conditions,  the  bracket  contents  in  the  above  equation  represent  the 
effect  of  the  medium  on  the  acoustic  emission  field,  which  we  may  call  the 
medium  transfer  function ,  which  my  be  denoted  by 
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The  last  fraction  in  (16)  includes  ,  which  is  the  mechanical  admittance 
of  the  object  at  the  frequency  of  the  acoustic  emission  (Aw),  and  we  denote 
this  by  Vau/-  It  is  convenient  to  combine  this  term  with  the  next  term  (2rrb2) 
in  (16).  as  Q =  2tt62Yau;  =  27r/>2 /Zaw,  which  is  the  total  acoustic  outflow 
by  the  object  per  unit  force  (acoustic  outflow  is  the  volume  of  the  medium, 
e.g.,  the  fluid,  in  front  of  the  object  surface  that  is  displaced  per  unit  time 
due  to  object  motion).  Function  Q Au>  represents  the  object  characteristics  at 
the  acoustic  frequency.  We  may  thus  rewrite  (16)  in  a  more  compact  form  as 

Pauj  =  Pc2Hau{1)Qau>Fa*  *  (18) 

Equation  (18)  indicates  that  the  acoustic  emission  pressure  is  proportional 
to  the  following:  (a)  the  radiation  force,  which  itself  is  proportional  to  the 
square  of  ultrasound  pressure  and  the  ultrasound  characteristics  of  the  ob¬ 
ject,  dr.  in  the  projected  area  S';  (b)  the  acoustic  outflow  by  this  object,  Qau, 
representing  the  object  size  b  and  its  mechanical  admittance  at  the  acoustic 
frequency,  Ya^:  and  (c)  the  transfer  function  of  the  medium  at  the  acous¬ 
tic  frequency,  H Av(l)-  Note  the  difference  between  the  projection  area  S  and 
the  vibrating  area  nb2.  The  projection  area  determines  the  extent  of  the  force 
applied  to  the  object  (1).  The  vibrating  area,  however,  influences  the  total 
acoustic  outflow  in  the  medium  caused  by  object  vibration. 

3.3.3  Applications 

Evaluating  the  characteristics  of  an  object  (or  a  medium)  by  listening  to 
its  sound  is  a  traditional  approach  that  has  been  used  for  many  purposes. 
Qualitative  evaluation  of  a  crystal  glass  by  tapping  on  it  is  a  simple  example. 
Vibro-acoustography  implements  the  same  approach  but  in  microscale,  and 
in  a  wav  that  could  be  applied  to  tissues. 

Equation  (18)  presents  a  general  relationship  between  the  mechanical  pa¬ 
rameters  of  the  objects,  surrounding  medium,  and  the  acoustic  emission  field 
resulting  from  object  vibration.  By  measuring  the  acoustic  emission  field,  it 
would  be  possible,  in  principle,  to  estimate  some  of  the  object  or  medium  pa¬ 
rameters,  either  in  an  absolute  or  in  a  relative  sense.  For  example,  one  can  use 
vibro-acoustography  to  measure  the  resonance  frequency  of  an  object,  and 
from  that  information  it  is  possible  to  estimate  some  viscoelastic  parameters 
of  the  object  or  the  medium  (Fig.  8).  This  method  has  been  used  to  mea¬ 
sure  the  Young’s  modulus  of  a  metallic  rod  [18],  In  another  experiment,  by 
measuring  the  resonance  frequency  of  a  known  resonator  in  a  liquid  medium, 
the  viscosity  of  the  fluid  has  been  estimated  with  good  accuracy  [19].  These 
applications  are  not  necessarily  medical,  but  the  principle  could  be  used  for 
evaluation  of  soft  tissues,  blood,  bone,  etc. 

In  medical  applications,  one  can  use  vibro-acoustography  to  obtain  images 
of  the  human  body  for  diagnostic  purposes.  In  such  applications,  the  image 
may  not  represent  a  physical  quantity,  rather  it  provides  a  means  to  visualize 
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Fig.  8.  Vibro-acoustography  system.  The  object  is  vibrated  by  the  radiation  force 
of  the  confocal,  two-beam  transducer.  Resulting  vibrations  are  detected  by  the 
microphone  (or  hydrophone)  and  mapped  into  the  image 


object  details.  To  demonstrate  the  capability  of  vibro-acoustography,  we  used 
this  method  to  image  a  specimen  of  the  human  iliac  artery  that  included 
calcifications.  A  3-MHz  confocal  transducer  with  dimensions  a\  =  14.8  mm, 
a!2  =  16.8  mm,  and  a 2  =  22.5  mm  and  a  focal  length  of  70  mm,  was  used 
for  this  experiment.  The  difference  frequency  was  set  at  44  kHz.  Figure  9 
shows  a  photo  of  the  artery  that  is  cut  open.  The  areas  with  calcification 
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Fig.  9.  Vibro- acoustic  image  of  an  ex-vivo  human  iliac  artery  near  the  bifurcation. 
On  the  left  is  a  photo  of  the  cut-open  artery  specimen.  Calcified  areas  are  marked 
with  letter  “C”.  Vibro- acoustic  image  of  the  specimen  is  shown  on  the  right.  Calci¬ 
fied  areas  stand  out  as  bright  spots.  This  image  was  obtained  with  a  3-MHz  confocal 
transducer  as  described  in  the  text.  The  difference  frequency  was  set  at  44  kHz.  The 
gray-level  scale  is  in  arbitrary  units 


272  Mostafa  Fatemi  and  James  F.  Green  leaf 


are  marked  by  the  letter  “C” .  The  vil  no-acoustic  image  of  the  same  artery 
is  shown  on  the  right.  Calcium  deposits,  which  are  formed  in  large  (a  few 
mm  to  a  few  cm  long)  and  stiff  (compare  to  the  arterial  wall)  plates,  are 
efficient  acoustic  radiators,  producing  strong  acoustic  emission  when  they 
are  exposed  to  the  radiation  force  of  the  incident  ultrasound.  Hence,  they 
stand  out  in  the  vibro-acoustography  image.  The  vibro-acoustic  image  has 
high  resolution,  no  speckle,  and  a  high  signal-to-noise  ratio.  These  features 
allow  vibro-acoustography  to  delineate  calcium  deposits  with  high  definition 
and  contrast. 

4  Capabilities  and  Limitations 

Although  the  methods  described  above  tend  to  estimate  the  mechanical  prop¬ 
erties  of  the  object,  the  quantities  they  measure  are  different.  The  quantity 
measured  in  the  transient  method  is  somewhat  arbitrary,  and  not  a  direct 
measure  of  any  physical  quantity.  The  merit  of  this  method  is  that  tissue 
hardness  is  measured  from  relative  displacement  values,  without  the  need  to 
measure  the  applied  force.  The  shear-wave  method  directly  measures  a  phys¬ 
ical  parameter  of  the  material  in  an  absolute  sense.  The  acoustic  emission 
field  measured  in  vibro-acoustography  is  a  function  of  several  physical  pa¬ 
rameters  that  relate  to  the  object  and  the  surrounding  medium.  It  is  possible 
to  measure  a  single  parameter  only  if  one  has  enough  knowledge  about  other 
parameters  involved.  For  example,  one  can  measure  the  shear  viscosity  of  a 
liquid  using  a  known  resonator  as  an  object  in  the  liquid  [19],  or,  for  a  given 
the  geometry,  the  Young’s  modulus  of  a  metallic  object  can  bo  measured  [18] 
by  vibro-acoustography. 

Detection  sensitivity  is  critical  to  the  success  of  any  of  the  above  methods. 
Greater  ultrasound  intensity  would  be  needed  if  the  detection  sensitivity  is 
poor.  As  a  result,  the  peak  ultrasound  intensity  in  the  focal  area  may  be 
limited  by  the  safe  in  vivo  limit.  Ultrasound  Doppler  methods,  at  conventional 
ultrasound  frequencies,  probably  can  detect  displacements  on  the  order  of  a 
few  micrometers.  The  sensitivity  of  the  MRT  technique  is  on  the  order  of 
lOOnrn  [2].  Vibro-acoustography  has  been  shown  to  detect  motions  as  small 
as  a  few  nanometers  [8].  The  high  sensitivity  of  this  method  is  a  result  of 
the  fact  that  small  motions  of  the  object  can  produce  an  acoustic  emission 
pressure  field  that  is  easily  detectable  by  a  sensitive  hydrophone. 

Image  spatial  resolution  is  an  important  parameter  when  assessing  the 
capability  of  an  imaging  method  in  delineating  fine  st  ructures  in  the  object. 
Vibro-acoustography  can  be  used  for  high-quality  imaging.  The  spatial  reso¬ 
lution  is  proportional  to  the  width  of  the  stress-field  main  lobe.  For  example, 
for  the  3-MHz  ultrasound  transducer  presented  in  Fig.  7,  the  spatial  reso¬ 
lution  is  about  700pm.  Imaging  complex  objects  using  the  transient  or  the 
shear-wave  methods  has  not  been  fully  explored  in  the  literature. 
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Performance  of  each  method  is  partly  influenced  by  viscoelastic  properties 
of  the  material.  Both  the  transient  and  shear-wave  methods  work  best  if  the 
material  under  test  supports  the  shear  wave,  has  enough  attenuation  to  al¬ 
low  the  build-up  of  enough  radiation  force,  and  is  compliant  enough  to  allow 
appreciable  displacement.  Application  of  these  methods  for  hard  materials, 
such  as  bone  and  calcifications,  would  be  difficult  or  impractical  because  (a) 
stiff  materials  have  low  compliance,  hence  their  displacement  in  response  to 
the  force  is  relatively  small  and  difficult  to  detect  and  (b)  the  shear  wave 
travels  at  high  speed  in  such  materials.  In  vibro-acoustography  hard  materi¬ 
als  produce  structural  vibrations  which  are  often  stronger  than  those  of  soft 
tissue.  Such  vibrations  result  in  strong  and  easily  detectable  acoustic  emis¬ 
sions.  The  bright  appearance  of  the  arterial  calcifications  in  the  experiment 
presented  in  [8]  is  a  demonstration  of  this  phenomenon.  Vibro-acoustography 
can  be  used  to  detect  particles  in  materials  that  do  not  support  shear  waves, 
for  example,  detecting  gas  bubbles  in  liquids.  The  acoustic  emission  resulting 
from  particle  vibration  includes  a  compressional-wave  component  that  can 
travel  in  the  surrounding  medium  including  the  liquid.  Tissue  attenuation 
for  compressional  waves  at  low  frequencies  is  small;  hence  such  waves  can  be 
detected  by  the  hydrophone  from  a  distance.  Shear  waves  in  a  liquid  medium 
decay  very  rapidly  and  are  difficult  to  detect  by  most  methods. 

The  ultimate  goal  of  viscoelastic  parameter  imaging  methods  is  to  mea¬ 
sure  these  parameters  in  the  human  body.  In  vivo  application  of  the  transient 
method  is  not  clear  because  it  requires  high  ultrasound  power  to  produce  a 
detectable  displacement.  Application  of  the  shear- wave  method  in  the  hu¬ 
man  body  largely  depends  on  the  detection  method  used.  Optical  methods 
are  not  likely  to  be  usable  in  vivo.  MRI  detection  methods  are  generally 
complex,  especially  when  the  technique  is  to  be  linked  to  the  ultrasound  sys¬ 
tem  in  the  magnetic  field.  Ultrasound  Doppler  is  a  more  practical  detection 
method  for  shear- wave  detection;  however  its  sensitivity  may  not  be  sufficient 
for  ultrasound  intensities  within  the  safe  limit  [7] .  Vibro-acoustography  uses 
a  hydrophone  for  detection  which  is  simple  to  operate  in  clinical  settings. 
It  however  requires  an  acoustically  quiet  environment  for  proper  detection. 
The  biological  noise  spectrum  seems  to  be  below  1  kHz,  which  can  be  easily 
filtered  out  if  the  operation  frequency  is  above  this  limit  [8]. 

Before  any  of  the  methods  discussed  here  are  used  in  the  human  body, 
one  must  be  certain  that  the  ultrasound  exposure  does  not  harm  the  sub¬ 
ject.  Vibro-acoustography  uses  either  cw  or  tone-burst  ultrasound.  Therefore, 
the  continuous  ultrasound  intensity  (spatial  peak  temporal  average  intensity) 
limit  must  be  observed  for  the  safety  reasons.  This  limit  according  to  the  FDA 
is  720  mW/  cm2.  However,  because  of  the  high  sensitivity  of  hydrophone  de¬ 
tectors,  it  is  possible  to  detect  very  small  levels  of  the  acoustic  emission  while 
using  low  transmittance  power.  For  example,  it  is  shown  that  a  submillimeter 
object  can  be  detected  at  ultrasound  intensities  far  below  the  FDA  limit  [9]. 
The  shear- wave  method  presented  in  [6]  uses  high-intensity  pulsed  ultrasound 
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for  excitation  and  either  the  optical  or  MRI  methods  for  detection.  It  is  ar¬ 
gued  that,  although  the  intensity  is  high,  the  total  exposure  is  within  the 
FDA  limits.  However,  in  vivo  use  of  such  detection  methods  would  be  diffi¬ 
cult.  Using  Doppler  ultrasound  for  detection  is  less  sensitive;  hence  it  requires 
higher  peak  ultrasound  intensity  to  produce  enough  radiation  force  on  the 
tissue.  This  intensity  level  may  exceed  the  safe  spatial  peak  pulse  average 
intensity  limit. 


5  Summary 

The  radiation  force  of  ultrasound  is  a  noninvasive  means  for  introducing  a 
highly  localized  vibrating  force  inside  tissue.  The  three  methods  described 
here  utilize  this  tool  to  evaluate  the  dynamic  viscoelastic  properties  of  tissue. 
The  transient  method  measures  minute  transient  tissue  deformation  versus 
time.  The  shear-wave  method  measures  the  amplitude  and  velocity  of  the 
shear  waves  resulting  from  excitation  by  a  radiation  force.  Finally,  vibro- 
acoustography  measures  the  acoustic  field  resulting  from  object  vibration  at 
a  specified  frequency. 
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Appendix  7 

Contrast  Coherent  Imaging  (CCI)  at  low  mechanical  index  was  used.  On  a  semiquantitative  scale  the  capability,  to :  . 

structures  and  possible  modifications  of  parenchyma  during  the  basal  arterial,  portal,  and  late  phases  using  color  or  power  Dopp  er  wi 
and  without  echocontrast  administration,  was  evaluated.  RESULTS  AND  CONCLUSIONS:  Color/power  Doppler ds/dely  used  for  t  e 
study  of  focal  liver  lesions.  Its  limitations  in  terms  of  sensitivity  prevents  the  possibility  of  showing  the  microcire&lation.  This  study 
demostrates  that  low  MI-CCI  and  contrast  enhancement  allow  the  visualization  of  both  macrovascular  andpdrenchymal  circulation  c  ar- 

acteristics. 
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sound  to  induce  low-frequency  elastic  waves  in  tissue  and  detects  its  response.  These  waves  are  detected  by  an  audio  hydrophone  and  can 
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Vibro-acoustic  images  of  tissue  samples  provided  clear  and  well  defined  display  of  microcalcifications.  The  number  and  positio 


S:25 


J  Ultrasound  Med  20:S1-S154,  March  2001 


American  Institute  of  Ultrasound  in  Medicine  Proceedings 


microcalcifications  matched  with  those  seen  in  mammography.  Microcalcifications  as  small  as  100  microns  in  diameter  could  he  detected 
hy  vibro-acoustography.  The  presence  of  microcalcifications  was  validated  by  histological  study  of  the  samples.  CONCLUSIONS:  Vibro- 
acoustography  can  be  used  for  imaging  small  breast  microcalcifications.  The  non-invasive  nature  of  vibro-acoustography  promises  in  vivo 
applications,  especially  in  cases  where  mammography  cannot  be  used  either  because  of  high  breast  density  or  hazards  of  x-rav  to  the 
patient. 
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Abstract  Over  the  past  few  years  there  has  been  an  increasing  interest  in  using  the  radiation  force 
of  ultrasound  for  evaluating,  characterizing  and  imaging  biological  tissues.  Of  particular  interest 
are  those  methods  that  measure  the  dynamic  properties  of  tissue  at  low  frequencies.  In  this  paper 
we  present  dynamic  radiation  force  methods  for  probing  tissue  as  a  new  field,  discuss  the  inter¬ 
relationship  of  several  methods  within  this  field  and  compare  their  features.  The  techniques  in  this 
field  can  be  categorized  into  three  groups:  transient  methods,  shear-wave  measurement  methods 
and  a  recently  developed  method  called  vibro-acoustography.  The  last  method  is  the  focus  of  this 
paper.  After  briefly  describing  the  key  concepts  of  the  first  two  methods,  we  will  present  a  detailed 
description  of  vibro-acoustography.  Finally,  we  will  compare  the  capabilities  and  limitations  of 
these  methods. 


1.  Introduction 

It  is  well  known  that  changes  in  the  elasticity  of  soft  tissues  are  often  related  to  pathology. 
Traditionally,  physicians  use  palpation  as  a  simple  method  for  estimating  the  mechanical 
properties  of  tissue.  In  this  method,  a  static  force  is  applied  and  a  crude  estimation  of  tissue 
elasticity  is  obtained  through  the  sense  of  touch.  In  palpation,  the  force  is  applied  on  the  body 
surface  and  the  result  is  a  collective  response  of  all  the  tissues  below.  Clinicians  can  sense 
abnormalities  if  the  response  to  palpation  of  the  suspicious  tissue  is  sufficiently  different  from 
that  of  normal  tissue.  However,  if  the  abnormality  lies  deep  in  the  body,  or  if  is  too  small  to 
be  resolved  by  touch,  then  the  palpation  method  fails.  The  dynamic  response  of  soft  tissue 
to  a  force  is  also  valuable  in  medical  diagnosis.  For  instance,  rebound  of  tissue  upon  sudden 
release  of  pressure  exerted  by  the  physician’s  finger  on  the  skin  provides  useful  diagnostic 
information  about  the  issue. 

Quantitative  measurement  of  the  mechanical  properties  of  tissues  and  their  display  in 
a  raster  format  is  the  aim  of  a  class  of  techniques  generally  called  elasticity  imaging,  or 
elastography.  The  general  approach  is  to  measure  tissue  motion  caused  by  an  external  (or,  in 
some  methods,  internal)  force  or  displacement  and  use  it  to  reconstruct  the  elastic  parameters 
of  the  tissue.  The  excitation  stress  can  be  either  static  or  dynamic  (vibration).  Dynamic 
excitation  is  of  particular  interest  because  it  provides  more  comprehensive  information  about 
tissue  properties  in  a  spectrum  of  frequencies;  alternatively,  the  transient  behaviour  of  the  tissue 
could  be  deduced  from  the  measurements.  In  most  elasticity  imaging  methods,  ultrasound  is 
used  to  detect  the  motion  or  displacement  resulting  from  the  applied  stress.  Magnetic  resonance 
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elastography  is  a  recently  developed  method  (Muthupillai  etal  1995)  that  employs  a  mechanical 
actuator  to  vibrate  the  body  surface  and  then  measures  the  strain  waves  with  a  phase  sensitive 
magnetic  resonance  imaging  (MRI)  machine. 

The  majority  of  elasticity  imaging  methods  are  based  on  an  external  source  of  force. 
In  these  methods  the  object  is  pressed  by  a  known  amount  of  force  or  displacement,  and  the 
resulting  internal  deformations  are  measured  by  means  of  pulse-echo  ultrasound.  The  elasticity 
of  the  region  of  interest  is  then  calculated  based  on  the  resulting  deformation  in  relation  to  the 
magnitude  of  the  applied  force  (or  displacement).  Normally,  the  region  of  interest  rests  deep 
in  the  body  and  away  from  the  source  of  the  force.  The  force  that  is  actually  exerted  on  the 
region  of  interest  depends  on  the  elastic  properties  of  the  tissues  located  between  the  source 
and  the  region  of  interest.  Hence,  the  deformation,  and  the  estimated  elasticity  of  the  region 
of  interest,  are  subject  to  variables  of  the  intervening  tissues. 

An  alternative  strategy  is  to  apply  a  localized  stress  directly  in  the  region  of  interest.  One 
way  to  accomplish  this  is  to  use  the  radiation  pressure  of  ultrasound.  Acoustic  radiation  force 
is  the  time  average  force  exerted  by  an  acoustic  field  on  an  object.  This  force  is  produced  by  a 
change  in  the  energy  density  of  an  incident  acoustic  field  (Chu  and  Apfel  1982),  for  example, 
due  to  absorption  or  reflection.  The  use  of  ultrasound  radiation  force  for  evaluating  tissue 
properties  has  several  benefits,  for  example: 

(a)  Acoustic  (ultrasound)  energy  is  a  non-invasive  means  of  exerting  force. 

(b)  Existing  ultrasound  technology  and  devices  can  be  readily  modified  for  this  purpose,  thus 
eliminating  the  need  for  developing  a  new  technology. 

(c)  Radiation  force  can  be  generated  remotely  inside  tissue  without  disturbing  its  superficial 
layers. 

(d)  The  radiation  stress  field  can  be  highly  localized,  thus  allowing  for  precise  positioning  of 
the  excitation  point. 

(e)  Radiation  force  can  be  produced  in  a  wide  range  of  frequencies  or  temporal  shapes. 

These  features  make  radiation  force  methods  more  attractive  than  other,  mostly  mechanical, 
excitation  methods  used  in  elasticity  imaging. 

Tissue  probing  with  the  radiation  force  of  ultrasound  can  be  accomplished  by  a  variety  of 
techniques,  depending  on  the  excitation  and  detection  methods  used.  As  in  elasticity  imaging 
methods  with  mechanical  excitation,  the  radiation  force  methods  can  use  either  a  static  or 
dynamic  stress. 

It  is  the  opinion  of  the  authors  that  using  a  dynamic  radiation  force  to  remotely  probe 
tissue  has  certain  unique  characteristics  and  capabilities  that  can  open  a  new  chapter  in  the 
field  of  tissue  characterization  and  imaging.  We  believe  that  the  techniques  based  on  this 
phenomenon  have  the  potential  to  provide  new  information  about  tissue  that  has  never  been 
available  from  other  methods.  We  also  believe  that  it  is  important  to  recognize  and  treat  this 
set  of  techniques  as  a  new  field.  In  particular,  it  is  insightful  to  set  this  new  field  apart  from 
conventional  ultrasound  tissue  characterization/imaging.  A  major  difference  being  the  fact 
that  using  the  dynamic  radiation  stress  allows  one  to  analyse  the  object  based  on  its  structural 
vibration  properties  as  opposed  to  its  ultrasonic  parameters. 

Our  intention  in  writing  this  paper  is  to  present  the  dynamic  radiation  force  methods  for 
probing  tissue  as  a  new  field  by  discussing  the  inter-relationship  of  several  methods  within  this 
field  and  comparing  their  features.  To  the  best  of  our  knowledge,  such  a  discussion  has  not 
been  presented  in  the  literature  before. 

The  dynamic  radiation  force  methods  may  be  categorized  as: 

(a)  Transient  methods,  where  an  impulsive  radiation  force  is  used  and  the  transient  response 
of  the  tissue  is  detected  by  Doppler  ultrasound  (Sugimoto  et  al  1990). 
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(b)  Shear-wave  methods,  where  an  impulsive  or  oscillating  radiation  is  applied  to  the  tissue 
and  the  resulting  shear  wave  is  detected  by  ultrasound  or  other  methods  (Andreev  et  al 
1997,  Sarvazyan  etal  1998,  Walker  1999). 

(c)  Vibro-acoustography,  a  method  recently  developed  by  the  authors,  where  a  localized 
oscillating  radiation  force  is  applied  to  the  tissue  and  the  acoustic  response  of  the  tissue 
is  detected  by  a  hydrophone  or  microphone  (Fatemi  and  Greenleaf  1998,  1999a). 

In  this  paper  we  systematically  discuss  the  features  and  capabilities  that  are  offered  by  the 
dynamic  radiation  force  techniques.  We  describe  the  general  approach  used  in  the  transient, 
shear-wave  and  vibro-acoustography  methods.  We  pay  particular  attention  to  the  last  method, 
present  its  theory  and  discuss  its  features.  Then,  we  compare  the  capabilities  and  limitation  of 
all  three  methods  and  discuss  their  possible  applications. 


2,  Radiation  force  methods 


The  acoustic  radiation  force  is  a  time  average  force  exerted  by  an  acoustic  field  on  an  object. 
This  force  is  an  example  of  a  universal  phenomenon  in  any  wave  motion  that  introduces  some 
type  of  unidirectional  force  on  absorbing  or  reflecting  targets  in  the  wave  path.  For  a  review 
of  this  topic  the  reader  may  refer  to  Chu  and  Apfel  (1982). 

Consider  a  plane  ultrasound  beam  interacting  with  a  planar  object  of  zero  thickness  and 
arbitrary  shape  and  boundary  impedance  that  scatters  and  absorbs.  The  radiation  force  vector, 
F,  arising  from  this  interaction  has  a  component  in  the  beam  direction  and  another  transverse 
to  it.  The  magnitude  of  this  force,  is  proportional  to  the  average  energy  density  of  the  incident 
wave  ( E )  at  the  object,  where  (>  represents  the  time  average,  and  S  the  area  of  the  projected 
portion  of  the  object  (Westervelt  1951): 

F  =  drS{E).  (1) 


Here  dr  is  the  vector  drag  coefficient  with  a  component  in  the  incident  beam  direction  and 
another  transverse  to  it.  The  coefficient  dr  is  defined  per  unit  incident  energy  density  and  unit 
projected  area.  For  a  planar  object,  the  magnitude  of  dr  is  numerically  equal  to  the  force  on 
the  unit  area  of  the  object  per  unit  energy  density.  Physically,  the  drag  coefficient  represents 
the  scattering  and  absorbing  properties  of  the  object,  and  is  given  by  (Westervelt  1951) 


i,  =pi(na  +  rij  -  j  yconctsds'j+q^  J  ysinofjdS 


(2) 


where  p  and  q  are  the  unit  vectors  in  the  beam  direction  and  normal  to  it  respectively.  The 
quantities  na  and  ns  are  the  total  absorbed  and  scattered  powers  respectively,  and  y  is  the 
scattered  intensity,  all  expressed  per  unit  incident  intensity.  Also,  as  is  the  angle  between  the 
incident  and  the  scattered  intensity  and  dS  is  the  scalar  area  element.  The  drag  coefficient 
can  also  be  interpreted  as  the  ratio  of  the  radiation  force  magnitude  on  a  given  object  to  the 
corresponding  value  if  the  object  were  replaced  by  a  totally  absorbing  object  of  similar  size. 
For  simplicity,  we  assume  a  planar  objected  oriented  perpendicular  to  the  beam  axis.  In  this 
case,  the  transverse  component  vanishes,  thus  the  drag  coefficient  (force)  will  have  only  a 
component  normal  to  the  target  surface  which  we  denote  by  scalar  dr(F ). 

To  produce  a  time-varying  radiation  force,  the  intensity  of  the  incident  beam  can  be  shaped 
in  various  ways.  For  example,  a  short  ultrasound  pulse  can  produce  a  transient  pulsed  radiation 
force,  and  a  sinusoidally  modulated  beam  can  result  in  a  sinusoidally  varying  force. 
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2.1.  Transient  methods 

In  the  transient  methods  the  radiation  force  of  ultrasound  is  used  to  make  a  minute  deformation 
in  the  tissue.  The  transient  recoil  of  the  tissue  resulting  from  this  deformation  is  measured  and 
used  for  the  evaluation  of  tissue  elastic  properties. 

A  method  for  measuring  tissue  hardness,  presented  by  Sugimoto  et  al  (1990),  uses  the 
radiation  force  of  a  single  focused  ultrasound  beam.  Ideally,  hardness  may  be  represented  by 
the  spring  constant  of  the  object,  which  is  the  ratio  of  the  applied  force  to  the  displacement. 
This  method  relies  on  evaluating  the  hardness  based  on  changes  in  relative  values.  In  this 
method  an  ultrasound  pulse  is  used  to  generate  a  short-duration  radiation  force  which  produces 
localized  deformation  of  the  tissue.  Immediately  after  the  force  pulse,  the  resulting  transient 
deformation  of  tissue  is  measured  as  a  function  of  time  with  Doppler  ultrasound  using  a 
separate  transducer.  The  deformation  includes  an  initial  rapid  squeeze  of  the  tissue,  followed 
by  a  relaxation  and  possibly  a  rebound.  This  displacement  is  a  function  of  tissue  viscoelastic 
parameters,  as  well  as  the  applied  force.  Sugimoto  et  al  (1990)  argue  that  because  measuring 
the  internal  radiation  force  in  an  absolute  sense  in  vivo  is  difficult,  it  is  advantageous  to  derive 
a  quantity  that  is  representative  of  relative  tissue  hardness.  To  derive  a  single  relative  quantity 
from  the  deformation  data,  the  relaxation  part  of  the  function  is  approximated  by  a  sum  of 
several  exponential  curves,  and  the  sum  of  the  first-order  derivatives  of  such  an  exponential  is 
calculated.  It  is  shown  that  this  quantity  is  correlated  with  the  spring  constant  of  the  tissue  and 
thus  may  be  used  as  a  measure  of  tissue  hardness. 


2.2.  Shear^ wave  methods 

Shear  modulus  is  related  to  the  hardness  or  elasticity  of  the  material.  It  is  known  that  the  shear 
modulus  of  various  soft  tissues  ranges  over  several  orders  of  magnitude,  while  bulk  modulus,  a 
parameter  that  is  associated  with  the  conventional  pulse-echo  ultrasound  compressional  wave 
speed,  varies  significantly  less  than  an  order  of  magnitude  (Goss  et  al  1978,  Frizzell  and 
Carstensen  1976).  These  features  indicate  that  shear  modulus  may  be  a  better  parameter  for 
tissue  characterization  than  bulk  modulus.  One  way  to  induce  localized  shear  waves  inside 
tissue  is  to  use  the  radiation  force  of  focused  ultrasound  (Andreev  et  al  1997).  In  a  method 
called  shear- wave  elasticity  imaging  (SEWI)  (Sarvazyan  et  al  1998),  an  amplitude-modulated 
single-focused  ultrasound  beam  is  used  to  induce  a  localized  radiation  stress  inside  the  soft 
tissue.  Localization  of  the  stress  field  is  the  key  to  success  of  the  method.  To  achieve  a 
high  degree  of  localization,  the  method  used  a  focused  ultrasound  beam.  It  is  shown  that  the 
radiation  stress  exerted  within  a  dissipative  medium  peaks  near  the  focal  region  of  a  highly 
focused  transducer.  It  has  also  been  suggested  that  localization  can  be  improved  by  designing 
the  transducer  and  selecting  the  beam  parameters  such  that  a  nonlinear  shock  wave  is  produced 
in  the  focal  region,  increasing  the  magnitude  of  the  stress  field  in  the  vicinity  of  the  focal  region, 
thus  augmenting  localization  of  the  stress  field  (Sarvazyan  etal  1998). 

Modulation  of  the  ultrasound  beam  can  be  in  the  form  of  an  oscillating  wave  or  short  pulse. 
The  resulting  radiation  force  elicits  a  shear  wave  propagating  in  the  radial  direction  with  respect 
to  the  beam  axis,  with  particle  motion  parallel  to  the  beam  axis.  Shear  waves  in  soft  tissue 
travel  at  a  very  low  speed,  typically  around  a  few  metres  per  second  (Frizzell  and  Carstensen 
1976,  Sarvazyan  et  al  1998),  thus  the  corresponding  wavelength  is  much  shorter  than  that 
of  compressional  waves  for  the  same  frequency.  Shear  waves  are  also  highly  attenuated  in 
soft  tissue,  with  an  attenuation  coefficient  two  or  three  orders  of  magnitude  higher  than  that 
of  compressional  waves  (Frizzell  and  Carstensen  1976,  Sarvazyan  et  al  1998).  Because  of 
the  high  attenuation  of  shear  waves,  it  is  possible  to  induce  them  in  a  very  limited  region  in 
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the  vicinity  of  the  focal  point  of  the  ultrasound  beam,  hence  avoiding  the  influence  of  tissue 
boundaries.  In  SWEI,  a  short  ultrasound  pulse  is  delivered  to  the  tissue,  causing  a  radial  shear 
oscillation  around  the  beam  axis.  By  measuring  the  amplitude  and  the  temporal  characteristics 
of  this  wave,  shear  parameters  of  the  tissue,  such  as  shear  modulus  and  shear  viscosity,  can 
be  calculated.  For  example,  the  time  required  for  the  wave  front  to  propagate  from  one  point 
to  another  can  be  used  to  calculate  the  shear- wave  speed,  and  consequently  the  shear- wave 
modulus  fi,  as 

n  =  pc]  (3) 

where  p  and  ct  are  the  density  and  the  shear-wave  velocity  respectively. 

Shear  waves  may  be  detected  optically.  In  this  method  a  laser  source  and  a  photodetector 
are  used  to  detect  the  displacement  of  particles  due  to  the  shear  wave  in  a  transparent  phantom 
(Sarvazyan  et  al  1998).  Because  this  method  requires  a  transparent  medium,  its  application 
in  vivo  is  difficult.  Phase-sensitive  magnetic  resonance  imaging  (Muthupillai  et  al  1995, 
Sarvazyan  etal  1998)  is  an  alternative  method  that  can  be  used  to  measure  the  two-dimensional 
distribution  of  particle  displacement  in  a  given  direction  versus  time  in  a  material.  In  an 
experiment  presented  by  Sarvazyan  et  al  (1998),  an  ultrasound  pulse  of  3.6  ms  duration 
produced  by  a  70  mm  diameter  transducer  focused  at  100  mm  was  transmitted  within  a 
cylindrical  rubber  phantom.  The  displacement  was  measured  by  2.0  T  MRI  system  at  two 
different  times  after  the  acoustic  pulse  was  applied.  The  position  of  the  peak  displacement  at 
these  time  points  was  used  to  estimate  the  shear  velocity,  which  was  shown  to  be  consistent  with 
the  independently  measured  value.  Shear  waves  can  also  be  detected  by  Doppler  ultrasound 
(Andreev  et  al  1997).  It  has  been  shown  (Walker  1999)  that  to  achieve  the  appreciable 
displacement  needed  for  Doppler  detection,  most  soft  tissues  require  high  ultrasound  intensities 
which  might  be  beyond  the  recommended  maximum  intensity  level  for  diagnostic  ultrasound 
(FDA  1997). 

2.3.  Vibno-acoustography 

2.3. 1.  Method.  This  technique  produces  a  map  of  the  mechanical  response  of  an  object  to  a 
dynamic  force  applied  at  each  point.  The  method  utilizes  ultrasound  radiation  force  to  remotely 
exert  a  localized  oscillating  stress  field  at  a  desired  frequency  within  (or  on  the  surface  of)  an 
object  and  records  the  resultant  acoustic  response  (Fatemi  and  Greenleaf  1998,  1999a).  This 
acoustic  response,  which  is  normally  in  the  low  kHz  range,  is  a  function  of  the  viscoelastic 
properties  of  the  object  and  can  be  used  to  produce  an  image  of  the  object. 

To  produce  an  oscillating  radiation  force  the  intensity  of  the  incident  ultrasound  must 
be  amplitude  modulated  at  the  desired  low  frequency.  Using  a  single  amplitude  modulated 
beam  seems  to  be  the  simplest  means  to  achieve  this.  However,  such  a  beam  will  exert 
a  radiation  force  on  any  object  that  is  present  along  the  beam  path,  producing  undesirable 
acoustic  emission.  To  confine  the  oscillating  radiation  stress  to  the  desired  region,  vibro- 
acoustography  uses  two  unmodulated  CW  beams  at  slightly  different  frequencies,  propagating 
along  separate  paths.  The  beams  are  arranged  to  cross  each  other  at  their  respective  foci,  and 
thus  produce  a  modulated  field  at  a  confined,  small  cross-sectional  region.  The  single-  and 
double-beam  arrangements  are  illustrated  in  figure  1. 

2.3.2.  System  diagram  and  image  formation.  Figure  2  illustrates  a  vibro-acoustography 
system.  The  elements  of  the  annular  array  transducer  are  driven  by  two  CW  signals  at 
frequencies  fx  =  (o\/2n  and  f\  +  A/  =  (a>j  4*  Aa>)/27T.  The  object  to  be  imaged  is  placed 
at  the  joint  focal  plane  of  the  transducer  elements  (the  scanning  plane).  The  ultrasound  field 
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Coofocft)  beams  Single  beam 

Figure  1.  Generation  of  a  modulated  ultrasound  field  by  single  and  confocal  beams.  The  single 
beam  produces  a  field  that  is  modulated  throughout  the  beam  path.  This  field  can  produce  an 
oscillating  radiation  force  on  any  object  in  the  beam  path.  The  confocal  annular  array  transducer 
generates  two  CW  beams  propagating  separate  paths.  These  beams  interfere  near  the  focal  zone, 
producing  a  modulated  field  only  in  that  region.  Hence,  the  oscillating  radiation  force  can  only  be 
generated  on  objects  in  the  beam-crossing  region. 


Figure  2.  Vibro-acoustography  system.  The  two  elements  of  the  confocal  annular  array  are  excited 
by  separate  CW  signals  at  two  different  frequencies.  The  object  is  placed  at  the  focal  plane  of  the 
transducer.  The  ultrasound  field  produces  a  radiation  force  at  the  difference  frequency  A /  on  the 
object.  This  force  vibrates  the  object,  which  in  turn  produces  a  sound  field  in  the  medium.  This 
sound  field  is  received  by  an  audio  hydrophone  placed  near  the  object.  To  form  an  image,  the 
beams  are  scanned  in  a  raster  motion  across  the  scanning  place,  and  the  resulting  acoustic  emission 
is  recorded  from  each  point.  The  brightness  of  each  corresponding  point  on  the  image  plane  is 
determined  by  the  amplitude  (or  phase)  of  the  recorded  acoustic  emission. 


produces  a  localized  radiation  stress  field  at  the  focal  point  on  the  object  at  frequency  A / . 
Depending  on  the  elastic  properties  of  the  object,  the  radiation  force  may  cause  a  portion  of 
the  object,  or  the  entire  object,  to  vibrate  at  this  frequency.  The  acoustic  emission  resulting 
from  object  vibration  is  received  by  a  hydrophone  which  is  located  nearby.  Normally  the 
wavelength  of  vibration  is  large  compared  with  the  object  size,  hence  the  acoustic  emission 
field  is  almost  omnidirectional.  Therefore,  hydrophone  position  is  not  a  critical  parameter  in 
this  measurement.  The  filter  is  used  to  eliminate  noise  and  other  interfering  signals.  To  form 
an  image,  the  focal  point  of  the  transducer  is  moved  across  the  scanning  plane  on  the  object  in 
a  raster  pattern.  The  acoustic  emission  is  received  at  each  position,  and  an  image  is  formed  by 
displaying  the  amplitude  (or  phase)  of  such  signals  at  corresponding  positions  on  the  image 
plane.  The  spatial  resolution  of  this  imaging  method  is  determined  by  the  ultrasound  beam  width 
at  the  focal  plane,  which  is  normally  of  the  order  of  the  incident  ultrasound  wavelength. 

2.3.3.  Theory.  In  the  beam-forming  method  to  be  described  the  amplitude  modulated  field 
is  obtained  by  the  interference  of  two  unmodulated  ultrasound  beams.  Radially  symmetric 
interfering  beams  are  obtained  when  two  coaxial,  confocal  transducers  are  used  (Fatemi  and 
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Greenleaf  1999a).  For  this  purpose,  elements  of  a  two-element  spherically  focused  annular 
array  (consisting  of  a  central  disc  with  radius  ax  and  an  outer  ring  with  the  inner  radius  of  a'2 
and  outer  radius  of  a2  are  excited  by  separate  CW  signals  at  frequencies  a)\  =  coo  —  A<o/2  and 
a)2  =  (oq  +  Aco/2.  We  assume  that  the  beams  are  propagating  in  the  +z  direction  with  the  joint 
focal  point  at  z  =  0.  The  resultant  field  on  the  z  =  0  plane  may  be  written  as 


pit)  =  Pi(r)cos(wit  +  fair))  +  P2(r)cos(a>2t  +  fair))  (4) 

where  r  —  y/x2  +  y 2  is  the  radial  distance.  The  amplitude  functions  are  (Morse  and  Ingard 
1968,  Kino  1987) 


(5) 

-M£)} 

(6) 

where  Uw  is  the  velocity  amplitude  at  the  ith  transducer  element  surface,  and  A.,  =  2n/o>i  for 
i  =  1, 2  are  the  ultrasound  wavelengths.  The  phase  functions 

Mr)  =  -  (7) 

AjZo 


for  i  =  1,  2,  are  conveniently  set  to  be  zero  at  the  origin.  Also,  jinc(X)  =  Ji(2nX)/7tX , 
where  J\  (.)  is  the  first-order  Bessel  function  of  the  first  kind. 

The  instantaneous  energy  is  E  =  p2(t)/pc2.  Replacing  p(t)  from  equation  (4),  this 
energy  will  have  a  time-independent  component,  a  component  at  the  difference  frequency 
Ao)  —  <o2  —  a)\  which  result  from  the  cross  product  of  the  two  pressure  fields,  and  high- 
frequency  components  at  a)\  and  o> 2  and  their  harmonics.  The  energy  component  at  the 
difference  frequency  is 


^A<y(0  — 


P\(r0)P2(ro)  cos(Aqj;  +  AVr(r0)) 

pc 1 


(8) 


where  A^(r0)  =  Mro)  -  Mro)- 

Now,  we  define  a  unit  point  target  with  an  area  of  dx  dy  at  position  (*o,  yo)  on  the  focal 
plane,  and  with  a  drag  coefficient  dr(x 0,  yo)  such  that:  dr(x 0,  yo)djrdy  —  1  on  the  target 
and  zero  elsewhere.  This  equation  is  merely  used  as  a  mathematical  model.  In  this  case,  the 
projected  area  can  be  considered  to  be  S  —  d*  dy.  Therefore,  if  the  projected  area  is  unity, 
then  dr(x ,  y)  =  1,  which  according  to  equation  (2)  corresponds  to  a  totally  absorptive  object. 

Referring  to  equation  (1),  and  replacing  drS  with  unity  and  (E)  by  e&a »(0  of  equation  (8), 
then  we  can  write  the  low-frequency  component  of  the  radiation  force  on  the  unit  point  target 
as 

/a*>(*o,  yo;  0  =  -\p\ (ro)F2(r0)  cos(A cot  +  A ir(r0))  (9) 

pcL 

where  arguments  xo  and  yo  are  added  to  denote  the  position  of  the  point  target,  and 
r0  =  po  +  yjj.  Referring  to  equations  (5)  and  (6),  the  complex  amplitude  of  the  stress 
field  can  be  found  as: 


,  tra}  ..  (  roa\  \  f  na]  ..  (  r§a2  \ 

■h 


x  exp 


2  AXzq 


) 


(10) 
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where  AX  =  2ncj Aco  is  the  wavelength  associated  with  A co.  The  above  equation  indicates 
that  the  radiation  stress  is  concentrated  at  the  focal  point  and  decays  down  quickly  with  distance 
as  jinc(.)2. 

As  an  example,  we  consider  a  confocal  transducer  with  dimensions  ax  —  14.8  mm, 
a'2  =  16.8  mm,  a2  =  22.5  mm,  and  focal  length  70  mm.  Also,  we  assume  that  the  centre 
frequency  is  3  MHz,  and  the  difference  frequency  is  7.3  kHz.  The  radiation  stress  at  the  focal 
plane  of  this  transducer  is  plotted  in  figure  3. 

In  medical  applications  the  maximum  ultrasonic  intensity  is  regulated  for  safety  reasons. 
It  is  useful,  therefore,  to  write  the  stress  field  in  terms  of  the  peak  ultrasonic  intensity  at  the 
focal  point.  The  long-term  average  of  ultrasonic  intensity  at  the  focal  point  can  be  written  as 


/(0)  = 


P2(  0)  +  P22(  0) 
2  pc 


(ID 


where  Pi(0)  and  P2(0)  can  be  found  from  equations  (5)  and  (6)  as 

/(0)  =  y  [t/02i (j^)  +  (12) 

Assuming  Uoi  =  C/02  —  C/o,  we  can  write  the  focal  plane  stress  field  in  terms  of  1(0): 

r  <  s  27  (°) _ fg?AlgO _ ,_J  M 

&G>Xo,yo  c  (7ro2/A.,Zo)2  +  (jt/A.2Zo)2(«2  -  «?)zJ  V^-iZo 


The  first  fraction  in  the  above  equation  represents  the  static  radiation  force  produced  by  the 
two  beams  on  a  total  absorber,  the  second  fraction  is  a  constant  factor,  and  the  rest  represents 
the  spatial  distribution  of  the  stress  field  on  the  focal  plant.  If  Ao>  then  we  may  replace 
Ai  and  k2  by  Xo,  and  simplify  the  expression  for  the  stress  field.  Under  these  conditions,  the 
stress  at  the  focal  point  is 


FAw(0,  0) 


2/(0)  a}(al  -  a?) 
c  a\  +  (al-a%)r 


(14) 


The  fraction  on  the  right  represents  the  effect  of  transducer  dimension  on  the  stress  field.  For  the 
transducer  used  in  the  previous  example,  the  value  of  this  fraction  is  0.4999.  (This  fraction  can 
be  also  calculated  for  a  signal-element  transducer  of  the  same  diameter,  excited  by  an  amplitude 
modulated  signal.  For  this  purpose,  we  may  let  =  a2  and  a2  —  0.  For  these  values,  the  last 
fraction  in  equation  (14)  is  0.5.)  The  resulting  radiation  stress  (assuming  c  =  1500  m  s"1  for 
water)  is  FAw( 0, 0)  =  6.67  x  10“4/(0)  N  nr2.  If  we  let  7(0)  =  7200  W  nr2,  which  is  the 
maximum  spatial  peak,  time  average  intensity  suggested  by  the  FDA  for  in  vivo  applications, 
(FDA  1997),  the  resulting  stress  at  the  focal  point  is  FAaJ( 0, 0)  =  4.80  N  m“2.  Referring  to 
figure  3,  we  note  that  this  stress  field  is  applied  only  in  a  small  region  around  the  focal  point 
to  the  object. 


2.3.4.  Acoustic  emission.  To  explain  the  acoustic  emission  we  consider  an  ‘object’  within 
a  homogeneous  infinite  medium.  This  model  allows  us  to  separate  the  roles  played  by  the 
parameters  of  the  object  and  the  surrounding  medium,  and  can  be  fitted  to  various  applications. 
When  an  oscillating  stress  field  is  applied  to  the  object,  the  object  vibrates  at  the  frequency 
of  the  stress  field.  Vibrational  energy  of  the  object  is  partly  transferred  to  the  surrounding 
medium,  which  is  called  the  acoustic  emission  field.  Here,  we  calculate  the  acoustic  emission 
of  an  object  resulting  from  applying  the  radiation  stress. 
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Figure  3.  Three-dimensional  plot  of  the  stress  field  amplitude  of  a  confocal  annular  array  transducer 
at  the  focal  plane  at  a  point  in  time.  Transducer  dimensions  are:  a\  —  14.8  mm,  ct’2  —  16.8  mm, 
«2  =  22.5  mm,  and  the  focal  length  is  70  mm.  The  centre  frequency  is  3  MHz  and  the  difference 
frequency  is  7.3  kHz. 


To  explain  the  physics  in  more  generalized  terms,  for  the  moment  we  consider  an  object 
in  the  form  of  a  flat  disc  facing  the  beam.  Here,  we  assume  that  the  vibrating  object  has  a 
circular  cross  section  of  radius  b  and  uniformly  vibrates  back  and  forth  like  a  piston.  This 
choice  allows  us  to  illustrate  the  concept  in  a  simple  form.  We  also  consider  an  area  S  ^  nb2 
of  the  piston  surface  to  be  projected  normally  by  the  beam.  We  can  always  return  to  our 
elementary  point  object  by  reducing  the  area  of  this  disc  to  dx  d y.  Similar  solutions  can  be 
carried  out  for  objects  of  other  forms.  The  theory  can  be  also  extended  to  include  arbitrary 
vibrating-part  shapes  and  non-uniform  displacement  of  the  object.  The  total  radiation  force 
on  this  object  FA(U,  can  be  found  by  integrating  the  radiation  stress  over  the  area  of  the  object. 
This  force  vibrates  the  target  object  at  frequency  A co.  The  steady-state  normal  velocity 
amplitude  of  a  piston  at  frequency  A coy  f/Aa),  due  to  a  harmonic  force  FAw,  can  be  written 


=  ^  (15) 

where  ZAtt>  is  the  mechanical  impedance  of  the  object  at  Ao>.  The  mechanical  impedance  of 
the  object  has  two  components,  one  resulting  from  the  inertia,  friction  and  elasticity  of  the 
object  itself,  and  the  other  resulting  from  the  loading  effect  of  the  surrounding  medium  on  the 
vibrating  object.  The  mechanical  impedance  can  be  interpreted  as  a  measure  of  object  rigidity 
and  how  much  it  yields  to  the  applied  force.  For  example,  for  a  rigid  object ,  Zaw  is  high,  and 
hence  the  object  resist  the  force. 
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Knowing  we  can  calculate  the  pressure  field  it  produces  in  the  medium.  We  assume 
that  the  acoustic  emission  signal  propagates  in  a  free  and  homogeneous  medium.  The  farfield 
acoustic  pressure  due  to  a  piston  source  of  radius  b  set  in  a  planar  boundary  of  infinite  extent 
is  given  by  (Morse  and  Ingard  1968) 


P*a>  =  -jA  cop 


exp(jAa>//c) 
4 nl 


/ 2J\ [(A cob/c)  sin &]  cos  #  \ 

\  (A cob/c)  sin  #  cos  ft  +  Pb  ) 


(2 7tb2UAJ 


(16) 


where  /  is  the  distance  from  the  observation  point  to  the  centre  of  the  piston,  &  is  the  angle 
between  this  line  and  the  piston  axis  and  Pb  is  the  specific  acoustic  admittance  of  the  boundary 
surface.  (The  specific  acoustic  admittance  is  Pb  =  Pc/zb >  where  Zb,  the  acoustic  impedance 
of  the  boundary,  represents  the  ratio  between  the  pressure  and  normal  fluid  velocity  at  a  point 
of  the  object.)  The  factor  of  two  in  the  parentheses  comes  from  the  presence  of  the  boundary 
wall.  It  would  be  replaced  by  unity  if  the  boundary  wall  were  not  present  (Morse  and  Ingard 
1968). 

The  acoustic  emission  field  resulting  from  object  vibration  can  be  written  in  terms  of 
object  mechanical  impedance  by  combining  equations  (15)  and  (16)  as 


Pcuo  =  pc 2 


A(o  exp(j  Atol/c) 

2J\[(A(ob/c)  sin  #]  cos  &  \1 

j  c 2  4#/ 

^  (A cob/c)  sin  #  cos  &  +  pB  )  J 

2: rb2— — 
F 


(17) 


For  wavelengths  which  are  long  compared  with  the  object  size,  i.e.  when  bAco/c  — ►  0,  the 
term  in  the  large  round  brackets  approaches  a  constant  because  2  J(X)/X  — >  1  for  X  ->  0. 
Hence,  ignoring  the  first  fraction  in  the  large  parentheses,  we  may  consider  the  contents  of  the 
large  square  brackets  to  be  an  object-independent  function  (the  specific  acoustic  admittance 
Pb  relates  to  the  surrounding  boundary  surface).  Under  these  conditions,  the  large  square 
brackets  in  the  above  equation  represents  the  effect  of  the  medium  on  the  acoustic  emission 
field,  which  we  may  call  the  medium  transfer  function ,  and  denote  it  by 


„  ,fx  ,  A(o  expQAcol/c)  costf  /10, 

HAw(v  -  j  5  A  I  Q  Q  (to) 

c 2  4  nl  cos  #  +  Pb 

The  last  fraction  in  equation  (17)  includes  l/ZAft>,  which  is  the  mechanical  admittance  of 
the  object  at  the  frequency  of  the  acoustic  emission  (Aa>),  and  we  denote  this  by  Y^. 
It  is  convenient  to  combine  this  term  with  the  next  term  (2nb2)  in  equation  (17),  as 
Gao>  =  2nb2Y&(t)  =  2nb2/Z&(t),  which  is  the  total  acoustic  outflow  by  the  object  per  unit 
force  (acoustic  outflow  is  the  volume  of  the  medium,  e.g.  the  fluid,  in  front  of  the  object 
surface  that  is  displaced  per  unit  time  due  to  object  motion).  Function  Q^w  represents  the 
object  characteristics  at  the  acoustic  frequency.  We  may  thus  rewrite  equation  (17)  in  a  more 
compact  form  as 

PAa>  =  Pc2HAco(1)QaJfZ.  (19) 


Equation  (19)  indicates  that  the  acoustic  emission  pressure  is  proportional  to: 

(a)  The  radiation  force,  which  itself  is  proportional  to  the  square  of  ultrasound  pressure  and 
the  ultrasound  characteristics  of  the  object,  dr ,  and  the  projected  area  5. 

(b)  The  acoustic  outflow  by  this  object,  Q  aw,  representing  the  object  size  b  and  its  mechanical 
admittance  at  the  acoustic  frequency,  Y 

(c)  The  transfer  function  of  the  medium  at  the  acoustic  frequency, 

Note  the  difference  between  the  projection  area  S  and  the  vibrating  area  nb2.  The 
projection  area  determines  the  extent  of  the  force  applied  to  the  object  (equation  (1)).  The 
vibrating  area,  however,  influences  the  total  acoustic  outflow  in  the  medium  caused  by  object 
vibration. 
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23,5.  Measuring  the  acoustic  emission.  A  hydrophone  is  used  to  measured  the  acoustic 
emission  field.  Equation  (18)  indicates  that  for  wavelengths  which  are  long  compared  with 
the  object  size,  the  acoustic  emission  field  pattern  resembles  that  of  a  dipole  with  an  angle 
dependency  represented  by  the  last  fraction  to  the  right  of  this  equation.  This  angle  dependency 
is  normally  weak,  and  hence  the  position  of  the  receiving  hydrophone  is  not  critical.  Also, 
because  the  attenuation  of  low-frequency  acoustic  waves  in  tissue  and  liquids  is  normally  low, 
one  may  place  the  hydrophone  at  a  distance  from  the  object  without  losing  too  much  intensity. 
The  simplified  diagram  presented  in  figure  2  describes  how  to  measure  the  amplitude  of  the 
acoustic  emission  field.  To  measure  the  phase  of  the  acoustic  emission,  a  phase  reference  needs 
to  be  set.  One  way  to  accomplish  this  is  to  generate  a  reference  signal  by  electronic  down¬ 
mixing  of  the  two  CW  excitation  signals  to  produce  a  beat  signal  at  the  difference  frequency 
A/.  Then,  using  this  reference  signal,  the  phase  of  the  acoustic  emission  signal  (hydrophone 
output)  can  be  calculated  by  conventional  phase  detection  methods. 


2.3.6.  Applications.  Evaluation  of  the  characteristics  of  an  object  (or  a  medium)  by  listening 
to  its  sound  is  a  traditional  approach  that  has  been  used  for  many  purposes.  Qualitative 
evaluation  of  a  crystal  glass  by  tapping  on  it  is  a  simple  example.  Vibro-acoustography 
implements  the  same  approach  but  on  a  micro-scale,  and  in  a  way  that  could  be  applied  to 
tissues. 

Equation  (19)  represents  a  general  relationship  between  the  mechanical  parameters  of  the 
object,  surrounding  medium,  and  the  acoustic  emission  field  resulting  from  object  vibration. 
By  measuring  the  acoustic  emission  field,  it  would  be  possible,  in  principle,  to  estimate  some 
of  the  object  or  medium  parameters,  either  in  an  absolute  or  in  a  relative  sense.  For  example, 
one  can  use  vibro-acoustography  to  measure  the  resonance  frequency  of  an  object,  and  from 
that  information  it  is  possible  to  estimate  some  viscoelastic  parameters  of  the  object  or  the 
medium.  This  method  has  been  used  to  measure  the  Young’s  modulus  of  a  metallic  rod 
(Fatemi  and  Greenleaf  1999b).  In  another  experiment,  by  measuring  the  resonance  frequency 
of  a  known  resonator  in  a  liquid  medium,  the  viscosity  of  the  fluid  has  been  estimated  with 
good  accuracy  (Fatemi  and  Greenleaf  1999c).  These  applications  are  not  necessarily  medical, 
but  the  principle  could  be  used  for  evaluation  of  soft  tissues,  blood,  bone,  etc. 

In  medical  applications,  one  can  use  vibro-acoustography  to  obtain  an  image  of  the  object 
for  diagnostic  purposes.  In  such  applications,  the  image  may  not  represent  a  physical  quantity, 
rather  it  provides  a  means  to  visualize  object  details. 


2.3. 7.  Experimental  results.  In  the  following  experiment  vibro-acoustography  has  been  used 
to  image  excised  human  carotid  arteries  that  include  calcifications.  Characteristics  of  the 
ultrasound  transducer  used  here  are  those  described  in  the  example  presented  in  section  2.3.3. 
The  object  consists  of  two  excised  human  carotid  arteries  secured  on  a  thin  sheet  of  latex.  This 
sheet  is  almost  transparent  to  the  ultrasound  and  x-ray.  Figure  4(a)  shows  an  x-ray  image  of 
these  arteries.  The  left  artery  is  from  a  young  person  and  shows  no  calcification.  The  right  artery 
is  from  an  older  person  and  includes  some  calcium  deposit  which  can  be  seen  as  a  bright  region 
near  the  bifurcation.  The  arteries  were  immersed  in  a  water  tank  for  scanning.  The  difference 
frequency  A/  in  this  experiment  is  7  kHz.  The  lead  number  ‘2’  was  placed  on  the  latex  sheet 
for  identification  purposes.  Figure  4(b)  shows  the  vibro-acoustography  image  (displaying  the 
amplitude  of  the  acoustic  emission  signal)  of  the  object.  In  this  image,  the  calcified  region  is 
seen  as  a  bright  region,  resembling  its  x-ray  image.  Calcification,  which  is  formed  from  stiff 
(compared  with  the  arterial  wall)  material,  is  an  efficient  acoustic  radiator,  producing  strong 
acoustic  emission  when  exposed  to  the  radiation  force  of  the  incident  ultrasound.  Hence,  it 
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Figure  4.  Vibro-acoustography  of  excised  human  carotid  arteries:  (a)  x-ray  image  showing  a 
normal  (left)  and  a  calcified  (right)  artery  placed  on  a  sheet  of  latex;  ( b )  vibro-acoustography  of  the 
arteries  in  ( a )  measurement  at  a  vibration  frequency  of  7  kHz.  The  calcification  is  seen  as  a  bright 
region  in  both  x-ray  and  vibro-acoustography  images.  Lead  number  ‘2’  is  used  for  identification. 


stands  out  in  the  vibro-acoustography  image.  This  image  has  a  high  spatial  resolution  (about 
0.7  mm),  low  speckle  and  high  signal-to-noise  ratio.  The  number  *2\  which  is  also  made  from 
a  hard  material,  is  seen  with  high  contrast.  This  experiment  demonstrates  the  utility  of  the 
method  in  delineating  stiff  material  against  a  soft  background. 


3.  Comparisons  of  the  methods 

In  this  paper  we  have  presented  three  approaches  for  using  a  localized  ultrasound  radiation 
force  to  probe  biological  materials.  Here  we  compare  the  features  and  capabilities  of  the 
methods. 


3.7.  Objectives 

The  transient  method  estimates  a  quantity  that  correlates  to  the  hardness  of  tissue.  While  it 
does  not  estimate  a  specific  physical  parameter,  the  estimated  parameter  is  quantitative  and 
may  be  useful  for  diagnosis. 

Shear-wave  methods  aim  specifically  at  the  shear-wave  characteristics  in  soft  tissue,  in 
particular,  estimating  the  shear-wave  velocity.  Shear-wave  velocity  in  soft  tissue  is  very  low 
(typically  in  the  order  of  a  few  m  s“!),  allowing  localized  measurement  of  the  parameter.  Also, 
because  shear-wave  attenuation  in  soft  tissue  is  high,  the  boundaries  of  the  object  have  minimal 
effect  on  the  measurement.  This  method  is  not  intended  for  use  in  hard  materials. 

Vibro-acoustography  has  a  two-fold  objective.  In  this  method,  the  acoustic  emission  signal 
is  dependent  on  both  microstructure  (local  drag  coefficient)  and  macrostructure  (the  mechanical 
admittance)  of  the  object.  The  acoustic  emission  signal  can  be  analysed  to  estimate  different 
characteristics  of  the  object.  In  some  applications,  both  the  local  parameters  and  those  that  are 
related  to  the  large-scale  structure  of  an  object  are  needed.  An  instructive  example  presented  in 
Fatemi  and  Greenleaf  (1999a)  demonstrates  that  vibro-acoustography  can  be  used  to  produce 
an  image  of  a  tuning  fork  displaying  its  detail  (microstructure)  while  measuring  its  resonance 
frequency  (a  parameter  related  to  the  macrostructure)  represented  by  colour. 
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3.2.  Excitation  method 

In  the  first  two  methods  the  radiation  force  is  produced  by  a  single  focused  transducer  which 
is  driven  by  an  amplitude  modulated  signal,  in  particular  pulse  shaped  (the  shear-wave  method 
can,  in  principle,  operate  using  a  single  sinusoidally  modulated  beam).  Vibro-acoustography, 
however,  uses  two  CW  excitation  signals  (or  long  tone  bursts  of  CW  signals).  The  major 
difference  in  these  excitation  techniques  is  the  distribution  of  the  radiation  stress  in  the  depth 
direction.  A  possible  drawback  of  using  a  single  amplitude  modulated  beam  method  is 
that  the  pulsed  ultrasound  exerts  a  radiation  stress  on  all  tissue  through  which  it  propagates. 
Also,  a  shock  can  be  produced  on  the  transducer  surface,  and  then  transmitted  to  the  tissue. 
Such  phenomena  result  in  extraneous  motion  of  the  object  that  could  introduce  errors  in 
the  measurements.  In  vibro-acoustography,  the  two  CW  beams  pass  through  separate  paths. 
Outside  of  their  interaction  at  the  focal  region,  they  can  only  produce  a  static  radiation  force 
on  the  object,  which  produces  no  sound.  The  beams  mix  near  the  focal  point  to  produce  a 
modulated  field  that  can  generate  an  oscillating  radiation  force  only  in  that  location.  The 
resulting  acoustic  field,  therefore,  is  directly  related  to  the  excitation  of  the  object  at  the  focal 
area.  This  is  an  important  feature,  because  the  acoustic  emission  from  other  parts  of  the  object 
or  the  transducer  should  not  be  allowed  to  interfere  with  the  acoustic  emission  from  the  focal 
area. 


3.3.  Detection  method 

The  transient  method  and  the  shear  wave  method  are  based  on  localized  detection  of 
displacement  of  material.  This  is  implemented  by  using  Doppler  ultrasound,  optical 
displacement  measurement  or  a  more  sophisticated  method  such  as  MRI.  Vibro-acoustography 
uses  an  audio  hydrophone  or  microphone  for  detection,  which  is  a  much  simpler  detection 
method  than  Doppler  or  MRI. 


3.4.  Sensitivity 

Detection  sensitivity  is  critical  to  the  success  of  any  of  the  above  methods.  A  greater  ultrasound 
intensity  would  be  needed  if  the  detection  sensitivity  is  poor.  Ultrasound  Doppler  methods,  at 
conventional  ultrasound  frequencies,  can  probably  detect  displacements  of  the  order  of  a  few 
micrometres.  The  sensitivity  of  the  MRI  technique  is  of  the  order  of  100  nm  (Muthupillai  et  al 
1995).  Vibro-acoustography  has  been  shown  to  detect  motions  as  small  as  a  few  nanometres 
(Fatemi  and  Greenleaf  1998).  The  high  sensitivity  of  this  method  is  a  result  of  the  fact  that 
small  motions  of  the  object  can  produce  an  acoustic  emission  pressure  field  that  is  easily 
detectable  by  a  sensitive  hydrophone. 

3.5.  Measured  quantities 

The  quantity  measured  in  the  transient  method  is  somewhat  arbitrary,  and  not  a  direct  measure 
of  any  physical  quantity.  The  merit  of  this  method  is  that  tissue  hardness  is  measured  from 
relative  displacement  values,  without  the  need  to  measure  the  applied  force.  The  shear  wave 
methods  directly  measure  a  physical  parameter  of  the  material  in  an  absolute  sense.  The 
acoustic  emission  field  measured  in  vibro-acoustography  is  a  function  of  several  physical 
parameters  that  relate  to  the  object  and  the  surrounding  medium.  It  is  possible  to  measure 
a  single  parameter  only  if  one  has  enough  knowledge  about  other  parameters  involved.  For 
example,  one  can  measure  the  shear  viscosity  of  a  liquid  using  a  known  resonator  as  an  object 
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in  the  liquid  (Fatemi  and  Greenleaf  1999c),  or,  for  a  given  the  geometry,  the  Young’s  modulus 
of  a  metallic  object  can  be  measured  (Fatemi  and  Greenleaf  1999b)  by  vibro-acoustography. 

3.6.  Exposure  safety 

The  shear-wave  method  presented  in  Sarvazyan  et  al  (1998)  uses  high-intensity  pulsed 
ultrasound  for  excitation  and  either  the  optical  or  MRI  methods  for  detection.  It  is  argued 
that,  although  the  intensity  is  high,  the  total  exposure  is  within  the  FDA  limits  (FDA  1997). 
These  limits  are  recommended  based  on  possible  thermal  and  mechanical  effects  of  ultrasound 
on  tissue.  However,  in  vivo  use  of  such  detection  methods  would  be  difficult.  Using  Doppler 
ultrasound  for  detection  is  less  sensitive,  and  hence  requires  a  high  peak  ultrasound  intensity  to 
produce  enough  displacement  of  the  tissue.  This  intensity  level  may  exceed  the  recommended 
spatial  peak  pulse  average  intensity  limit. 

Vibro-acoustography  uses  either  CW  or  tone  burst  ultrasound.  Therefore,  the  continuous 
ultrasound  intensity  (spatial  peak  time  average  intensity)  limit  must  be  observed.  This  limit, 
according  to  the  FDA,  is  720  mW  cm-2  (FDA  1997).  However,  because  of  the  high  sensitivity 
of  the  hydrophone  detector,  it  is  possible  to  detect  very  small  levels  of  acoustic  emission  while 
using  low  transmit  power.  For  example,  it  has  been  shown  that  a  submillimetre  object  can  be 
detected  at  ultrasound  intensities  far  below  the  FDA  limit  (Fatemi  and  Greenleaf  1999a). 

3.7.  In  vivo  applicability 

In  vivo  application  of  the  transient  method  is  not  likely  because  it  requires  high  ultrasound  power 
to  produce  a  detectable  displacement.  Application  of  the  shear-wave  method  in  the  human  body 
largely  depends  on  the  detection  method  used.  Optical  methods  are  probably  not  usable  in  vivo. 
MRI  detection  methods  are  generally  complex,  especially  when  the  operation  is  to  be  linked  to 
the  ultrasound  system  in  the  magnetic  field.  Ultrasound  Doppler  is  a  more  practical  detection 
method  for  shear-wave  detection;  however,  its  sensitivity  may  not  be  sufficient  for  ultrasound 
intensities  within  the  recommended  maximum  level  (Walker  1999).  Vibro-acoustography  uses 
a  hydrophone  for  detection,  which  is  simple  to  operate  in  clinical  settings.  It,  however,  requires 
an  acoustically  quiet  environment  for  proper  detection.  The  biological  noise  spectrum  seems 
to  be  below  1  kHz,  which  can  be  easily  filtered  out  if  the  operation  frequency  is  above  this 
limit  (Fatemi  and  Greenleaf  1999a). 

3.8.  Image  resolution 

Imaging  complex  objects  using  the  transient  or  the  shear-wave  methods  has  not  been  fully 
explored  in  the  literature.  Vibro-acoustography  can  be  used  for  high-quality  imaging.  The 
spatial  resolution  is  proportional  to  the  width  of  the  main  lobe  of  the  stress  field.  For  example, 
for  the  3  MHz  ultrasound  transducer  presented  in  figure  3,  the  spatial  resolution  is  about 
700  p.m. 

3.9.  Material  properties 

Both  the  transient  and  shear-wave  methods  are  based  on  local  deformation  of  soft  tissue  in 
the  beam  direction.  These  methods  work  best  if  the  material  under  test  supports  the  shear 
wave,  has  enough  attenuation  to  allow  the  build-up  of  enough  radiation  force,  and  is  compliant 
enough  to  allow  appreciable  displacement.  Application  of  these  methods  for  hard  material, 
such  as  bone  and  calcifications,  would  be  difficult  or  impractical  because  (a)  stiff  materials 
have  low  compliance,  hence,  their  displacement  in  response  to  the  force  is  relatively  small 
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and  difficult  to  detect  and  (b)  the  shear  wave  travels  at  high  speed  in  such  material.  In  vibro- 
acoustography  hard  materials  produce  structural  vibrations  which  are  often  stronger  than  those 
of  soft  tissue.  Such  vibrations  result  in  strong  and  easily  detectable  acoustic  emissions.  The 
bright  appearance  of  the  arterial  calcifications  in  the  experiment  presented  in  figure  4  is  a 
demonstration  of  this  phenomenon. 

In  contrast  to  the  other  two  methods,  vibro-acoustography  can  be  used  to  detect  particles 
in  materials  that  do  not  support  shear  waves,  for  example  detecting  gas  bubbles  in  liquids.  The 
acoustic  emission  resulting  from  particle  vibration  includes  a  compressional  wave  component 
that  can  travel  in  the  surrounding  medium  including  the  liquid.  Tissue  attenuation  for 
compressional  waves  at  low  frequencies  is  small,  hence  such  waves  can  be  detected  by  the 
hydrophone  from  a  distance.  Shear  waves  in  a  liquid  medium  decay  very  rapidly  and  are 
difficult  to  detect  by  most  methods. 

4.  Summary 

Recent  progress  on  techniques  for  evaluation  of  dynamic  viscoelastic  parameters  of  tissue, 
based  on  the  radiation  force  of  ultrasound,  has  prompted  us  to  view  these  techniques  as  an 
emerging  new  field.  We  recognize  the  capabilities  and  potentials  of  this  new  field  as  unique 
and  distinguishable  from  those  of  conventional  elasticity  or  ultrasound  imaging  techniques. 

Three  methods  for  probing  biological  tissue  using  the  dynamic  radiation  force  of 
ultrasound  were  presented.  The  transient  method  measures  a  minute  transient  tissue 
deformation  versus  time.  The  shear-wave  method  also  uses  a  transient  excitation,  but  measures 
the  resulting  shear-wave  amplitude  and  velocity.  The  last  method,  vibro-acoustography, 
measures  the  acoustic  field  resulting  from  the  vibration  of  an  object  at  a  specified  frequency. 
Capabilities  and  limitation  of  these  methods  for  probing  different  materials  and  their  potentials 
for  in  vivo  application  were  discussed. 
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Chapter  1 


Introduction 


1.1  General  background 


Predicting  and  understanding  the  behavior  of  material  when  they  are 
subject  to  mechanical  forces  is  the  basis  for  many  modem  engineering  practices, 
from  material  science  to  the  design  of  solid  propellant  rockets  [Timoshenko  et  al. 
1970].  In  the  early  development  of  the  field  of  mechanics,  it  often  involved 
applying  loads  to  the  system  until  the  system  failed,  and  then  the  failure 
mechanism  was  studied  and  used  to  infer  the  behavior  of  the  material  prior  to 
failure.  Then,  as  the  technology  advanced,  it  became  possible  to  study  the 
behavior  of  complex  material  systems  using  non-destructive  testing  procedures, 
e.g.  X-ray  analysis,  acoustic  behavior  and  photo-elastic  behavior.  The  results  of 
experimental  study  have  contributed  to  the  understanding  of  material  properties 
and  to  development  of  mathematical  models  that  help  predict  the  behavior  of 
more  complex  material  systems  [Cheney  1997]. 
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In  recent  years,  this  concept  has  been  introduced  into  medical  diagnosis 
and  produced  significant  amounts  of  research  interest.  The  primary  motivation 
for  imaging  tissue  mechanical  properties  stems  from  the  following  two 
important  facts. 


First,  it  is  well  known  that  changes  in  tissue  mechanical  properties  such  as 
the  elastic  modulus  can  be  sensitive  indicators  of  pathology.  Palpation,  the 
routine  physical  examination  process  used  by  physicians  to  distinguish  between 
normal  and  abnormal  tissues  is  a  method  for  qualitative  estimation  of  tissue 
elasticity.  The  clinical  success  of  palpation  as  a  diagnostic  tool  is  evident  from  the 
fact  that  the  American  Cancer  Society  recommends  breast  self-examination  as  a 
screening  tool  for  early  detection  of  breast  cancer  [Holleb  et  al.  1991].  In  some 
series,  30%  of  breast  cancers  are  detected  by  palpation  only  and  not  by 
mammography  [Chapuis  and  Hessler].  Changes  in  mechanical  properties  may  also 
indicate  the  physiologic  state  of  the  tissue,  e.g.,  relaxed  and  contracted  states  of 
muscle  [Levinson  et  al.  1995;  Duck  1990]. 


Second,  it  is  also  known  that  the  elastic  modulus  can  vary  by  several 
orders  of  magnitude  between  various  tissues.  This  is  in  marked  contrast  with 
parameters  provided  by  other  existing  imaging  modalities.  For  example,  the 
ability  of  conventional  B-mode  ultrasound  imaging  to  differentiate  various 
tissues  depends  principally  on  the  acoustic  impedance,  which  in  turn  depends 
upon  the  bulk  modulus  of  the  tissue  under  examination.  Figure  1.1  is  a  summery 
of  material  properties  for  various  biological  tissues  from  the  literature 
[Sarvazyan  et  al.  1998].  It  shows  that  for  soft  tissues,  the  range  of  variation  of 
bulk  modulus  is  very  small  compared  to  that  of  the  shear  modulus.  In  X-ray 
imaging,  the  variation  in  X-ray  attenuation  coefficients  of  soft  tissues  is  less  than 
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a  factor  of  five  and  in  MR  imaging,  the  variation  in  the  tissue  relaxation 
parameters  spans  barely  an  order  of  magnitude  [Muthupillai  1997]. 
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Figure  1.1  Summary  of  data  from  the  literature  concerning  the  variation 
of  the  shear  and  bulk  moduli  for  various  materials  and  body  tissues 
After  Sarvazyan  et  al.  (1998). 

For  the  two  reasons  described  above,  imaging  techniques  that  exploit  the 
shear  modulus,  or  stiffness,  of  an  object  have  great  potential  in  medical 
application.  This  also  explains  why  elasticity  imaging  [Gao  et  al.  1996;  Ophir  et  al. 
1996,  1999],  a  general  material  characterization  method  under  extensive 
investigation  in  recent  years,  has  been  primarily  aimed  for  tissue 
characterization.  While  substantial  variations  exist  in  how  the  details  are 
implemented,  all  elasticity  imaging  methods  share  the  following  basic  elements: 

•  a  mechanical  force  is  applied, 

•  material  motions  are  measured, 

•  some  derived-quantity  is  reported. 

Elasticity  imaging  methods  can  be  classified  according  to  how  these  three 
elements  are  implemented.  Mechanical  force  can  be  applied  in  a  quasi-static,  or 
dynamic  fashion.  In  addition,  forces  can  be  either  endogenous  (naturally 
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occurring)  or  exogenous  (artificially  applied).  Exogenous  methods  can  be  further 
distinguished  based  on  whether  the  force  is  applied  at  the  surface  of  the  material 
or  deep  within  the  object  using  radiation  force.  Material  motions  can  be 
measured  using  acoustic,  ultrasonic,  optical,  or  nuclear  magnetic  resonance 
methods.  Furthermore,  there  are  several  quantities  that  can  be  usefully  derived 
and  reported  from  the  measured  motion  including  vibration  amplitude  and 
phase,  strain,  modulus,  and  apparent  phase-velocity  (wave-speed).  Figure  1.2  is  a 
summary  of  this  classification  scheme. 

Most  elasticity  imaging  methods  use  distribution  of  motion  within  the 
object  as  display.  However,  it  is  only  a  semi-quantitative  mapping  of  the  true 
mechanical  properties  of  the  underlying  material.  When  the  object  being  imaged 
is  very  complex,  the  interpretation  of  the  strain  distribution  image  is  not 
straightforward  and  various  artifacts  can  result  [Ophir  et  al  1996].  Various  efforts 
were  made  to  reconstruct  the  true  elastic  modulus  of  the  object,  as  will  be 
described  in  Chapter  2.  But  this  generally  requires  knowledge  of  stress 
distribution  within  the  object,  which  is  difficult  to  obtain.  Some  methods  (i.e. 
MRE)  cleverly  use  estimation  of  local  shear  wave  speed  to  solve  for  local  shear 
modulus.  Since  the  shear  wave  speed  is  also  related  to  shear  viscosity,  stiffness 
obtained  with  such  methods  is  not  truly  elastic  modulus. 
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Figure  2.2  Various  mechanisms  used  in  elastography  for  (a)  application 
of  force,  (b)  measurement  of  motion,  and  (c)  display. 
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Vibro-acoustography  [Fatem  and  Greenleaf,  1999]  is  a  new  elasticity 
imaging  method,  which  images  the  acoustic  response  of  a  material  to  localized 
harmonic  excitation  produced  by  radiation  force  of  two  interfering  ultrasound 
beams.  The  details  about  this  promising  imaging  method  and  its  beneficial 
features  will  be  described  in  Chapter  2.  However,  the  acoustic  emission  at  each 
scanning  point  is  actually  a  combined  response  of  the  reflectivity,  stiffness,  and 
object  emittance  at  that  particular  location.  Therefore,  the  vibro-acoustic  image  is 
not  a  direct  presentation  of  a  single  mechanical  property,  but  only  a  semi- 
quantitative  mapping  of  the  "stiffness"  distribution  of  the  object.  As  explained 
above,  truly  quantitative  techniques  are  more  preferable. 

In  this  thesis  work,  we  want  to  explore  the  possibility  of  using  vibro- 
acoustography  for  quantitative  characterization.  The  hypothesis  of  this  thesis  is 
that  it  is  possible  to  derive  mechanical  properties,  such  as  the  complex  shear 
modulus,  of  an  object  by  evaluating  the  response  of  the  object  to  many  excitation 
frequencies.  For  this  purpose,  a  quantitative  model  is  proposed  for  a  spherical 
target  embedded  in  a  viscoelastic  homogenous  medium.  The  goal  is  to  estimate 
the  mechanical  properties  of  the  medium  by  evaluating  the  displacement- 
frequency  response  of  an  embedded  spherical  target. 


1.2  Summary  of  chapters 

A  brief  summary  of  the  contents  of  each  of  the  following  chapters  is  given 

below. 


Chapter  2  consists  of  a  conceptual  introduction  to  mechanical  property 
measurement  and  a  review  of  the  currently  available  elasticity  imaging  methods. 
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Then  vibro-acoustography  is  introduced  as  a  new  method  to  imaging  the 
mechanical  properties  of  the  object.  The  advantages  and  disadvantages  of  this 
imaging  method  are  also  analyzed.  Finally  the  motivation  and  specific  aims  of 
this  thesis  work  are  presented. 

Chapter  3  presents  the  theoretical  model  of  a  sphere  in  water.  The 
dynamic  radiation  force  on  the  sphere  due  to  interference  of  two  ultrasound 
beams  is  derived.  This  driving  force  is  used  with  the  radiation  impedance  and 
the  mechanical  impedance  of  the  sphere  to  calculate  the  vibrating  velocity  of  the 
sphere,  from  which  the  acoustic  emission  by  the  sphere  is  also  solved.  Results  of 
experiments  are  reported,  where  the  velocity  and  acoustic  emission  of  several 
spheres  are  measured  versus  vibration  frequency.  The  results  agree  with  the 
prediction  of  theory. 

Chapter  4  extends  the  model  of  Chapter  3  to  viscoelastic  medium.  The 
formulae  of  radiation  force  and  radiation  impedance  are  modified  to  consider  the 
effects  of  shear  wave  and  attenuation.  Theory  shows  that  the  resonance  of  the 
sphere's  vibration  magnitude  versus  oscillating  frequency  is  related  to  the  shear 
properties  of  the  medium  around  the  sphere.  Experimental  measurements  on 
spheres  in  gel  phantoms  confirm  the  theory.  The  velocity  profiles  are  used  to 
estimate  the  shear  modulus  of  the  gel,  which  are  also  verified  by  independent 
measurement  methods. 

Chapter  5  is  the  verification  of  the  theory  in  viscous  fluid.  The  velocity 
profile  of  the  sphere  shows  no  resonance,  just  like  the  case  in  water.  The  slope  of 
the  curve  is  related  to  the  viscosity  of  the  fluid.  Experiments  are  performed  in 
viscous  standard  fluids  with  known  viscosity.  The  shear  viscosity  estimation 
derived  from  theory  agrees  with  the  certificated  value  of  the  test  fluid.  This 


method  is  also  used  to  measure  the  viscosity  changes  of  glycerin  solutions  versus 
temperature. 

Chapter  6  demonstrates  the  feasibility  of  measuring  vibration  velocity  of  a 
target  by  the  ultrasound  Doppler  method.  In  the  previous  chapters,  the 
vibrations  of  spheres  are  measured  by  a  laser  vibrometer,  which  is  limited  to 
applications  in  transparent  medium.  Results  shown  in  this  chapter  eliminate  this 
restriction.  An  analytical  expression  of  the  echo  from  a  harmonically  vibrating 
target  is  provided,  from  which  a  demodulation  scheme  is  proposed. 
Experimental  results  in  gel  phantom  are  compared  to  measurements  from  the 
laser  vibrometer,  which  is  used  as  a  gold  standard.  A  number  of  possible 
improvements  upon  this  ultrasound  Doppler  method  are  also  discussed. 

Chapter  7  summarizes  the  thesis  work  presented  in  this  dissertation  and 
provides  a  discussion  of  a  few  related  issues.  A  future  perspective  of  this  work  is 
also  presented. 
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